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PREFACE 

This book has been written at the suggestion of some students* 
who have been unable to find a book on the subject suitable 
for their requirements, and also with a view to completing 
the series of text-books on Mechanics which I have produced 
in recent years. 

The book is only an introduction to a subject which, in the 
eighteenth and nineteenth centuries, occupied the interest 
and attention of many eminent mathematicians. Nowadays 
the theory of potential is usually studied in connection with 
electricity and magnetism, which is given priority as a subject 
of study because of its wider developments and importance 
in relation to everyday affairs, so that there is a tendency to 
overlook the theory of the attraction of gravitating solids or 
regard it as a subject to be 'crowded out'. But the subject 
is of importance in connection with dynamical astronomy 
and the Figure of the Earth, and Todhunter's history of the 
subject 'from Newton to Laplace' ran to nearly a thousand 
pages. 

In the early chapters I have not forgotten the needs of 
students reading for a pass degree, though the book as a whole 
is intended for an honours course. A good many of the 
theorems in pure mathematics specially required for the 
subject are collected in the first chap^e?,- but^sti^dents who 
prefer to do so may begin to read at Chapter II and only refer 
back to Chapter I as the need arises. The later chapters are 
on Green's theorem, Harmonic functions and the Attraction 
of ellipsoids. 

I have indicated the source of the examples when taken 
from examination papers, the abbreviations M. T., C., and P. 
denoting the Mathematical Tripos, a College or Intercollegiate 
Examination, and the Preliminary Examination in Mathe 1 
matics in the University of Cambridge, and I gratefully 
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acknowledge the kind permission of the Senate of the Uni- 
versity of Ldhdon to use questions from the examinations of 
that University, the source of which is also indicated. 

A line in each set of examples indicates a rough division, 
of those which are easier from those which are more difficult} 
of solution. 

The proofs have been read by my friend Dr S. Verblunsky 
and again I tender to him my thanks. The book has been 
greatly improved by his valuable criticism and useful sugges- 
tions. I also wish to thank the staff of the University Press 
for their careful composition and assistance in eliminating 
errors. 

A. S. K. 

CAMBBIDGB 
May 1940 
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Chapter I 
PRELIMINARY MATHEMATICS 

1-1. We propose in this chapter to give a brief account of 
some mathematical ideas and propositions which the reader 
will find of use in later chapters. 

1*2. Surface and volume integrals . * In the theory of attrac- 
tions it is necessary to make considerable use of surface and 
volume integrals, and the applications of the theory frequently 
involve the evaluation of such integrals. The reader is therefore 
advised tq acquire some familiarity with methods of integra- 
tion from a textbook on the subject. We propose here merely to 
explain what is implied when such symbols as 

/ (x, y y z) dS and / (x, y, z) dv 

are used to denote integration over a surface or through 
a volume, and then to prove some important propositions 
connecting surface and volume integrals. 

C b 
A definite integral of afunction of one variable, say f(x)dx, 

J a 

may be defined thus: let the interval from a to 6 on the cc-axis 
be divided into any number of sub-intervals S l9 S 2 , ... 8 n , and 
let f r denote the value off(x) at some point on S r . Then, pro- 

n 

vided that the limit as n - oo of 2 / r S r exists and is indepen- 

r=l^ 

dent of the method of division into sub-intervals and of the 
choice of the point on 8 r at which the value off(x) is taken, 
this limit is the definite integral off(x) from a to 6. It can be 
shewn that the conditions for the existence of the limit are 
satisfied if/(#) is a continuous function. 

In the same way we may define \f(x,y,z)dS over a given 
surface; let the given surface be divided into any number of 
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small parts 8 X , 8 2 , . . . 8 n and let f r dfenote the value of/ (#, y, z) 

* n 

at some point on 8 r , then the limit as n -+ oo of f r S r , provided 

rl 

the limit exists under the same conditions as aforesaid, is 
defined to be the integral /(#, y, z) dS over the given surface. 

Any difficulty as to the precise meaning to be attributed 
to 'area of a curved surface' may be avoided thus: after 
choosing the point on each sub -division 8 r of the surface at 
which the value off(x,y,z) is taken, project this element of 
surface on to the tangent plane at the chosen point, and take 
the plane projection of the element as the measure of S r in 
forming the sum. 

The integral \f(x,y,z)dv through a given volume may be 

defined in a similar way. 

The evaluation of a surface integral of course involves 
a double integration, and that of a volume integral a triple 
integration, and when it is desired to exJiibit the forms chosen 
for the elements of area or elements of volume such symbols 

as /(#, y, z)dxdy&ud /(#, y, z)d#efo/dz are used; but when 
there is no doubt as to what is implied a single integration 
symbol is employed as above for the representation of an 
integral along a line or over a surface or throughout a volume. 

1"3. Green's Theorem. Ifu, v, w are functions ofx, y, z which 
have continuous derivatives with respect to x, y, z respectively 
throughout a singly -connected region bounded by a closed surface 
8 free from singularities and I, m, n denote the direction cosines 
of the normal drawn outwards from a point on an element dS 
of the surface, then 



where the surface integral is taken over the boundary and the 
volume integral throughout the region enclosed. 
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where the volume integrals are taken throughout the region 
considered and the surface integrals over its boundary. 

1*4. Solid angles. The solid angle of a cone is measured 
by th# area intercepted by the cone on the surface of a sphere of 
unit radius having its centre at the vertex of the cone. 

The solid angle subtended at a point by a surface of any form 
is measured by the solid angle of the cone whose vertex is at 
the given point and whose base is the given surface. 




Let PP' be a small element of area dS which subtends a 
solid angle da) at 0. 

Let the normal to dS make an acute angle y with OP, and 
let OP = r. Then the cross-section at P of the cone which dS 
subtends at is dS cos y, and this cross-section and the small 
area da> intercepted on the unit sphere are similar figures, 
so that 



Whence da> = (dS cos y )/r 2 j 

or d$ = r 2 sec da) ) 



It follows that the area of a finite surface can be represented 
as an integral over a spherical surface, thus 



|V 2 secyda> 



with suitable limits of integration. 
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1*41. </o> in polar co-ordinates, da) is an element of the 
surface of a unit sphere. Let the 
element be PQJRS bounded by 
meridians and small circles, where 
the angular co-ordinates of P are 
0, <f). Then since the arc PS sub- 
tends an angle d$ at the centre of 
a circle of radius sin 8, therefore 
PS = sin 0d</> ; and PQ = dO, so that 

= sm9d9d<f>. 




1*42. Volumes in terms of solid angles. 

of small solid angle dcu and 
length It. Let PP', QQ' be 
cross-sections at distances r, 
r + dr from the vertex 0. 

Then the area PP' = r 2 rfco, 
and, neglecting an infini- 
tesimal of higher order, the 
volume of the frustum PQQ'P' 
is r 2 da}dr. 



Consider a cone 




Hence the volume of the cone = c 



> ( R r*dr 
Jo 



Hence if be an origin inside a region bounded by a con- 
tinuous surface such that lines from each meet the surface 
in one point only and we divide up the region into narrow 
cones with vertices at and R denotes the radius vector or 
length of the cone of solid angle do>, then 



will represent the volume of the region. 

1*5. Scalar functions of position and their gradients. 

Let V (x, y, z) be a single-valued continuous function of the 
position of a point in some region of space. Suppose that the 
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function F is not constant throughout any region, so that 

the equation 

F (#, y, z) = const. 

represents a surface. We assume that through each point of 
the region in which F is 
defined, there passes a sur- 
face F = const. We also as- 
sume that at every point P 
on this surface there is a 
definite normal PN and that 
the tangent plane at P varies continuously with the position 
of P on the surface. 

It follows from the definition of F that two surfaces 

V(x,y,z) = a and V(x,y,z) = b 

cannot intersect; for if they had a common point it would be 
a point at which F had more than one value, in contradiction 
to the hypothesis that F is a single-valued function. 
Consider two neighbouring surfaces 

Let P, P' be points on each and let the normal at P to 
F = a meet F = a -f da in N r For small values of da, PN will 
also be a normal to F = a -f da. 

Then using V P to denote the value of F at P, we h&ye 

Fp,-F P da V N -V P V N -V P PN 



PP' ~ PP'" PP' ~ PN 'PP' 

. _ f*f\Q H 

- PN cos "> 

where 6 is the angle NPP'. 

Now if PP' = ds and PN = dn, and we make da and therefore 
also ds and dn tend to zero, the limit of (V P > V P )/PP f is 
the rate of increase of F in the direction ds and is denoted by 

9F 

-~ ; and similarly the limit of (1^ V P )/PN is the rate of 

OS 
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increase of F in the normal direction dn and is denoted by 

SV ' 

^ , and we have 



(1). 



Thus we have proved that the space rate of increase of V 
in any direction ds is the component in that direction of its 
space rate of increase in the direction normal to the surface 
V const.; or that if we construct a vector of magnitude dV/dn 
in direction PN, then the component of this vector in any 
direction is the space rate of increase of F in that direction. 

This vector is called the gradient of F and written grad F. 

To recapitulate: F is ^continuous scalar function of position 
having a definite single value at each point of a certain region 
of space, and the gradient of F is defined in this way: through 
any point P in the region there passes a surface F = const., 
then a vector, normal to this surface at P, whose magnitude 
is the space rate of increase of F in this normal direction, is 
defined to be the gradient of V at P, and it has the property 
**iat its component in any direction gives the space rate of 
increase of F iif that direction. It is clear that the gradient 
measures the greatest rate of increase of F at a point. 

1-6. Vectorial methods. We shall denote vectors by letters 
in Clarendon type. The magnitude of P is denoted by | P |. 
We recall that vectors are compounded by the parallelogram 
law and that they obey the commutative and associative laws 
of addition, viz. that 

p+Q=Q + p, and that P + (Q + R)==(P+Q)4-R. 

The product of a vector P and a scalar m is a vector whose 
magnitude is the product of the magnitude of P and m, i.e. 
m | P | , whose direction is that of P or the opposite according 
as m is positive or negative. Multiplication by scalars follows 
the ordinary rules of algebra, viz. 

m (nP) n (mP) = (mn) P; 



and 
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A vector s is called a.unit vector when | s | = 1. Any vector 
whose direction is that of a unit vector s can be represented 
by ms, where m is a scalar. 

For convenience in representing vectors by rectangular 
components we introduce three fundamental unit vectors 
i, j, k whose directions are mutually perpendicular; these 
directions being those of the co-ordinate axes Ox, Oy, Oz. 
Then, if P x , P y , P z be the projections of the vector P on the 
axes, we have a vector equation 

............... (1). 



1-61. Scalar product. Let P, Q be two vectors and </> 
the angle between their positive directions, then the expression 

| P | | Q | cos< 

is called the scalar product of the vectors and is denoted by 
(PQ) or (P, Q), or simply PQ. 

The scalar product of two vectors at right angles is zero, 
so that for the fundamental vectors i, j, k we have 

ij = jk = ki = and ii = jj = kk=l. 
Hence, from 1-6 (1), 



and on multiplying out this gives 



1-62. Again if V denotes a function of x, y, z, it follows from 
1'5 that the projections on the axes of the gradient of V are 

dV/dx, dV/dy, dV/dz. Hence we have the relation 

.dV .dV . dV 
= i^---f]^--f k^-, 

ox oy oz 

and the equivalence of the operators grad and i^--fj^--fk 

when the latter is applied to a scalar function. This latter 
operator is known as the Hamiltonian operator* and denoted 
by V. When applied to a scalar function it turns it into a vector. 

* After Sir William Rowan Hamilton (1805-1865), Astronomer 
Royal for Ireland, greatest of Irish Mathematicians, inventor of 
Quaternions. 
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We shall see in 1*65 that when V is applied to a vector it turns 
it into a scalar. 

The reader who wishes for further information on the 
subject of vector algebra may find it in the author's Dynamics, 
Part II, Appendix. 

^ll'63. A vector field. If to every point of a given region 
there corresponds a definite vector A, which may vary in 
magnitude and direction from point to point, then the region 
is called a vector field, or the field of the vector A; e.g. 
gravitational field, electric field. 

vl'64. Flux of a vector. If a surface 8 be drawn in the 
field of a vector A and A n denotes the component of A normal 

to an element dS of the surface, then the integral L4 n d$ is 

called the flux of A through S. Since a surface has two sides 
the sense of the normal must be taken into account, and the 
sign of the flux is changed when the sense of the normal is 
changed. The flux of a vector through a surface is clearly 
ct scalar magnitude, in fact if we use the symbol dS to represent 
a vector of magnitude dS directed along the normal to the 
surface, then the scalar product of the vectors A and dS 
denoted by (AdS) is the same thing as the subject of integra- 
tion A n dS, so that the flux of the vector may also be repre- 
sented by I (AdS). 

\Xl-65. Divergence of a vector field. Let A denote a 
vector field which has no discontinuities throughout a given 
region. Let &v denote any small element of volume containing 

a point P in the region and let L4 n d$ denote the outward 
flux of A through the boundary of Sv, then the limit as 8v -> 
of I \A n dS I / 8 v is defined to be the divergence of A at the point 

P and denoted by div A. 

Let P be the point (x, y, z). Let A x , A y , A z denote rectan- 
gular components of the vector A. 
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Let Z, m, n be direction cosines of the normal to dS. 
Then A n = I A x + mA 



so that L4 n d$= \(lA x + mA y +nA z ) dS 



by Green's theorem (1-3), where the volume integral is taken 
through 8v. The value of the integral is the product of 8v and 
the mean value of the integrand, say the value at a point Q. 
As Sv->0, this mean value tends to the value at P. 



TT T A T J 9^-r <^4;/ <^4* 

Hence divA= hm ^ = -^4----^ + --^ (1) 

S->o <w d# oy dz 

and we note that this result is independent of the form of the 
element of volume 8v which surrounds the point P, and also 
independent of the way in which the element contracts to 
zero. 

Further, in the notation of 1-61, 1-62, we have 

.) 



dx dy dz 
= divA, 

so that when the operation V is performed on a vector A the 
result is the scalar function div A. 

Also if A = gradF (1*5), where V is a scalar function of 
*,y,*,sothat 97 W d_V 

x) y> z x > y> z > 

we have div grad V = div A 



dx dy dz 
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9 2 3 2 3 2 
The operator ^-2 + ^-2 + ^2 * s usu * a Uy denoted by V 2 , for 

it is the same as 

9 .3 



The reader will note that the operator V may be applied either to 
a scalar, turning it into a vector, or to a vector, turning it into a scalar; 
i.e. that V means 'grad' when applied to a scalar and it means 'div' 
when applied to a vector. But the operator 'grad ' only has a meaning 
when applied to a scalar function, and ' div ' only has a meaning when 
applied to a vector. Also the order of the operations is not inter- 
changeable, thus 



but grad div F is meaningless. 

Again, div grad A is meaningless, but 



166. If we assume that the divergence as defined in 1-65 is indepen- 
dent of the shape of Sv, it is easy to obtain the result ( 1 ) by taking for 
Sv a rectangular parallelepiped with edges parallel to the axes and of 
lengths Sic, y, 8z with P at its centre. The method is given in full in 
the author's Electricity and Magnetism, pp. 6, 7. We shall not reproduce 
it here, as in the next article we shall apply the same method to a 
more general case. 

1*67. Curvilinear co-ordinates. Gradient and Diver- 
gence. Let the equations 

jfi(a?,y,z) = a, f 2 (x,y,z) = P, f 3 (x,y,z) = y ...... (1), 

in which a, j3, y denote variable parameters, represent three 
families of orthogonal surfaces; i.e. surfaces which cut at 
right angles wherever they intersect one another. 

We assume that one and only one surface of each family 
passes through any point (x,y,z), so that the parameters 
a, j8, y m^y be regarded as curvilinear co-ordinates of the 
point whose cartesian 'co-ordinates are x, y, z. 

Let ds lt ds 2 , ds 3 denote elements of arc of the curves of 
intersection of the pairs of surfaces j3, y; y, a and a, j3. Then 
since a is the only parameter that changes along the curve 
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t 

j8, y, it follows that ds 1 must be a multiple of da. We may 
therefore suppose that 



where h l9 A 2 , h 3 are in general functions of the co-ordinates, 

Since the surfaces a, /?, y through any point are at right 
angles, a vector A at the point may be resolved into three 
components normal to these surfaces and denoted by 

J A n A 

"> --p> -^y 

Consider an element of volume dv bounded by surfaces 
a ^da, j3 |dj8, y |dy. Let the surface a cut this element of 
volume in the curvilinear rectangle 
PQRS whose centre is the point 
(a, )8, y). 

To find an expression for div A in 
curvilinear co-ordinates we have to 
find the flux of A out of dv. The 
component of A at normal to 
PQRS is A^. Hence, if we assume 
what can easily be provedf, that the 
average value of a function over a 
small plane area is the value at the centroid, then the flux 
across PQRS is A^ (area 




* If the values of h lf h z , h s are required, they may be found by 
differentiating each of equations (1) with regard to a, ft and y, and 
solving the resulting equations for dx/dot., dx/dfi, etc. Then 



v 

t If/ (x, y) be the function, S the area and (x, y) its centroid, we have 



. + higher powers of A?% 



since / x'dS = / y'dS = 0, 

where c is of order of the greatest linear dimension of S. 
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Now the parallel surfaces ABCD^EFGH are the surfaces 
a |da and a + ^da, and the outward flux of A across A BCD is 



where e x is an infinitesimal which -> as da. -> 0. 
Similarly the outward flux across the face EFGH is 



where e 2 has a like meaning to l . 

Therefore the contribution of this pair of faces to the flux 
of A out of dv is 



^ (h^A^) + e 



where -> as a -> 0. 

In like manner for the other pairs of faces of dv, so that 
the total outward flux of A is 

{^\ r\ ^\ 

g- (A a h 3 A J + ^ (Ag A!^) + g- (Aj Aa^y) + ^ 

where 17 -> as dv -> 0. 



snce vA= m 

dv -+ 

therefore 



Hence, since divA= lim -= , and dv h<J 

dv -+ Q dv 



...... (3). 

Further, when A is the gradient of a scalar function V, 
we have 

_ 3 Z f!Z a E = !?Z liE L!Z 

' ^' y~' ' ' ' 



and by substituting (4) in (3), we have divgrad V or 



, , _ 

9y\ A 3 

...... (5). 
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1-68. Polar and cylindrical co-ordinates* When using 
polar co-ordinates r, 0, </>, the families of orthogonal surfaces 
are concentric spheres r = a, cones = j8, and planes <f> = y , and 
it is easy to see that 

ds l) ds 2 , ds 3 ~dr, rdd, rsint 
so that h l9 h 2 , h 3 == 1, r, rsin0, 

and 

1 



V 2 F = 



r 2 sin l3r 
or V 2 F = 4^~^ 2 







3 2 F 



r a sin 2 0a< 2 
(1). 






When using cylindrical co-ordinates r, 0, z, the families 
of orthogonal surfaces are co-axial cylinders r = a, planes = j8, 
and planes z = y, and in this case we have 

ds l9 ds 2 , ds 3 ==dr, rd9, dz, 
so that h l9 A 2> A 3 ~l,r, 1, 



and 



or 



V 2 F = -5- 
r or 



.(2). 



Either of the results (1) and (2) can of course be obtained 
directly by applying the method of 1'67 to the appropriately 
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shaped element of volume, and we suggest as an exercise to 
the student \he independent verification of results (1) and (2) 
in this way. 

!?. If a vector field A is continuous throughout a region E 
bounded by a surface S\ then the outward flux of A through the 
surface S is equal to the volume integral of div A throughout It. 

For we may subdivide R into a large number n of elements 
of volume St^ , 8v% , . . . 8v n . The total flux of A out of all these 
elements will be the outward flux across $; because over all 
boundaries which are not parts of S the flux will be taken 
twice in opposite directions. 

Hence if S r denotes the boundary surface of 8v r , we have 



where c r -> when 8v r - 0. The latter step follows from the 
definition of div A as a limit. 

Therefore, by increasing n so that the Sv's tend to zero, 
we get r 

\A n dS or (AdS)= divAefo ......... (1). 

1-71. The theorem of 1-7 may be extended to apply to the case 
in which the region R contains a surface S' across which the vector 
field is discontinuous, i.e. there is a sudden change in A in crossing 
the surface S'. 

Suppose that S f divides R into two parts R l and R z in each of which 
we assume that A is continuous and let A ln , A 2n denote the normal 
components of A in JR l , R% close to ', taken positively towards S' in 
both regions. Then we can apply the theorem of 17 to each of the 
regions R , R z separately and add the results. This gives 



= JA n d8+ I <A ln + A tn )dS' ......... (1). 

J o J & 



div Ado 

JK 

Across any surface at which A is continuous, we have, of course, 
^ n ^ A ln and the result reduces to 17 (1). 

* See the author's Electricity and Magnetism, pp. 8, 9. 
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1-72. Green's Second Theorem. The theorem of 1-31 may 
be obtained in vector form from 17, thus: 

Let A = FgradF', where F, V are scalar functions which 
with their first and second derivatives are finite and continuous 
through the region bounded by the surface S. 

Then f ( A dS) = IV (grad F'dS), 

and div A = VA = V (FVF') 



= grad F grad V + V div grad V. 
Therefore 1-7 (1) is equivalent to 

("grad F grad V dv = JV (grad F'dS) - IV div grad V'dv 

...... (1), 

or if we use the Hamiltonian operator 

...... (2); 



fvFVF'cfo = IV(VF'dS)- 



and the symmetry of the left-hand side shews that F and F' 
might be interchanged on the right-hand side. The reader 
will have no difficulty in identifying (1) and (2) with 1-31 (2). 

1-8. Convergence of volume integrals. In the theory of 
attraction and kindred subjects we frequently have to consider 

integrals of the form - , where dv denotes an element of 

volume and /> a scalar function of position, and the region 
of integration R is bounded by a surface 8 which includes 
the origin O from which the radius vector r to a point of dv 
is measured. It follows that the integrand becomes infinite 
within the range of integration, and we have to investigate 
the circumstances under which the integral can be said to 
have a meaning. 

Draw a closed surface $ x inside S and surrounding O, then 
so far as the subject of integration depends on r it is finite 
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between &^ and S, and 



has a definite value, where the 



limits indicate that the integral is 
taken through the volume bounded 
by S and S v If as S l contracts upon 
the point the last integral tends to 
a definite limit which is independent 
of the shape of S l9 then 



lim 



? pdv 




is defined to be the value of the given integral through R 
and the integral is said to converge. If the integral does not 
converge, it is said to diverge. 

A condition for convergence is that corresponding to any 
small positive number e there must be a surface 8 : , defined as 
above, such that for all surfaces S 2 inside S l and surrounding 



< 



.(I). 



For if a limit A exists we can choose S l such that 



e, and 



L^-L^ 



and therefore - 

Is. r' 1 



g fa pdv 

? JSi ^ n 



Conversely, if condition (1) holds, then a limit exists, for 
by taking $ x small enough we can make the fluctuation of 
the value of the integral for different cavities inside S 1 as 
small as we please and therefore small without limit. 

Now let M denote the greatest value of p within a finite 
region round 0. Let rj denote the greatest value of r on /S 1 
and r 2 the least value of r on $ 2 . Let S be the sphere of radius 
T) and a the sphere of radius r 2 . Then 



(2); 
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but, as in 1*42, we may put dv = r 2 da>dr, so that the contribu- 
tion to the last integral of the cone of solid angle do>, when 
n< 3, is 



da> ^ 
J r 



or 



Consequently < ~ < 

^ J Jar n 3-nJ 

pdv 



or 



47719 

' 



,3-71. 



3-n 



and 



,3-n 



Hence, provided 7i<3, for any arbitrarily chosen positive 
number e however small we have only to take 77 less than 

_JL 
{(3 n)/4 : 'jTM} 3 ~ n in order to get a surface S x such that for 

all surfaces $ 2 inside S L 



<. 



Consequently the given integral is convergent if n < 3. 

We note that the form of the outer boundary $ does not 
come into consideration, because it is only in the vicinity of 
the origin that the integrand tends to become infinite. 

It can be shewn that for n ^ 3 the integral is divergent.* 

* For further information on the subject see J. G. Loat hem's Tract 
on Volume and Surface Integrals used in Physics, Camb. Univ. Press, 
1905. 
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GRAVITATIONAL ATTRACTION AND 
POTENTIAL. SIMPLE APPLICATIONS 

%!. The Law of Gravitation. The law discovered by 
Newton* is as follows: Every particle in the universe attracts 
every other particle with a force which is directly proportional 
to the product of the masses of the particles and inversely propor- 
tional to the square of the distance between them.'f 

Thus if m, m' denote the masses of two particles and r their 

VYivn ' 

distance apart, the force of attraction between them is y 5- ; 

where y is known as the gravitation constant, and measures 
the attraction of two particles of unit mass at unit distance 
apart. 

It is evident that some limitation must be imposed upon the 
dimensions of the particles in order to avoid a difficulty in 
defining the distance between them. For this purpose we 
may define a material particle as a body so small that, for 
the purposes of our investigation, the distance between its 
different parts may be neglected. 

We shall prove later that a uniform solid sphere attracts 
external particles as though its mass were concentrated at 
its centre; and this may be regarded as some justification for 

* Sir Isaac Newton ( 1 642- 1727), Lucasian Professor of Mathematics 
at Cambridge and Fellow of Trinity College. His great work Philo- 
sophiae Naturalis Principia Mathematica was published in 1687 by the 
Royal Society, then under the Presidency of Samuel Pepys. He made 
discoveries in Optics and invented his Methodus Fluxionum which 
constituted the foundation of the Calculus. In later life he left 
Cambridge for London, becoming President of the Royal Society and 
Master of the Mint. 

f An outline of the process of reasoning by deduction, from Kepler's 
laws of planetary motion through the propositions of Newton's 
Principia to the law of gravitation, is given in the author's Dynamics, 
Part I, pp. 167-8. 
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the hypothesis we propose to make, namely, that for the 
purpose of calculating attractions we may regard particles 
as quantitjesj^i^ 



2-11. Numerical value of y. An approximate value for y 
may be determined thus: Regarding the earth as a sphere 
of radius r and assuming it to attract a small body on its 
surface as though its mass were collected at its centre, then if 
E denotes the mass of the earth in grammes and M that of 
the body, EM 

y = weight of the body = Mg ...... (1). 



So that, if p denotes the mean density of the earth, E = frrpr 3 , 
and from (1), 

y=^ ........................ (2). 

4npr 

But, in O.G.S. units, 



<7 = 981 cm. sec.~, p = -7 approx. 
and \7tr = a quadrant of the earth's circumference 
= 10 9 cm. 

rr-u r 3x981 1 

Therefore y = = appr X ....... 



when C.G.S. units are used. 

The mean density of the earth has been determined by a 
series of experiments in which the attraction of the earth on 
a given body has been compared with the attraction on the 
same body of another body of known mass.* 

2-12. Astronomical units. We can avoid the continuous 
repetition of the constant y throughout our work by choosing 
units such that y = 1. Such units are called astronomical or 
theoretical units. 

We observe that the acceleration /produced by the attrac- 
tion of a particle of mass m on a particle at a distance r is 

ryn 

given by f=y~^\ so that y= 1 when/, m and r are all unity. 
* See Encyd. Brit, article * Gravitation*. 
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Hence the astronomical unit of massiis the mass of a particle 
which by its attraction produces unit acceleration at unit 
distance. 

We can find the astronomical unit of mass in grammes by 
taking the above formula for acceleration, which holds good 
in all systems of units, and putting r = 1 cm. and/= 1 cm. sec."" 2 , 
so that ym 1, or w= l/y = 15,500,000 grammes. 

In what follows we shall omit the constant y, save that in 
the collections of examples from examination papers we shall 
print the questions in the form in which they were set, in 
order that readers may become familiar with the use of y. 

wl2*2. Field of force. When we speak of the attraction of 
a body or system of particles at a point external to itself, 
we mean the force of attraction which the body or system 
would exert on a particle of unit mass placed at the point. 
There must be a definite value (including zero at a point of 
equilibrium) for this force at every point at which a particle 
can be placed. Thus we arrive at the conception of a field of 
force, or region of space with every point of which there is 
associated a force which is definite in magnitude and direction. 
Inasmuch as a force belongs to the class of physical magni- 
tudes known as vectors, the field of attraction of any system 
<jf bodies is an example of a vector field. 

^2*21. Attraction of a system of particles. It is part of 
the Newtonian theory of attrac- 
tion that the force exerted at a 
point by a system of particles 
is the vector sum of the forces 
exerted by each of the particles 
separately. This is called the 
principle of superposition of 
fields of force. 

Let particles of masses m l9 
w 2 , w 3 , ... be situated at points 
A l9 A z , A%, ... whose co-or- 
dinates referred to rectangular axes are a l9 b l9 c x ; a 2 > ^2> C 2 
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#s> &a> C 8 ; .... It is required to find expressions for the com- 
ponents X , Y, Z of the attraction of the system &f particles 
at any point P whose co-ordinates are x> y, z. 

Let r 1? r 2 , r 3 , ... denote the distance P-4 1? P^4 2 > ^^s? 
The attraction at P of Wj^ at A l is w^ 2 directed along 

i A * T . a, a; 61 y c^ z . 

x; and PA l has direction cosines -- , - - , , where 



Therefore the particle m x exerts at P a force whose com- 
ponents parallel to the axes are 



The other particles ra 2 , m 3 , . . . make like contributions, so 
that, by the principle of superposition of fields of force, we 
have 



...... (1), 

where s=l, 2, 3, ... and th-.tiinmation extends to all the 
attracting particles. 

^2'22. Potential. Let us now define a function V by the 
formula 

F=S ........................ (1), 



where r, 2 = (o.-) a + (6 8 -y) 2 + (c a -) 8 ......... (2), 

in the notation of 2*21. 

Thus defined, V is a function related to a system of attracting 
particles having a definite value at every point P of space 
external to the particles.. It is a function of the co-ordinates 
(x, y, z) of P and is clearly a single-valued function, in the 
sense that it cannot have more than one value at each point 
P; for it represents simply the sum of the masses of the separate 
particles divided by their respective distances from P. Further, 
though we have expressed V as a function of x, y, z, it clearly 
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represents a sum which does not depend on the particular 
system of sixes of reference. 
Now, by differentiation, we have 

^I == _S^ s 
dx r 8 * dx ' 

dr 
But, from (2), r tf g= -(a s -x), 

. . dV 

so that ^- 

dx 

and by 2-21 (1) this is equal to X. Similarly 

-Y and -* ................ (3). 

3y dz 

The function F defined by (1) is called the potential of 
the attracting particles, or the potential of the field of 
force. We have proved that the derivatives of V with 
regard to x, y, z give the components of attraction at P in 
the directions of the axes. And since the directions of the 
axes can be chosen arbitrarily, it follows that the space- 
derivative of V in any direction gives the component of 
attraction in that direction. This may be verified thus: let 
9/9 denote a differentiation in a direction ds whose direction 

7 9# 9?/ dz ._ 

cosines are L m, n or ~ , ~ , -=- . Inen 

os os ds 

dV dVdx dVdy dVdz 
_ _ __ _ _________ i ___ &_ _| __ _ _ 

9s dx ds dy ds dz ds 

9F 9F 9F 

-~- + ^-5- 
dy dz 



= force component in direction ds ...... (4) . 

In the language of vectors the forcie of attraction is the 
gradient of the potential (1*5), or, if R denotes the resultant 
attraction at a point at which the potential's F, 

R = gradF .../. ............... (5). 



2-23] POTENTIAL IN TERMS OF WORK 25 

The potential function V plays an important part in the 
theory of attractions, since, as we have just seen, if the 
potential of any given distribution of matter can be determined 
the attraction at any point can be found immediately by 
differentiating the potential. 

"2-23. Physical interpretation of potential. With the 
notation of 2*22 we have, for the total differential of the 
potential, 



^- -^- -^- 

. dx dy y dz 

= Xdx+Ydy + Zdz, 
from 2-22 (3). 

Hence, by integrating along any path from P to Q 9 we get 

. ......... (1). 



But the integral represents the work which the forces of 
attraction would perform upon a particle of unit mass as it 
moved along this path from P to Q. This gives us a measure 
of potential in terms of work per unit mass, in that the potential 
at any point Q exceeds the potential at any other point P by the 
work which the forces of attraction would perform upon a particle 
of unit mass as it moves along any path from P to Q. 

It is clear that the addition of a constant to the potential 
will not affect the values of the force components since these 
are obtained by differentiating the potential. Also, when the 
potential is determined by integration, from known force 
components, as in (1) above, the constaftt of integration may 
be so chosen as to make the potential vanish at tn infinite 
distance from the attracting matter, since this will $txjord with 
the definition of 2'22 (1). On this hypothesis we see that the 
potential at a given point due to a given attracting system u 
the work that would be done by the attractions of the system on 
a particle ofufiit mass as it moves along any path from an infinite 
distance up to the point considered. / 

We have seen that the definition of potential as S/m/r) 
leads to this expression for potential in terms of wo^k done 
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per unit ma^s, and it is easy to demonstrate the converse, thus : 

Let m be the mass of a typical 

particle of the system situated 

at the point A. Let PP' = ds be 

an element of any path from an 

infinite distance to the point Q, 

and let AP = r, and AP' = r + dr. 

Then so far as the field of force 

depends upon the particle m at 

A its value at P is m/r 2 directed along PA; and the work 

done by this force on the unit particle as it moves from P to 

P' is -z I r I ds. or r . Hence the total work done by 

r 2 \ ds / r 2 J 

the attraction of the particle m on a unit particle moving 
from an infinite distance to the point Q is 




J. 



m 

or -r-pr. 
AQ 



By the principle of superposition of fields, the total work 
done by the attractions of all the particles of the system is 
obtained by adding for their separate effects, so that 



AQ 

gives the potential at Q] where the formula represents the 
sum of the masses of the separate particles each divided by 
its distance from Q. The interchangeability of the two defini- 
tions of potential is thus completely established. 

N.^2'24. Dimensions . Gravitational potential must not be con- 
fused with potential energy. They have different physical 
dimensions; because the dimensions of potential energy are 
those of work, i.e. ML 2 T~ 2 , in terms of the fundamental units 
of mass, space and time; but the dimensions of gravitational 
potential as explained in 2*23 are those of work per unit mass, 
i.e. L 2 T~ 2 , or the square of a velocity ^ /^ 

In this connection it is important to remember that we are 
using astronomical units (2*12) and omitting the gravitational 
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constant y, and thougii this does not affect thp argument 
when potential is defined as work per unit mass, yet when we 
use the formula for potential of 2*22 (1) if we want to find the 
dimensions of V we must restore the constant y and write 

V = yS (m/r), because y has dimensions. By definition y g- 

represents a force and is therefore of dimensions MLT~ 2 , so 
that the dimensions of y are M- 1 L 3 T~ 2 , and the dimensions 

777 

of y are L 2 T~ 2 as above. 

' r 

It is important also to note that though potential energy 
decreases when work is done, gravitational potential increases. 

2*25. The definition of potential in terms of work (223) leads at 
once to the relation 

R = gradF 

of 2-22. Thus if P, P' are any two neighbouring points at a distance da 
apart, F, V + dV denote the potentials at P, P' and J? the component 
of attraction at P in direction PP', we have 

(V + dV) F = increment in potential from P to P' 

= work done by attractions on unit mass passing from 
Pto P' 



dV 
Therefore F= 9 and this being true for components of attraction 

OS 

in every direction round P it follows that 



V2-26. Equipotential surfaces and Lines of force. Re- 

garding the potential F of a given attracting system as a 
function of co-ordinates x, y, z, the equation 

F (x, y, z) = const. 

represents a surface over which the potential is coiistant. 
Such surfaces are called equipotential surfaces. It is 
evident from the definition of potential that only one such 
surface passes through any point of space, so that no two 
equipotential surfaces can intersect. Also since F has a con- 
stant value over such a surface no work would be done by 
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the attractions on a particle moving on such a surface, there- 
fore at every point the resultant attraction is normal to the 
equipotential surface through the point. This is also obvious 
from the relation R = grad V (1*5). 

A line so drawn that its tangent at each point is in the 
direction of the resultant attraction at the point is called 
a line of force. It is evident that a line of force is at right 
angles to the equipotential surfaces at all their points of 
intersection; and, conversely, a surface which cuts all lines 
of force at right angles must be an equipotential surface, 
because at no point on the surface is there a component of 
force tangential to the surface, so that no work would be done 
on a particle moving on the surface and there could therefore 
be no variation in the potential. 

v'a-S. Continuous bodies. We have now to pass from the 
attraction components and potential of a system of separate 
particles to the attraction components and potential of 
distributions of matter regarded as continuous bodies. '"We 
remark in the first place that such bodies are not really 
continuous but consist of separate molecules and ultimately 
of indivisible electrons* By the principle of superposition 
therefore there should be no difficulty in obtaining the attrac- 
tion components and potential of such a body provided that 
we have a means of summing the contributions of the separate 
particles. It is natural to look to integration to effect the 
summation, but when we consider what the process of inte- 
gration involves, we find that it does not fit the physical 
conditions of the problem precisely, and it is only by giving 
a special interpretation to our conception of 'body' that we 
can justify the use of integration. Thus it is usual ^represent 

the potential of a continuous body by a volume integral 



pdv/r, where dv denotes an element of volume of the body at 

a distance r from an external point P and p denotes the density 
of matter in dv. The process of integration implies that the 
density of the body is continuous; but this is not the case 
for a real body with a molecular structure. To justify the 



242] LONG CIRCULAR CYLINDER 33 

24 COB. for the attraction cf a long rod, we have 
attraction at P due to rod through QQ' 

along PC, 



attraction at P due to rod through RR' 

= pp = 2m sec </> dO along PR. 

Hence these two rods exert equal and opposite attractions at P; 
and by dividing up the whole cylinder into similar pairs of rods we 
see that its resultant attraction at any internal point* is zero. 

It follows that the potential must be constant at all points inside 
the cylinder sufficiently far from its ends, and is therefore equal to 
the potential at 0, and by 2-41 (5) the potential due to the rod QQ' is 

2mQQ' log 21 - 2mQQ' log OQ, 

where 21 is the length of the cylinder. Therefore, for the whole potential, 
F = 2Mlog2Z-2Mloga ..................... (1), 

where M = 2-rrma, a being the) radius. 

Again, for an external point, take P' the inverse of P, so that 
OP . OP' = a 2 , and, by similar triangles, OP'Q = OQP = < = OP'R. Then 
the attraction at P 7 of the rod through QQ' 

2mQQ' 2mPQdOsee<t> Jn , a . ,-. 

= "" p^T = P r Q = sec * ' of' g Q ' 

But the resultant attraction of the cylinder at P' is clearly along 
P'O, and by resolving the attraction of the rod thrbugh QQ' in this 

direction we get ^.^ dd. Similarly, other rods give like results and the 

whole attraction at P' is therefore sypT = /-)/J and we observe that 

this is the same as if the whole mass of the cylinder were condensed 
into a rod of equal mass along the axis of the cylinder. 
To find the potential, we may put r for OP' and write 

, dV__ _2M 
dr ~~ rV 

so that V = const. 2-M^log r, 

where, in order that F may take the form (1) when r = a, (3-7), the 
constant is 2M log 2/, so that 

V = 2M log 21- 2M log r. 

It is clear that the same results will be true for the attraction and 
potential of a uniform long solid cylinder at external points, for the 
cylinder can be divided up into coaxial thin cylinders for each of which 
the results are true. 



34 ATTRACTION AND POTENTIAL [2'5- 

5. Attraction of a uniform circular disc at a point on 
its axis. Let be the centre of the disc, P a point on its 
axis Oz Sit a distance z from 0, and 
let ra denote the surface density or 
mass per unit area. Divide the disc 
into concentric rings. Let OQ = x be 
the radius and QQ' = dx the breadth 
of one of these rings. The mass of the 
ring is 2-jrmx dx, and the attraction at 
P of each element is got by dividing 
the mass of the element by PQ 2 . But 
the resultant attraction of the ring is clearly in the direction 
PO, so that its magnitude is 




where 6 is the angle OPQ. But x = z tan 6, so if a is the angle 
which a radius of the disc subtends at P, we have for the whole 
attraction of the disc 



a 277-m2 2 tan 6 sec 2 cos 6 

_ .,,..,., / 

2 ~2 5~ 

n ^ sec t/ 



cosoc) ...(1). 



For an infinite plate we may put a = |TT, so that the attraction 
of an infinite plate is 27rm* at right angles to itself. But we 
may not infer that the attraction of a circular disc at a point 
on its axis close to itself is 2irm unless the thickness 9&the 
plate is negligible, so that the particles of a ring in the 
immediate vicinity of the attracted particle may be regarded 
as equidistant from it. 

A frustum of a right circular cone bounded by planes at 
right angles to its axis may be regarded as a pile of discs all 
subtending the same angle at the vertex, and the attraction 
at the vertex of the cone due to such a frustum is easily seen 
to be 27Tw(l-~cosa), where t denotes the thickness of the 
frustum, m the mass of unit volume, and a the semi-angle 
of the cone. 
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a thick frustum is got by adding the attraction of a number of thin 
ones, it is sufficient to prove the theorem for parallel thin slices of 
a cone of small solid angle. 

Let PQ be a slice of the cone of small thickness t, and let p be the 
density, O the vertex, OQ^r and 
6 the angle the normal to PQ 
makes with OQ. Then the area 
PQ r z sec0da} (1-4), and the 
mass of the frustum = pr 2 t sec daj . * r - 

JPtfe" attraction -at O of this 

frustum is therefore ptsec6da>, 

and this has the same value for all parallel frusta of the same thickness. 

(ii) Prove that the potential of a solid of uniform density p at an external 

point P can be represented by a surface integral %p I cos 6dS over the surface 

of the solid, where 9 is the angle between the inward normal to dS and the 
line joining dS to P. 

Let a cone of small solid angle da> p 

and vertex P cut the surface of the solid 

in elements of area dS l9 dS 2 at A, B 9 ^, 

where the inward normals make angles 
e l , 2 with the line BAP. Let 




The mass of an element of volume of 
the cone at a distance r from P is, as in 
1*42, pr*da>dr. Hence the mass of the 
cone between A and B produces at P a 
potential , df 

I 

z cos 6 2 + dS l cos ^). 
If we take the sum for all cones which intersect the solid wo shall get 
I cos 6dS integrated over the surface. 
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EXAMPLES 



1. A sphere of 150 kilogrammes, placed with its centre J30 cm. 
vertically below that of another sphore, is found to cause an apparent 
increase of the latter's weight by the 1'14/10 8 part. What value does- 
this imply for the constant of gravitation ? [M. T. 1918] 

j 2," Three equal uniform rods, each of length 2a and of mass m pet 
unit length, form an equilateral triangle. Shew that the attraction 
on a particle of unit mass at a point on the straight line through the 
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centre of t^e triangle at right angles to its plane, and at a distance z 
from it, is A, 

1-, , ,,... [C- 1906] 



3. Prove that any plane closed rectilinear polygon of thin uniform 
attracting wire, all of whose sides touch one and the same circle, 
exerts no attraction at the centre of the circle. 

[London Univ. 1927] 

4. A uniform wire of mass M has the form of a semicircular arc 
of radius a together with the diameter of the semicircle. Find the 
attraction of the wire upon a particle of unit mass placed at the other 
end of the diameter through the middle point of the arc of the semi- 
circle. [London Univ. 1933] 



Prove that the attraction of a thin uniform cylindrical shell of 
radius a and length I at a point on its axis at a distance b from one 
end and I b from the other (b < %l) is 

M T __ 1 _____ 1 



where M is the mass of the shell. 

[The ends of the shell are open and circular.] 

[London Univ. 1926] 
**"* \ 
( 6* Shew that the attraction of a uniform solid hemisphere at the 

centre of its plane base is f yM/a 2 , where M is its mass and a is its radius. 

7. Prove that the potential of a circle of uniform surface density a 
and radius a at a point on its circumference is 4yaa. 

[London Univ. 1932] 

8. Prove that the potential of a plane annulus bounded by concentric 
circles at a point on its axis is 27ry/> (l f I), where p is the mass per unit 
area and I, V are the distances of the point from the inner and outer 
edges. 

9. A uniform circular lamina is divided into two segments by 
a chord which subtends an angle 2a ( < TT) at the centre. Prove that 
the attraction at the centre of the circle due to the minor segment is 

2yor [log tan ( JTT + |a) sin a], 
a being the surface density of the lamina. [London Univ. 1931] 

10. Prove that, if the surface density of a circular disc at any point 
distant r from the centre is Ar, then the attraction at a point on the 
axis at a distance d is 

a 2o_) 



Try Ad | log 



[London Univ. 1932] 
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11. Four uniform rods, each of the same mass per unit length, 
form a square A BCD. Any point P is taken inside the Square and 
E, F, G, H are the feet of the perpendiculars from P upon AB, EC, 
CD, DA respectively. Prove that the potential at P due to the square 
A BCD is the same as that due to four rods EF, FG, OH, HE each of 
the same mass per unit length as that of the rods forming the square. 

[London Univ. 1935] 

12. Two uniform rods A B and CD, of length a and mass m per unit 
length, are placed so that CD lies along a right bisector of AB, and C 
and D are distant a and 2a respectively from the rod AB. Show that 
the attraction of either rod on the other is of magnitude 

2ym 2 {log e ( 1 + V5) + loge 2 - log. ( 1 + V 17)}/a 2 . 

[London Univ. 1938] 

13. Find the attraction of a uniform circular ring, of mass M 
and radius a, upon a thin uniform rod, of mass M ' and length 21, 
lying along the axis of the ring with its middle point at a distance 
x ( < I) from the centre of the ring. 

If the ring is fixed and the rod released from rest when its centre is 
at a small distance x from the centre of the ring and moves under the 
mutual attraction only, shew that the period of small oscillation is 

27r(a 2 + P)*/(yAf )*. [London Univ. 1935] 

14. Find the attraction of a uniform thin shell, in the form of the 
portion of a circular cylinder between two planes perpendicular to 
the axis, on a particle at the axis; and prove that a particle slightly 
displaced on the axis from the position of equilibrium will perform small 
oscillations of period 2rrM-t (a 2 -f c 2 )*, where 2a is the length and c 

the radius of the cylinder, and M its mass. [C. 1892 

v 

15. Prove that the potential due to a thin rod AB, of uniform line 
density m, at any point P in the line AB produced is ymlog (PA/PB). 
Two such rods, each of length 2a, are placed so that they lie in the 
same straight line with their nearer ends at a distance 2a apart. 
Prove that a particle, released from rest at a point in the line of the 
rods at a small distance x from the point midway between them, begins 
to move with acceleration 16yrax/9a 2 , approximately. 

[London Univ. 1936 

_6. Prove that the resultant attraction at P of a straight rod AB 
of uniform line density m is 2yme/PG directed along PG, where e is 
the eccentricity of the ellipse of foci A , B passing through P, and G 
is the point in which the normal to the ellipse at P cuts AB. 

17. The ends of a thin uniform rod are at the foci of a hyperbola, 
and a similar rod is placed parallel to it with its ends on the curve. 
Shew that the attraction between the two rods is inversely proportional 
to their distance apart. [C. 1908] 
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18. Shew that the potential of a closed hemispherical shell, made 

of uniform Ihin material of surface density or, at a point on the axis 

of symmetry at a distance a from its plane face remote from its curved 

surface is 2-nyaa, where a is the radius of hemisphere. 

/ [London Univ. 1937] 

\S." f " 

w*19. AB, A'B' are two parallel rods of line densities m, ra' and are 
placed with the line joining their middle points perpendicular to the 
rods and of length c. Shew that the attraction between the rods is 

^^-(AB'-AA'). [C. 1914] 

,SQ. Two uniform rods AB, AC, each of mass m and length a, are 
freely jointed at A. Shew that the couple required to keep them at 
right angles against their mutual attraction is (2 \/ 2) ym?/a. 

[London Univ.] 

\S nr i 

^2L Two uniform thin rods AB and CD are of linear density A 
and A' and length 2a and 2c respectively. The midpoint E of AB and 
the midpoint F of CD are at a distance b apart, and AB, CD and EF 
are mutually perpendicular. Shew that the attraction between the 
rods is 



where y is the gravitational constant. [P. 1934] 

j f jti& Shew that the attraction of a uniform solid circular cylinder, 
of* radius a and density p, at an external point P on the axis of the 
cylinder is 



where a, ft are the angles subtended at P by the radii of the nearer and 
farther ends of the cylinder. [London Univ. 1925] 

23. A finite solid, of uniform density p, is bounded by the paraboloid 
x 2 + y 2 = 4az and the plane z = a. Prove that the intensity of attraction 
at the origin is 

27rypa[l - V 5 + 4 log. J ( 1 + Vs)]. 

[London Univ. 1938] 

24. A quantity of homogeneous matter of density p is in the form 
of a portion of a paraboloid of revolution bounded by a plane perpen- 
dicular to the axis. Prove that its attraction at its vertex tends, as c 
the length of the axis increases, to the limit 4-jrpa log c/ae, where 4a 
is the length of the latus rectum of the generating parabola. 

H 25. A uniform solid, of density p, is in the form of the finite portion 
of a paraboloid of revolution of latus rectum 4 cut off by a plane 
perpendicular to the axis at a distance 2a from the vertex. Prove that 
the attraction at the centre of the circular end is 

[London Univ. 1933] 
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26. Shew that for a long rod of length 21 and line density m, the 
component attractions parallel and perpendicular to the rod are 
approximately 

and 



where the origin is at the middle point of the rod and the axis of x 
along it, and terms of order equal to or higher than (d/J) 4 are neglected, 
where d is the distance from the middle point. [M. T. 1919] 



27. A homogeneous prism, infinite in length, whose cross-section 
is an equilateral triangle ABC attracts a particle at A; prove that the 
resultant attraction is 47rAf/3a, where M is the mass of a unit length 
of the prism and a is the length of a side of the triangle. [C. 1908] 

28. ABC is a normal cross-section of a uniform triangular prism 
of infinite length and of mass m per unit length. Show that the resultant 
attraction at the point A is compounded of 4ymA/a perpendicular 
to BC, and 4ymlog(6/c)/a parallel to BC, the angle A being measured 
in radians. [London Univ. 1925] 

20. Two infinitely long uniform circular cylinders, each of mass m 
per unit length, are placed with their axes parallel. Prove that the 
potential, F, at a point, outside both cylinders, whose distances from 
their axes are r, r', is given by 

V const. 2ym log (rr f ) . 

Further, prove that the lines of force outside the cylinders are arcs of 
rectangular hyperbolas. [London Univ. 1936] 

'30/ A tetrahedron is made from a tiniform sheet of thin metal. 
Shew that, if the law of attraction were the inverse cube, a particle 
would be in equilibrium if placed at the centre of the inscribed sphere. 

ftl* Two thin straight uniform rods AB, CD are pivoted together 
at their middle points. Prove that the attraction between them 
reduces to a couple of moment 2mm' (AC~ BC) cosec a, where m and 
m' arc their line densities and a is the angle between them. [C. 1903] 

32. The cross-section of a long uniform solid cylinder is a semicircle. 
Shew that the direction of the resultant attraction at a point on an 
edge far from the ends of the cylinder makes an angle tail" 1 (2/;r) 
with the plane face. [C. 1928] 

33. Shew that a uniform thin square lamina of mass M, whose 
side is of length I, exerts an attraction at a point P on the perpen- 
dicular to the lamina through one corner, whose component parallel 
to the lamina is 

V2 y M, r(fe4.p)*(l + (fe + p)t}-[ 
P 10g L h{l + (h* + 2l*)l} J> 
where h is the distance from P to the lamina. [Lpndon Univ. 1931] 
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34. A plane sheet of matter of uniform surface density is bounded 
by two infinite parallel straight edges and extends to infinity in both 
directions. Find the components of the attraction parallel to and at 
right angles to its plane; and shew that the loci of the points for which 
these are constant are two sets of orthogonal cylinders. [C. 1898] 

35. Shew that, if a uniform straight rod CD were free to slide in 
a smooth groove parallel to a fixed similar rod A B, its time of oscillation 
about its equilibrium position would be 



I CD \BD AD) } ' 

where m is the mass of unit length of tho rod AB 9 and AD, BD are 
distances in tho equilibrium position. [C. 1926] 

36. Determine the potential of a solid uniform right circular cone 
at the centre of its base. [London Univ. 1931] 

37. A vertical solid right circular cylinder of height h, radius a 
and density p, bounded by plane ends perpendicular to the axis, is 
divided by a plane through the axis into two parts. Shew that the 
horizontal attraction of either part at the centre of the base is 



38. Shew that two finite surfaces having the same surface density, 
which are geometrically inverse to one another, exert equal attractions 
at the centre of inversion. [C. 1906] 

39. If the law of attraction of a thin uniform straight rod is that of 
the inverse cube, shew that the equipotential surfaces are generated 
by the revolution of a family of curves, whoso polar equation is 

(r 4 + a 4 - 2r 2 a 2 cos 20)* = 2ar sin 6 . cosec (cr sin 9) 

about tho initial line, 2a being tho length of tho rod and c a parameter. 

[C. 1907] 

40. The mass per unit length at a point P of a rod A CB, of which A 
and B are the ends and C a particular point, varies as CP/\/(AP . PB). 
Prove that the resultant attraction of the rod at any point at which 
AC, CB subtend equal angles acts through C. [M. T. 1902] 

41 . An infinitely long thin rod, of uniform line density ra, and a thin 
uniform rod of mass M and finite length 2a are so placed that the 
shortest distance between them is of length b and bisects the finite 
rod. Prove that the resultant attraction of either rod on the other 
is equivalent to a wrench of intensity 

2yMra0/(a sin 0) 
and pitch b cot (tan <j> <f>)l<f>, 

where is the angle between the rods and <^ = tan~ 1 [(a sin 0)/6]. 

[London Univ. 1938] 
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42. A solid uniform hemisphere is of mass M and radjus a. Prove 
that its attraction at any point P of its axis on the convex side of the 
curved surface is equal to 

M/ R* + 2r*-WR\ 

7" 2 \ O? ) ' 

where r and R are the distances of P from the centre and rim of the 
base. [C. 1902] 

43. Shew that, if a denotes the radius and p the density of a hemi- 
spherical hill and R the radius and o- the mean density of the earth, 
the difference of the latitudes as measured by a plumb line of stations 
on the north and south sides of the base of the hill is increased by the 
attraction of tho hill by approximately ap/Rcr. 

44. Consider the effoct on the plumb lino of a deep crevasse running 
east and west, and shew that if the density of the earth's crust is 
greater than two -thirds of the earth's moan density the north edge 
of the crevasse will have a lower apparent latitude than the south edge. 

45. Prove that at the odgo of a very long straight depression of 
rectangular section, whose depth h is small compared with its breadth 
c, the plumb line will bo deflected through an angle approximately 
oqual to 2fjL?i(l+logc/h)/g from its mean position, where p, is tho 
attraction at unit distance of unit volume -of the superficial strata 
supposed condensed to a point. [M. T. 1887] 

^46. Tho surface density of a circular disc of radius a is a ( I r 2 /a 2 )* 
at a distance r from the centre of the disc. Prove that the potential 
at a point on tho axis of the disc at a distance a cot (f> from its centre is 

Try cr a ( tf) cosec 2 < cot <f>) . 

[London Univ. 1927] 

47. Prove that a lamina, of uniform density o- in tho form of the 
segment of a circle, attracts a point P on the remaining arc of the circle 
with a force tending towards D whoso magnitude is 



where AB is the base of tho segment, D tho middle point of its arc and 
a the radius of the circle. [C. 1906] 

48. Prove that a solid regular tetrahedron of uniform density p 
exerts at a vertex an attraction of intensity 

yph (3 cos" 1 J TT), 

where h is the perpendicular distance of a vertex from the opposite 
face. [London Univ. 1936] 
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49. Prove, that a thin triangular lamina ABC of uniform surface 
density a exerts, at a point on the line through A perpendicular to 
the plane of the triangle and at a distance k from A, an attraction the 
intensity of whose component perpendicular to the lamina is 

k cos B T k cos O 

- sm "- sm 



where p is the perpendicular distance of A from BC. 

[London Univ. 1939] 

50. If a uniform lamina is the portion of the plane xy containing the 
origin and bounded by the hyperbola x 2 /a 2 y*/b z = 1 , shew that 
the z component of its attraction at any point on the ellipse 

x z /(a 2 4- b 2 ) + z 2 /b 2 1 , y = is -77 m PA where h is the z co-ordinate 
of the point. [C. 1900] 
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attract with a force Mass/r 2 . Hence it is easy to see that the 
required attraction is 



J'51. Example. Find the law of density when the attraction at an 
internal point of a sphere of given mass is constant and equal to its value 
at the surface* 

Let a be the radius and p the mean density so that the mass is fyrrpa 9 . 
Then with the notation of 3-5 



so that 

Differentiating with regard to r, we get 
/(r)r 2 = J/>ar, 
and/(r) = o gives the law of density. 



^SE . Uniform shell bounded by non-intersecting spheres . 

When the boundary spheres are not concentric, let A, B be 

their centres and a, b their radii. Let p 

be the density of the matter between 

the spheres, and let r, r' be the distances 

of any point P from A and B. The field 

of force due to the given shell is the 

difference of the fields due to two ujii- 

form spheres of radii a, b and centres 

A, B. Hence, using the results of 3-31, 

the potential has the following values: 




(a 3 6 3 \ 
-- _ _ I j 

(26 3 \ 
3a 2 r 2 --- 7 ~ ) '> 



and in the hollow F 3 = f irp (3a 2 - r 2 - 36 2 + r' 2 ). 

Also, using the results of 33, the resultant attraction has 

components as follows: when P is outside the shell ^p-^ 



4-2 
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6 3 

along PA and f^-^ along JSP; in the substance ^npr along 






and f 770-7^ along JSP; and in the hollow &npr along P-4 
r * 

and f 7T/>r' along JSP. And in this last case the forces compound 
into a resultant force |?rp BA parallel to BA ; so that the force 
inside such a hollow shell is constant jn magnitude and 
direction, i.e. a uniform field of force 

3*6. Long uniform solid cylinder, we may regard a solid 
cylinder as composed of thin cylindrical shells having the same 
axis. We have seen in 2*42 that the attraction of a uniform 
long thin circular cylindrical shell is zero at internal points, 
and at external points is the same as if the mass were con- 
centrated uniformly along the axis of the cylinder. For the 
attraction at an internal point P, all the shells external to P 
exert no attraction at P, and the resultant attraction at P 
is therefore due to the shells nearer to the axis than P, so that 
if M be the mass of unit length of the cylinder and a its radius, 
and r the distance of P from the axis, the cylinder attracts 
at P like a rod of mass Mr 2 /a 2 per unit length along the axis; 
i.e. the attraction at P is 



For the potential we have 



_ 
dr ~ a 2 "' 

so that V = const. - M r 2 /a 2 , 

where r < a. But from 2-42 the potential at an external point 

is F = 21flog2Z-2Jflogf, 

where 21 is the length of the cylinder. And in order that these 
expressions may give the same value for V on the surface 
r a (3-7), the internal potential must be 

V = 2M log 21 - 2M log a + M (1 - r 2 /a 2 ). 

3*7. Variation in the attraction in crossing a surface on which 
there exists a thin layer of attracting matter. Let P x , P 2 be two 
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points on opposite sides of the surface close to one another 
on the same normal to the surface. Let 
m be the mass per unit area of a small 
circular disc of the surface between P t 
and P 2 . 

The points may be regarded as so 
close to the disc that, in accordance 
with the result of 2-5, it exerts equal attractions 27rw at P t 
and P 2 both directed towards the surface. 

The attraction of the rest of the matten on the surface (or 
elsewhere) may be represented in the neighbourhood of P t 
and P 2 by a normal component X taken in the sense PiP a 
and a component Y in some direction tangential to the surface. 
And the values of these components at P l , P 2 will differ by 
infinitesimals of the order PP%. For if X be the value at 

o -\r 

P x and P x P 2 = 8n, the value at P 2 is X + ^ Sn + higher powers, 

C/Ti 

and the infinitesimal may be neglected if we neglect the 
thickness of the layer of matter. 

We now compound the attraction of the small disc and the 
attraction of the rest of the matter, and let X l9 Y l and 
X%, Y 2 denote the components of the whole attraction at 
P l , P 2 , taking X l , X 2 in the sense P l P 2 , and we have 



Zirm, Y^Y 
and JT a = Z-2iiw, F 2 =F. 

Therefore X 2 - Z x = - 47rra, Y 2 =Y l ......... ,..(1), 

or the tangential component of attraction is unaltered and 
the normal component is diminished by 47rra in crossing the 
surface. 

As regards the potential, it is clear from 2-51 that a thin 
disc of matter produces the same potential at points on 
opposite sides close to itself, and the potentials of the rest 
of the matter at P x , P 2 only differ by an infinitesimal, which 
may be neglected if we neglect the thickness of the layer of 
matter. We conclude that the potential is continuous across 
such a surface as we are considering, in the sense that there 
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is no jumj in its numerical value. This is illustrated by the 
case of the thin spherical shell of 3*21. The internal and 
external potentials are represented by different functions, but 
they have the same value on the sphere r = a. 

Reverting now to the result (1), if we use V l9 V% to denote 
the potential functions on opposite sides of the surface, and 
dn to denote an element of the normal taken positively from 
P l to P 2 , we have, at the surface, 



and 



3V dV 

U'% W\ 

dn dn 



.(2). 



We remark that m denotes the mass per unit area or the 
'surface-density' of the layer, and the thickness of the layer 
is neglected; this implies that a finite mass is condensed into 
a volume of no thickness, so that the volume density of the 
matter is infinite. 



l. Attraction of a thin uniform spherical shell at a 
point of itself. The attraction at a point of itself of a thin 
layer of matter depends on the shape of the gap in the surface 
in which the attracted particle is placed. We may define the 
principal value of the attraction as the limiting value of the 
attraction at the centre of a circular 
hole when the radius tends to zero. 
We proceed to calculate this for a 
spherical shell. Let m be the mass 
per unit area of the shell, and omitting 
a small circular element of the shell 
surrounding a point P, consider the 
attraction of the rest of the shell at P. 
Let an element QQ' of area dS subtend 
a solid angle da) at P. The attraction of this element at P is 
mdS/PQ* along PQ, and resolving this along PO, which is 
clearly the direction of the resultant attraction, we get 
mdScosOPQ/PQ 2 , or mdScosOQP/PQ*, which is equal to 
mdco. If now we allow the gap in the shell round P to shrink 
to vanishing point, we see that we have to take the sum 
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for all cones on one side of the tangent plane at P, so 



that the resultant attraction is 

N3*8. Mutual work of two Attracting systems. Let there 
be a system of particles of masse&rai , m 2 , m 3 , . . . at the points 
A l , A 2 , A 3 , . . . and a second system of masses m/, m 2 ', w 3 ',. . . 
at the points B l , JB 2 , J5 3 , . . . . 

If we suppose that initially the particles of the second 
system are scattered at an infinite distance and are brought 
from thence up to their final positions B l , J3 2 , B 3 , ... by the 
attractions othe first system, the work done by these attrac- 
tions is represented by 



h ... or 
where V l9 V 2 , V 3) ... are the potentials of the first system at 

Similarly the work done by the attractions of the second 
system in bringing the particles of the first system from a state 
of infinite diffusion up to their final positions A l9 A z , A%, ...is 
represented by 

Vi'nii + V 2 'm 2 + V 3 'm 3 + ... or S V'm 9 
where F^, V 2 ' y Fj', ... are the potentials of the second system 

If we denote the distances A l B l , A 1 B 2 ^ ... A p B q , ... by 
r n , r 12 , ... r pq , ..., so that 

* _L. ^ 1 3 . 

^ T "1 r * j 



T7 " _i_ 2 _i_ 3 

V 2 " "T ' I 

^12 ^22 r 32 



it is evident that each of the expressions SFw', SF'm is equal 
to the sum Sram'/f extended to the product of masses of 
every pair of particles, one of each system, divided by their 
distance, so that 

SFm'^SF'm ..................... (1^ 
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For continuous masses we have in like manner 



(Vdm'= (V 



dm (2), 

where dm, dm' denote elements of mass of the systems and 
V, V the potentials of the first system at dm' and of the 
second system at dm. 

These expressions may be called the mutual work of the 
two systems, or the exhaustion of their potential energy. 
For potential energy can be considered to exist in matter in 
a state of infinite diffusion, and work is done and potential 
energy lost when the matter assumes a more condensed form, 
whether under the attractions of another system or under the 
attractions of its component parts. When we speak of the 
mutual potential energy of two attracting systems we should 
therefore prefix a negative sign to either of the expressions 
(1) or (2). 

3*81. Lost potential energy of a gravitating system. 

Let the system consist of particles of masses m l , w 2 , w 3 , . . . at 

points A lf A 2 ,A 3 , Let us find the work done by the mutual 

attractions of the particles in assuming these positions starting 
from a state of infinite diffusion. Let m l be in its final position 
A^\ the work done by the attraction of the particle m on the 
particle m 2 as it comes from an infinite distance to A 2 is 
flHra 2 /r 12 . Then these two particles attract the third and do 
work m 1 m 3 /r 13 4-m 2 w 3 /r 23 in bringing m 3 to A 3 and so on; so 
that the total work done may be represented by ^m p m q lr pqj 
where the summation extends to every pair of particles. 

T . TT ^9 m * 

Let j7 =: _2 + _J+.., > 



'21 '23 



Then V l is the potential at A l of w 2 , m z . . . ; V 2 is the potential 
at A 2 of m 1 , w 3 . . . ; and so on. 
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It follows that the total work done is equal, to SwF, 
the factor J being introduced because the sum SmF contains 
every term such as ^i^2/ r i2 twi ce - 

This expression JSraF represents the work done by the 
mutual attractions of the system of particles or the exhaustion 
of their potential energy as they come together from a state of 
infinite diffusion. 

If the system forms a continuous body, its lost potential 
energy can be found by dividing it into elements dm and 
taking the limit of JSFdra. Here F is, the potential at the 
position c$ dm 0f the rest of the body, and its limit is, by 
definition, the potential of the body at the point to which dm 
contracts. Thus the lost potential energy is 

Vdm, 

where F is the potential of the body at a point of itself and 
dm is an element of the body at the point. The potential 
energy of the body may be denoted by the same expression 
with a negative sign. 

It used to be supposed that the duration of the sun's heat might be 
calculated by assuming that, as the matter forming the sun condensed 
from a state of infinite diffusion, the mechanical potential energy lost 
was converted into heat, and then estimating how long it would take 
for this heat to be lost by radiation. 

39. Examples, (i) The density of a sphere varies as the depth below 
the surface ; shew that the resultant attraction is greatest at a depth equal 
to J of the radius, and that the value there is $ of the value at the surface. 

Let a be the radius, and let the density p at a distance r from the 
centre be given by p = A (a r). 

The attraction at a distance R from the centre, when R < a, is due 
to the sphere of radius R and its value is 



The greatest value of F occurs where R = %a, i.e. at a depth equal 
to Ja, and this greatest value is ^TrAa 2 ; whereas the value at the surface 
R = a is seen to be 



58 ATTRACTION AND POTENTIAL [3'9~ 

A 
(ii) Prove'that, if the law of potential be e~ r/A , the potential ofauniform 

thin spherical shell of radius a at an external point is the same as that of 
a particle of mass sinh T / r that of the shell placed at the centre. 

[M. T. 1918] 

Let P be an external point at 
a distance c from the centre O. 

Consider a narrow zone of the 
shell of breadth QQ' at a distance 
r from P. If 6 denotes the angle 
QOP, the area of the zone is 
2ira 2 sin 6 dO ; so that if m is the 
mass of unit area, the potential 
at P due to this zone is 




and the whole potential is 

"T- 

sin Od8. 



/"T 

J o 



o r 
But r 2 = a 2 + c 2 2ac cos 0, so that rdr ac sin 6d6, and 

F = 2irm - A [ C * a e-+fidr = 27rm - A\ (e-- fl >' A - e-< c+0)/A ) 
C J c-a c 

= 47rm - A Ae~ cyA sinh T . 
c A 

Also the potential at P of a mass 47rwa 2 at O is, with the same law of 

A 

potential, 4-Trwa 2 e~ c/A ; hence the result follows. 
c 

The reader will notice that this method of direct integration can be 
applied for finding the potential inside or outside the shell with the 
ordinary law of potential. 

(iii) Prove that the attraction of a uniform elliptic disc at a focus is 
2rrm{l ^/(l e 2 )}/e if the vanishing cavity round the focus is circular 
with the focus as centre, and zero if the cavity is bounded by a similar 
and similarly situated ellipse with the focus as centre of similarity. 

Let l/rl e cos 6 be the equation of the boundary ellipse and r = a 
that of the small circular cavity round the focus. 

The attraction is easily seen to bo given by 



ri/ 
Ja 



l-ecos0) C os0 ,., 
i- mrdBdr, 



and integrating with regard to r, this 

= m| cos 6 log T^J 7r\dO. 

jo fe a(l-ecos0) 
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Then, integrating by parts, this 

r / i I H 27r f^eain 

= m sin log 7= ^ -f m I ^ 

L 6 a(l-ecos0)Jo Jo l- 

_ /"27TJ i i_ e 2 i 

"~~ J 1 e ( 1 e cos 0) j 



But if the inner boundary , which gives the lower limit of integration 
with regard to r, be Z'/r = 1 e cos 0, the integral of with respect 

to r is cos 6 log Z/r, and integrating with respect to 6 the result is zero. 

(iv) Prove that the gravitational potential energy of a thin uniform 
circular disc of mass m and radius a is 
C 8ym 2 /37ta, where C is independent of a, 
and y is the constant of gravitation. 

[M. T. 1913] 

If a denotes the surface density, the 
potential at any point P on the edge of such 
a disc is easily seen to be 



fjTT 

= y a I 2a cos d9 = 4y era. 

J ~~ i 71 " 




Now suppose that the particles which make up the disc are initially 
in a state of diffusion at an infinite distance and that C denotes their 
potential energy in this state. Let the particles be assembled in con- 
centric rings. When the radius of the disc is r, the potential at its edge 
is, from above, 4ycrr, so that the work done in assembling a narrow 
ring of breadth dr is 4yor.27rordr. It follows that the work done by 
its own attractions in bringing together the whole disc is 



f 
/ 



f 

But w = 7r<ra 2 , so that the work done = 8ym 2 /37ra, and the potential 
energy remaining = C 8ym 2 /37ra. 

EXAMPLES 

1. Taking the earth as a uniform solid sphere of 4000 miles radius, 
shew that a pendulum clock would lose 10-8 sec. per day if taken to 
the bottom of a mine 1 mile in depth. [M. T. 1921] 

2. A small smooth cylindrical hole is drilled through a solid uniform 
sphere. Shew that a particle set free at the surface will oscillate in 
a period independent of the direction of the hole. 
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3 . The density of an infinite circular cylinder of radius a at a distance 
r from its axis is proportional to (a 2 r 2 )*. Shew that the attraction 
at an internal point is *, 

* 



where M is the mass of unit length of the cylinder, and y is the gravita- 
tional constant. [London Univ. 1926] 

4. A sphere is described on a radius of another sphere as diameter. 
If the latter be a homogeneous sphere of mass M and radius a, shew 
that the resultant attraction on the portion included in the smaller 
sphere is &M*/a*. [C. 1914] 

5. A solid sphere of radius a is such that its density at any point is 
proportional to the nth power of the distance of that point from the 
centre of the sphere. Find the potential at any external point, and 
shew that the potential at an internal point at a distance r from the 

centre is - y ft ( -- --T-O ) , where M is the total mass of the sphere. 
n + 2\ a a w+8 / ^ 

[London Univ. 1931] 

6. A uniform solid sphere of mass M is cut in two by a diametral 
plane; shew that the resultant attraction between the halves is 
$5yM 2 a~ 2 , where a is the radius of the sphere, and y the constant of 
gravitation. [C. 1906] 

7. Prove that the force required to separate the two parts of a solid 
uniform sphere divided in any manner is yMM ' . QG'IcP, where M, M' 
are the masses, and 6r, G f the centres of gravity of the two parts 
and a is the radius of the sphere. {C. 1905] 

8. P is a point on a uniform sphere of mass M and radius a. The 
sphere is divided into two parts by a plane perpendicular to the 
diameter through P. Show that, if the radius of the section by the 
plane subtends an angle a at P, the attraction at P of the portion 
of the sphere which contains P is 3yM cos 2 a ( 1 f cos a)/a 2 . 

[London Univ. 1937] 

9. A solid right circular cylinder is of density p and length I, and 
the radius of each end is a. Prove that the potential, due to it, at the 
ceAtre of each end is 

vy? (I V( 2 + Z 2 ) - 1 2 + 2 sinh- 1 (I/a)}, 
where y is the constant of gravitation. 

If in the above cylinder I = 8a/3, shew fchat the potential at the 
centre of the cylinder is 

27rypa 2 ($ + log e 3). [London Univ. 1935] 

10. Prove that the potential of a uniform circular cylinder, of 
length 2a, radius of base a, and mass M at the middle point of its axis, is 

yM { V 2 - 1 + log. ( V 2 + 1 )}/a. 

r London Univ. 1938] 
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11. P is a point outside a thin uniform spherical shell at a distance 
c from the centre. Prove that the sphere of centre P and radius c 
cuts the shell into two portions whose potentials at P are equal. 

[London Univ. 1933] 

12. A thin uniform spherical shell of density p has an internal radius 
a and external radius a -f a, where (a/a) 2 is negligible. Find the attrac- 
tion of the shell at a distance a + ka. from the centre, where 0<fc< 1. 
Hence shew that, if the outer radius contracts so that the shell becomes 
a thin layer of surface density or and radius a, the attraction of this 
layer at a point in itself is ^irak. How does this illustrate the discon- 
tinuity of normal force in crossing a surface? [C. 1928] 

13. Prove that if a solid sphere of uniform density p is divided into 
two parts by a plane which cuts the surface in a circle of radius 6, the 
Attraction between the two parts is |7r 2 yp 2 6 4 . [London Univ. 1936] 

14. A solid homogeneous sphere, of radius R and mass M, is divided 
into two portions by a plane at a perpendicular distance %R from the 
centre. Shew that the mutual attraction between the two portions 

27 M 2 
of the sphere is y ~ -~= . [London Univ. 1933] 

ZOD jrt 

15. From an infinite plate of thickness 2a of uniform homogeneous 
matter a spherical portion touching the two plane faces is cut away. 
Prove that the attraction of the plate at a point where the sphere 
touches either of its faces is f of its former value. 

[London Univ. 1927] 

16. A uniform solid sphere of mass M is placed near an infinite 
plate whoso surface density is uniform and equal to or. Prove that the 
sphere attracts the plate with a force 2-rryMcr. [London Univ. 1932] 

17. Shew that, if M be the mass of a sphere of radius a, the loss of 
gravitational potential energy in assembling the particles from a state 
of diffusion at an infinite distance is $yM 2 /a. [C. 1926] 

18. If the radius of a solid gravitating sphere be a, the density p, 
and an internal point P be at a distance b from the centre 0, prove that 
the difference of the potentials at P due to the two portions into which 
the sphere is divided by a plane through P perpendicular to OP is 

4 ^{a 8 ~(a 2 -& 2 )*}. [C. 1908] 

19. The density at distance r from the centre, in a sphere of attract- 
ing matter, is p Ar, where p Q is the density at the centre, an4 A is 
constant. Find the potential at any internal point, and prove that the 
potential at the centre is 



where a is the radiiis of the sphere, p l the density at the surface of the 
sphere, and y the constant of gravitation. [London Univ.] 
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20. If A B is a diameter of a spherical shell of uniform surface density, 
find the position of a plane perpendicular to AB which divides the 
shell into two portions such that their potentials at A are equal. 

If the surface density of the shell is a, find the attraction at its 
centre due to each portion of the shell. [London Univ. 1938] 

21. Shew that the initial rate of decrease of g in descending a mine 
shaft would be equal to g/a if the density (p) of the earth (radius a) 
were uniform. But if the oarth had a spherical nucleus of different 
density and of radius 6, the density of this nucleus must be 



where Xg/a is the initial rate of decrease of g in descending the shaft. 

[C. 1930] 

22. A sphere of radius a, mass M, and density varying directly as 
the distance from the centre, is built up of matter brought in small 
amounts from a state of complete dispersion at infinity. Shew that the 
work W made available during the process is given by "FP = fyM 2 /a, 
where y is the constant of gravitation. 

Shew that, if the matter is now redistributed so as to form a sphere 
of the same radius but of uniform density, a further amount of work 
is made available. [P. 1936] 



23. Find the total work done by the gravitational forces in collecting 
from infinite dispersion six equal particles of mass m and placing them 
at the corners of a regular hexagon of side 2a. [London Univ. 1932] 

24. Six particles, each of mass m, are placed symmetrically on the 
circumference of a circle of radius r and are allowed to move from rest 
under their mutual attractions. Find the velocity of each particle 
when they are on a circle of radius ^r. [London Univ. 1933] 



25. Assuming the earth to bo a homogeneous sphere of density a 
a^d radius a, surrounded by a spherical shell of density p and thickness 
6 a, shew that the weight of a body given by a spring balance is 
unaltered at first on descending vertically from the surface if 
p = 2a 3 cr/(2a 3 + 6 8 ). The rotation of the earth may be neglected. 

[M. T. 1914] 

26. A solid sphere of radius a and uniform density is surrounded 
by a concentric shell of internal radius a and external radius b in which 
the density is a function of the distance from the centre. Find the law 
of density if the force of attraction is constant in the substance of the 
shell, and prove that the mass of the shell is to the mass of the sphere 
as b 2 - a? is to a 2 . [C. 1932] 
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27. The density of a sphere of radius a and mean density p is a 
function of the distance from the centre such that the attraction at 
any point inside the sphere is proportional to the square of the distance 
from the centre. Express the density at a distance r from the centre 
in terms of p, a and r, and prove that the potential at an interna 

point at a distance r from the centre is s ^ (4a 3 - r 8 ). [P. 1933] 

u Qi 

28. If be the centre of a solid uniform sphere and P an external 
point, shew that the sphere on OP as diameter cuts off a portion whose 

1 M 

attraction at P is equal to ^ n P2 > where M is the whole mass. 

2 U [C. 1905] 

29. Find the potential at any point in the interior of a sphere of 
radius a, the density of the sphere at a distance r from the centre being 



Shew that, if the sphere be divided into two parts by a plane through 
the centre, the mutual pressure between the parts is 



[London Univ. 1934] 

30. A portion of a uniform solid paraboloid of revolution is cut 
off by a plane at right angles to the axis OX. Prove that the attraction 
at the focus S will be zero when sin(^P$J")= l/V e P being a point 
on the rim. [M. T. 1905] 

31. If p be the density of a homogeneous solid sphere and a the 
radius, prove that the force in astronomical units exerted by the rest 
of the sphere on a conical portion of itself bounded by a right circular 
cone of angle 2a with its vertex at the centre is ^7r 2 p 2 a* sin 2 a. 

[C. 1915] 

32. A solid of uniform density p is in the form of the spheroid 
obtained by rotating an ellipse, of latus rectum 21 and eccentricity e, 
about the major axis. Prove that the potential and the intensity of 
the attraction at a focus of the generating ellipse are 

2< rr ypl*/( 1 - e 2 ) and (2wy/>Z/e a ) [log {( 1 + e)/( 1 - e)} - 2e] / 

respectively. [London Univ. 1936] 

33. Find the potential and attraction, at any point, of a homo- 
geneous sphere with a spherical excentric cavity, and deduce the 
corresponding values for a thin spherical shell whose small thickness 
at any point is proportional to its distance from a plane. [C. 1899] 

34. Two uniform solid spheres are of density p and radius a, and 
the distance between their centres is c ( > 2a). Prove that they attract 
each other with a force 167r 2 yp 2 a 6 /9c 2 , 

where y is the constant of gravitation. [London Univ. 1926] 
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35. A solid sphere of radius a and mass M is of uniform density. 
Prove that the attraction on an octant of the sphere due to the rest 
of the sphere has components parallel to the edges of the octant each 
equal to ^yM 2 /a 2 , where y is the gravitational constant. [C. 1914] 

36. Prove that the work done by the gravitational forces when 
a uniform spherical shell, of mass M and of external and internal 
radii a and 6 respectively, is changed to a solid sphere of the same 
density is 

5 - 2a 5 - 36* + 5a 2 6 3 ], 



where c 3 = a 3 - 6 3 . [London Univ. 1939] 

37. The intensity of gravitational attraction at any point on the 
surface of a certain uniform solid sphere, of radius a and centre O, is g. 
A spherical cavity of radius b is made with its centre at a point C, 
where 0(7 = c, (b + c<a). Prove that the magnitudes of the intensities 
at the points where the line OC (produced) meets the surface of the 
sphere are 



and that at the ends of any diameter through O perpendicular to OC 
they are 



[London Univ. 1935] 

38. A uniform spherical shell of gravitating matter, in which the 
volume of the matter is equal to that of the hollow it encloses, contracts 
under its own gravitation into a uniform solid sphere of the same 
density. Prove that the work done by the gravitating forces is 
3(2$ I)yjlf 2 /10a, where M is the mass and a the radius of the solid 
sphere. [London Univ.] 

39. The density of a sphere, of mass M and radius a, at a distance r 
from the centre is proportional to a r. Prove that the potential at 
this distance is yM (2a 8 2ar 2 + /* 3 )/ 4 , and that the potential energy 
exhausted in collecting the elements of the sphere is 26yM 2 /35a, where 
y is the gravitational constant. [London Univ. 1926] 

40. Find the negative potential energy W for a spherical distribution 
of mass M and radius a in which tho density varies uniformly from p 
at the surface to p at the centre, and deduce that, if p^ = 0, 



while if /> = 0, pi^O this value of W is reduced in the ratio 10/13. 
O denotes the gravitation constant. [M. T. 1928] 

41. Prove that, if a segment of a uniform thin spherical shell has 
an angular radius a, its attractions on the two sides, at points close to 
the pole of the segment, are 27rya(l sina), where a is the surface 
density. [M. T. 1922] 
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42. Defining the principal value of the attraction at a point of an 
attracting surface as the limit of the attraction at the centre of a 
circular hole of which the radius tends to zero, shew that its tangential 
components are continuous with those of the attraction near the 
surface, whMe its normal component is the mean of the limiting normal 
components on either side. 

Shew that the attraction at the focus of a small elliptical hole of 
eccentricity e in a thin uniform spherical shell makes an angle 



tan 



" 1 



with the radius of the shell. [C. 1929] 

43. Prove that the work done by a gravitating system in attracting 
a new body from infinity is equal to I p ( V + V) dr, where the integral 

is taken through the space occupied by the new body, and F, V are 
the gravitational potentials before and after the new body arrives. 

Calculate the work done in filling up a spherical hollow of radius 6 
in a uniform sphere of radius a, the distance between the centres being 
c and the density p. [C. 1914] 

44. Find an expression for the mutual gravitational potential energy 
01 a uniform sphere of mass M and a uniform rod of length 2a and mass 
m when they are placed in any relative position. 

Shew that this is constant when the centre of tho sphere moves on 
the spheroid formed by rotating an ellipse whose foci are the ends of 
the rod about the rod. [M. T. 1923] 

45. A uniform rod BC of mass M is free to turn about its middle 
point O which is fixed, and a uniform sphere of mass m has its centre 
at A. Shew that the couple on^the rod tending to turn it into the 

line OA is \yMrn ( j - tan Ja ) , where b, c denote the lengths of 

\ oc J 

AC, AB and a the angle BAG. [M. T. 1930] 

46. The profile of a long straight mountain range of uniform 
density p is an arc of a circle less than a semicircle. Prove that if h 
is tho height of the range and 6 its breadth at its base, a plumb line at 
its foot will be deflected towards the range through an angle approxi- 
mately equal to 3 b 2 fi 

-'* 



where R and p are the radius and mean density of the earth. 

[London Univ. 1927] 

47. Shew that, if the potential of a particle of mass m at distance r 
is m<f>" (r)/r, then the potential of a uniform solid sphere of mass M 
and radius a at an external point at distance r from the centre is 

- )] - W (r + a) - * (r - a)]}, 
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where the functions <f>' (x) t </>"(%) , </>'"(%) are successive derivatives 
of </>(x). 

Consider in particular $'" (x) = ye~~ A *, where A^O. Shew that the 
external field of the sphere is the same as that of a single particle at 
its centre; and that the mass of this particle is 

3M (Aa cosh Aa sinh Aa)/A 8 a 8 
if A>0,orMif A = 0. [M. T. 1924] 

48. Shew that the potential of a uniform circular disc at a point of 
,self is 



r 

Jo 



fo vv 

where r is the distance of the point from the centre and a is the mass 
per unit area. 

Deduce that the work done in removing all the particles to infinity is 
|tra a a 8 . [C. 1924] 

49. If the attraction between two elements of matter of mass m 
and m' is mm' f(R) 9 where R is the distance between them, shew that 
the inward acceleration due to a uniform thin spherical shell of mass M 
and internal radius a, at an internal point distant x(<a) from the 
centre, is F, where 



If F is zero for all values of a and x, subject to the condition x<a, 
shew that the only possible form of f(R) is AjR z , where A is some 
constant. [M. T. 193" 
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this expression is equal to mda); -f or according as 6 is acute or 
obtuse, i.e. -f or according as the cone proceeding from. O is entering 
or leaving the region bounded by S. It follows that when O is outside 
the region (as in the figure), since the cone leaves the region as often 
as it enters it, therefore the total contribution to the outward normal 
flux arising from the intersections of this cone and the surface is zero. 
But when lies inside S the cone leaves the region and if it re-enters, 
it finally emerges, so that it leaves the region once more often 
than it enters it, and the total contribution in this case is mdw. 
Hence by taking cones in all directions round 0, we see that a particle 
of mass m contributes zero to the outward normal flux when it lies 
outside the surface, and 47rw when it lies within the surface. And 
by summing for all particles of the system we obtain the required 
result, viz. total outward normal flux = 4?r (mass enclosed). 

4*22. When the field is due to the attraction of a continuous body 
and the surface S intersects the body, 
we may regard the body as divided into 
distinct parts of masses M t and M , inside 
and outside S, and suppose that these are 
separated by an infinitesimal distance so 
that MI lies wholly within S and M 
wholly without it. The outward normal 
flux of force from S is then 4?rM f as 
before, and this will continue to bo true 
as the distance between M l and M 
decreases without limit. 

The result is not affected by the ' matter 
on the surface *S" because wo are con- 
sidering a body of finite volume density 

and, as explained in 37, it requires a finite surface density of matter 
(i.e. an infinite volume density) to affect the normal component of 
force across the surface. 

There is, however, an exceptional case in which the theorem needs 
modification, namely when the matter is a distribution of finite surface 
density on the surface S. Reverting to 421 we observe that if the 
point O lies on the surface S, then only cones on one side of the tangent 
plane at S intersect the surface, so that the contribution of the mass m 
to the total outward flux is 27rw, and the same being true for all 
matter lying on the surface, the total outward flux due to a total mass 
M is - 27rM. 

It is easy to verify, for a uniform distribution of mass M on the 
surface of a sphere S, that the total outward flux of force 




and 



across a concentric sphere inside S is zero (3'2), 

across S is - 2<rrM (371), 
across a concentric sphere outside S is 4wM (3*2). 
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4*23. Applications of Gauss's Theorem to spheres and 
cylinders. 

(a) Spherical shell. Let M be the mass of a uniform spherical shell 
of radius a. Apply Gauss's theorem to a concentric sphere S of radius r. 
By symmetry the outward normal force R has the same value at all 
points of S. 

When r<a, fig. (i), S encloses no matter, so that 

7rr 2 R 0, and R = at all points inside the shell. 
When r > a, fig. (ii), iS encloses a mass M , so that 
47rr 2 J? = - 477.M, and R = - M/r 2 at all points outside the shell (3-2). 





(6) Solid sphere of uniform density p. To find the attraction at an 
internal point, apply Gauss's theorem to a concentric sphere S of 
radius r < a, fig. (i). S now encloses a mass ^rrpr 3 , so that 

477^^ = - 47r . $777>r 3 and R = - %npr (3-3). 

The attraction at an external point may be found as for a uniform 
spherical shell. 

(c) Long uniform cylinder. Let the figures now represent cross- 
sections of a uniform long cylinder of radius a and mass M per unit 
length. At points sufficiently distant from the ends of the cylinder 
we may neglect force parallel to the axis, and as above let R denote 
radial force. Then we may take for S a unit length of a coaxial cylinder 
of radius r and ignore any force across its plane ends. 

When r<a, fig. (i), S encloses no matter, so that 



0, and R = at all points inside the cylinder. 
When r>a, fig. (ii), S encloses a mass M, so that 
2irrR = 4?rM, and R = 2M/r at all points outside the cylinder (2*42). 

4*24. Tubes of force. If the lines offeree (2-26) be drawn 
through every point of a small closed curve drawn on an 
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equipotential surface, these lines form a tube of which the 
equipotential surface is a cross-section, called a tube of 
force. 

Apply Gauss's theorem to the region bounded by a portion 
of a tube of force cut off be- 
tween two equipotential sur- 
faces and not enclosing any 
matter. If R l9 R 2 denote the 
force of the field at the ends of 
the tube and a> ly o> 2 the areas 
of the ends, since there is no 
flux of force through the sides 
of the tube and no matter 
inside it, therefore 




I 



R 



or the force varies inversely as the area of the cross-section 
of the tube. 

Example. If the field of force be due to the attraction of a uniform 
infinite plane sheet of matter of mass m 
per unit area, it is evident that the equi- 
potential surfaces are planes parallel to 
the given plane and the tubes of force are 
straight, so that the force is constant. If 
jR denotes the force directed away from x rp 
the plane and we apply Gauss's theorem R ^ 
to a portion of a tube of force of cross - 
section o> extending to both sides of the 
plane, since by symmetry the force has 
the same value on both sides of tho 
plane 

4*3. Theorems, on the potential, (i) The potential cannot 
have a maximum or a minimum value at any point of space 
unoccupied by matter. 

For if the potential had a maximum value at a point 
we could surround by a small sphere S such that at every 
point of S the potential would be less than at 0. The outward 
normal force R n would therefore be negative at all points of 
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8, and R n dS taken over 8 could not vanish, as it must if 8 

contains no attracting matter. Similarly for a minimum. 

As a consequence of this theorem we note that at a point 
in free space, if the potential be not constant in the neighbour- 
hood of the point it increases in some directions and decreases 
in others, so that, although a particle placed at the point 
free from constraint might be in equilibrium its equilibrium 
could neither be stable, nor unstable for all displacements. 

This latter theorem is due to Earnshaw* 

(ii) // the potential be constant over a closed surface contain- 
ing none of the attracting matter, it must be constant throughout 
the interior. For otherwise it would have a maximum or 
minimum value in the interior, which is impossible by (i). 

* (iii) Gauss's Mean Value Theorem. The mean value over 
a spherical surface of the potential of any attracting system is 
equal to the potential at the centre of the sphere due to the mass 
outside the sphere, plus the mass inside the sphere divided by 
the radius. 

Let a be the radius of the sphere, m the mass of a particle 
at an external point A and m 1 
the mass of a particle at an 
internal point A'. - 

The potential at any point Q 
on the sphere is 

m m' 




The mean value of V over the sphere is - ^ VdS, where 



dS denotes an element of area at 
This 



r d Sf 

But ~j-^ is the potential at A due to a spherical layer of 
J AQ 

unit surface density and therefore = 4c7Ta*/OA. 



P 

Similarly, I - 
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/J Q 

-JTT\ ^ t^ 6 potential at A' due to the same layer 

and therefore = 4?ra 2 /a = 4?ra. 

TT . , , r rr m m' 

Hence the mean value of V ~ -^-7 + 

OA a 

In like manner the mean value of the potential due to any 

attracting sy stem = 2 -^-j + - - , where the first sum repre- 
\JA a 

sents the potential at O due to the matter outside the sphere, 
and thus the theorem is proved. 

(iv) // the potential be constant throughput any region T, 
unoccupied by matter, it has the same constant value through 
all space that can be reached from T without passing through 
matter (Gauss). 

Let V be the value of the potential throughout T. If the 
potential differs from F in neighbouring space there must be 
some regions contiguous to T where the potential is greater 
than F and others where it is less. In the neighbourhood 
of a place where the potential is greater than F draw a sphere 
with its centre inside I 7 . Then on the part of the sphere 
inside T the potential is F and on the part outside it is greater 
than F, so that the mean value of the potential over the 
sphere exceeds F, the value at the centre of the sphere, 
which is impossible by the last theorem since the sphere con- 
tains no matter. Hence the region T of constant potential 
must be extended to include the whole of the sphere, and 
so on step by step through all space that does not include 
attracting matter. 

* (v) In any case of symmetry about an axis, if the potential 
be constant through any finite distance along the axis, however 
short, not within the attracting matter, it has the same constant 
value through all space that can be reached from this portion 
of the axis without passing through matter* Without loss of 
generality we may suppose that the potential is zero at all 
points of a segment of the axis. Let P be a point of the 

* Thomson and Tait, Treatise on Natural Philosophy, Art 498, Oxford 
1867. The proof given above is due to Dr S. Verblunsky. 
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segment. On account of the axial symmetry, there is no force 
in any direction through P perpendicular to the axis. Hence 
the axis is normal at P to a surface on which V is constant ; 
i.e. F = Pp = 0. Since the argument applies to any point P of 
the segment, it follows that F = throughout a region. The 
theorem therefore follows from (iv). 

4*31. Comparison theorems, (i) // two different distribu- 
tions of matter have equal potentials over any dosed surface not 
including any attracting matter, they have equal potentials 
throughout the space enclosed and through all external space 
y;hich can be reached without traversing any of the attracting 
matter. 

Let the attraction of one of the distributions be changed 
to repulsion. This changes the sign of its potential and makes 
the potential due to the two distributions taken together zero 
over the given surface. Therefore by 43 (ii) the potential is 
zero throughout the space enclosed, and by 4-3 (iv) zero 
throughout all external space which can be reached without 
passing through matter. Then changing back from repulsion 
to attraction, it follows that the potentials of the two distribu- 
tions are equal. 

(ii) // two different distributions of matter produce equal 
potentials over any surface enclosing both, they produce equal 
potentials throughout all space external to the surface. 

Let the attraction of one of the distributions be changed 
to repulsion. The two distributions combined then produce 
zero potential over the given surface, and they also produce 
zero potential over the infinite sphere. Therefore their com- 
bined potential is zero throughout the space between the given 
surface and the infinite sphere, for this space contains none 
of the matter (4*3 (ii)). Reverting to the original distributions, 
it follows that they have equal potentials throughout the 
same region. 

(iii) // two finite distributions of matter have the same equi- 
potential surfaces outside both distributions, their attractions at 
every external point have the same direction and are proportional 
to the masses of the distributions. 
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The attractions at any external point are in the same 
direction because they are along the inward normal to the 
equipotential surface through the point. 

Let F, V denote the potentials of the two distributions. 
Since V is constant whenever V is constant, there must be 
a functional relation between them; say V =f(V). Therefore 

dV '' dV 



Adding the similar results obtained by differentiating with 
regard to y and 2, we get 



But, outside the matter, 



Hence, unless F is constant everywhere to make the second 
factor vanish, we have f <v\ c\ 
giving, on integration, 



where A, B are constants. But F and V vanish at infinity 
in the ratio of the masses M , M ' of the distributions, so that 

M' 

B = and A = M'fM. Hence V = -=-= F, and the result follows. 

4*32. Condition that a family of surfaces is a possible 
family of equipotential surfaces in free space. A family 
of surfaces in general is riot a family of equipotential surfaces 
unless its equation satisfies a certain condition. To find the 
condition that the equation 

/(#,2/,z) = const. 
may represent a family of equipotential surfaces. 
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If the potential V is constant whenever /(x, y, z) is constant, 
there must be a functional relation between V and/ (x, y, z) ; say 



- 



By adding similar results in y and z, we get 



But in free space V 2 F = 0; therefore 

f(/j , ,. 

i = - j>-77r = a function of /, say x (/) 



ay ey 



/a/\ 2 /a/\ 2 /a/\ 

(dx) + (fy) + \dz) 



(1). 



This is the necessary condition, and when it is satisfied the 
potential F can be expressed in terms of f(x, y, z) thus: 
We have F =<(/), where 



Therefore 

or f (/H 

It follows that V=<f>(f) = A IV'*^' d/+ 5 ............ (2). 

4*33. Example. Shew that a family of right circular cones with a 
common axis and vertex is a possible family of equipotential surfaces, 
and find the potential function. 

Taking the axis of z for the common axis, the equation of the family 
of cones is 2 . ^2 

f(x,y,z)= -y =const ................... (1). 
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The condition 4-32 (1) in this case becomes 



22 
* a 



__ 2-f3/ . 
**~ 2l' 



and the condition is satisfied since the expression on the right has 
reduced to a function of/. Hence 

t'U) i i _ 



and, by integration, 

) = const. 



Put / = tan 2 6, and we get d^ = : j . 

Hence V = <^ (/) = 27 log tan 

It is clear that 6 is the half-angle of a cone, and V is constant when 6 
is constant. 

The truth of this proposition is evident at once if we consider 
Laplace's equation in polar co-ordinates, namely 



r 2 sin 2 00< 2 ' 

and whether it can have a solution which makes F constant over the 
surface of a right circular cone. There is clearly such a solution if V 
is independent of r and </> and satisfies 

JV A 



where A is a constant, and this leads to 

F = .41ogtanJ0+B, 
as before. 

4-4. When the potential is given, to find the distribution 
of matter. We have seen that potential functions may become 
infinite at points or lines where matter is concentrated, and 
that otherwise potential functions are finite and continuous, 
even when crossing surfaces on which matter is concentrated 
(3*7); but in crossing such a surface there is a discontinuity 
in the gradient of the potential. Thus a closed surface on 
which matter is concentrated divides space into regions in 
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which the potential function takes different forms, though 
these forms take equal values on the surface itself. 
When the potential is given at all points of space, Poisson's 

equation, in the form />= V 2 7, serves to determine the 

volume density of matter wherever a finite volume density 
exists. If the potential is given by different functions V ly V 2 
on opposite sides of a surface S, this implies that there is 
on 8 matter of surface density a determined by the relation 



where the direction of the normal dn is from 1 to 2. 

If the potential becomes infinite at a point, the amount of 
matter concentrated there can be found by applying Gauss's 
theorem to a small sphere having its centre at the point. 

If the potential becomes infinite at all points of a line, 
the line density of the matter can be found by applying 
Gauss's theorem to a short cylinder whose axis is the line and 
whose radius tends to zero. 

4*41. Examples, (i) The potential outside a certain cylindrical 
boundary is zero ; inside it is V = a; 3 3xy* ax 2 -f 3ay 2 . Find the 
distribution of matter, 

We have first to find the boundary. Since the potential is continuous 
across the boundary and zero outside, the boundary must be given 

by 

A 



or (x a}(x 

i.e. the section is an equilateral 
triangle OA B of height a. Then 

2 F, ~ 2 F, 

g 1 = Qx 2a, "Q~T = ~~ ^ x "^~ > 

and ~^Y ^ ^ 

so that inside the prism 



-7T 7T 

and outside p = since F 2 = 0. 
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A* r> AT> l 

At PonAB cr= -~ T - 
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= -; (a 2 - 3^/2) = . 



a= AP.PB. 

47T 



At Qon^lO (x= 

= L 

a 4ir a 



2 dx 2 

c-3 

f' 



J ^ 
+ 3 



and similarly on OB. This shews that a solid prism of uniform density 
a I tr would produce the same external field as a distribution of matter 
of surface density SAP . PB/^ir on each of the faces of the prism, 

(ii) Use the theorems of Laplace and Gauss to verify that the function 
V m log 7 ^ is the potential due to a uniform line density m on 

the straight line joining the points A, B from which r, r' are measured. 

Take the origin at the middle point of AB and the axis Ox along BA, 
and let f denote distance from Ox. 




BOA x 

We notice that V is a finite and continuous function except on the 
line between B and A where it becomes infinite. 



We have 
and F = 



r 2 = (x - a) 2 



r /2 = (x + a) 2 + J 2 , 



6-2 
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Hence ^-=^~- ( ^ + l) ( + 

dx r' + x-fa\ r' / r+x-a\ r ) 



m m 



&V m(aH-a) { m(x-a) 



'r' r + x-a'f 



x a r x + a\ 
r 



Similarlv ' 
similarly, 



or -g = 

1 d f Y dV\ m(x-f-ci) m(x a) 



_m/r' x a 
~T\ r' 



By adding (1) and (3) we get 

a 2 F 1 d 



shewing that, since V is symmetrical about Ox, it satisfies Laplace's 
equation in cylindrical co-ordinates <, #. This holds good at every 
point save at points between A and B on the line AB where V is' 
infinite, so there is no volume density of matter; and since there is no 
discontinuity in the form of V there is no surface density. Hence the 
matter is confined to the line AB. It remains to shew that the line 
density is uniform. 



6 

K 

Apply Gauss's theorem to a small cylinder of length e and radius f 
with its centre at any point P on AB and its axis along AB. The flux 
across the curved surface is 

^V 

~ - = 2 



/ x + a , x a\ /n . 
( -- 7 4- U from (2), 



~> 47rem as 
and the flux across the plane ends is 7r 2 *- 2 - , which tends to zero 

with . It follows from Gauss's theorem that the mass on any element 
of AB of length c is me, so that the line density is uniform and equal 
to m. 

4*5. Simple applications of Laplace's equation. 

(i) Field due to a uniform plane sheet. Taking the axis of z at right 
angles to the sheet, V is clearly a function of z only, so that Laplace's 
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equation reduces to d 2 F/dz 2 = 0, which gives on integration dV/dz = A, 
where A is the constant value of the force in the direction in which z 
increases; and, by 25 or 424, this is 2frw, where m is the mass per 
unit area of the sheet, so that 

dv 
-j- = 

dz 

and F = 2-rrW2 + const. (251). 

(ii) Symmetry about a point. When the field is symmetrical in 
such a way that F is a function of r only, in three dimensions, we can use 
the transformation of Laplace's equation into 
polar co-ordinates given inl68 , putting d V/86 = 
and dV/d(j) = 9 or we can obtain the required 
form directly by applying Gauss's theorem to 
the space between concentric spheres of radii r 
and r + dr. Assuming that this space contains 
no matter and that the force at distance r is 
dV/dr, the flux of force out of the region across 

the surface of the inner sphere is 47rr 2 7 - . 

^ dr 

The outward flux across the sphere of radius r -f dr is therefore 

4 * dv , d (4 * d V\j 
47TT 2 r~ -f ,~ 477T 2 -j- dr 

dr dr\ dr J 

to the first power of dr, and the vanishing of the total outward flux gives 




dr 

XT u * * * dv A dV A 

Hence, by integration, r 2 -, =^L, or = -g, 



and v=- + B ........................... (1), 

where A and B are arbitrary constants to be determined to suit each 
particular problem. For example, consider the field due to a uniform 
thin spherical shell of radius a and mass m per unit area. Denote the 
potentials inside and outside the shell by V l and F 2 . The conditions 
to be satisfied by F x and F 2 are given in 4- 13, namely 

where r < a, V 2 Fi = 0, so that from ( 1), V l = - -- 1 + B^\ 



wherer>a, 

T7 T7 A dV * 

~-~ -, 



, T7 T7 A dV * dV l 

where r = a, V l = F a and ~-~ -, = 
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The last two conditions give 



(3). 



We also know that, for large values of r, V tends to zero like 4:7rma 2 /r, 
so that B 2 = and ^4 2 = 47rma 2 . Hence from (3) ^ = 0, and from (2) 
B t = 4:irma; so that F x = 4?rma and F 2 = 



' (iii) Symmetry about an axis. When the field is such that V is 
a function of r only, where r denotes distance from an axis; we may 
either quote the transformation of Laplace's equation into cylindrical 
co-ordinates given in 1-68, putting dV/dd=Q and dV/dz = 0, or we may 
apply Gauss's theorem to unit length of the space bounded by cylinders 
of radii r and r + dr. In this case, the flux out of this region across the 

dV 
cylinder of radius r is 2?rr -,- , so that the corresponding flux across 

the cylinder of radius r -f dr is 2rrr _- -f -.- ( ZTTT , ) dr to the first power 
of dr, and since there is no matter in this region 
d ( dV\ A dV A 

-j-[ r 'j-) Q> or r -j~A, 

dr\drj dr 

so that V = Alogr + B, 

where A, B are arbitrary constants to be determined to suit each 
particular problem. 

If the matter be a long uniform thin cylinder of radius a and mass m 
per unit area, we may take for the internal and external potentials 



When r = a, we have V l = F 2 , so that 

.................. (1); 



, , 2 x ,, , 

and when r = a, ~~ --, - = 47rm, so that 
dr dr 



a a 



Also applying Gauss's theorem, as in 423 (c), to a cylinder of radius 
r(>a), we have 



~ 
dr 



so that ^4 2 = 4irma. 



Then from (2) A l = 0, and from ( 1 ) ^ = 47rma log a -f 5 2 , shewing 
that the forms for V^ and F 2 are in agreement with the results of 2*42. 
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>. Potential of a body at a distant point. MacGullagh's 
Formula.* Let G be the centre of gravity of the body and P 
a point whose distance R from G is large compared with the 
dimensions of the body. Let ra be the mass of an element 
of the body, situated at Q, where GQ = r and the angle 
QGP = 0. 




The potential at P is given by 



'PQ' 



y m/ 
~~ E\ 



r 3r 2 cos 2 



___. _ __ ..., 

But if M denotes the mass of the body, A, B, G its principal 
moments of inertia at G and / its moment of inertia about 
GP, we have 

2m = M 9 Srar cos = (since G is the centre of gravity), 

2Z t mr 2 = A + B + C and Smr 2 sin 2 = /. 

_ r T7 M A + B+C-3I 

Therefore F== "^" f -- ^ - ............... ( X ) 

is an approximation to the value of the potential. 

4*61. Attraction components at a distant point. Taking 
axes through the centre of gravity, let x, y, z be the co-ordinates 
of the distant point and r its distance. In 4*6 this means that 

* James MacCullagh (1809-1847). Irish mathematician and 
physicist. 
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GP is now of length r and its direction cosines are xfr, y/r, z/r, 

so that I=(A 

and 



+ B-2C)z*} ...... (1), 

to this order of approximation. 

The attraction components at P are therefore given by 

9F__.Mx (B + C-2A)x 
~~~~~ + 



(2), 

and similar expressions for F and Z. 

Since action and reaction are equal and opposite it follows 
that a unit particle at P exerts on the body M forces equal and 
opposite to X, F, Z acting at P. To measure their effect on 
the body it is convenient to transfer them parallel to them- 
selves to act at O by introducing suitable couples. Thus the 
reaction on the body is represented by a force at G with 
components X , F, Z, and a couple with components 

-v ..**(C-B)y*\ 



3(A-C)zx 



. 3(B-A)xy 



.(3). 



462. As an example of 461 (3), if the body be the earth regarded 
as an oblate spheroid, and a distant body be regarded as a particle of 
mass M ' at distance r and angular distance 8 above the equator, we 
may without loss of generality suppose the distant body to be in the 
plane zx, then the couple exerted by the distant body on the earth 
tending to cause rotation about a diameter (the y axis) is 

3M'(A-C)xz BM'(A-C) . 

- v__ l f or 2- - 



and this has its greatest value when the declination 8 is JTT. 



4-7] MUTUAL POTENTIAL ENERGY 89 

4*63. Mutual potential energy of two bodies at a great 
distance apart. Let M , M 1 

be the masses and G, G' the 
centres of gravity of the 
bodies, and let GO' =r. Then 
if m'is the mass of an element 
of the second body at P and 
V the potential of the first 
body at P, the mutual potential energy is, by 3-8, equal to 




v , T7 _ 
- Sm'F = - Sm' 



M 



+ 



f Aa 
from 4-6, 



where I P denotes the moment of inertia of the body M about 
GP\ and this 

m' A + B + C- 37 



correct to the order -=, where / is the moment of inertia of the 
r 3 

body M about GG'. 

m' 
But S -~p is the potential at G of the body M ', and by 4*6 



. ., D/ ^ ., , 
, where A', B r , C l , /' bear 



,, 
is equal to + 

the same relation to the body M' that A, B, C, I do to the 
body M . 

Hence we get for the mutual potential energy 



f \M r A' + B' + C'-ZI' 

I \T + ~ ^ ~ 

VMM' 
~~~ 



J 



J' 



4*7. Gentrobaric bodies. When the action of terrestrial 
or other gravity on a rigid body is reducible to a single force 
always acting through a point fixed relatively to the body, 
that point is called the centre of gravity of the body and the 
body is said to be centrobaric. 



90 
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It follows that uniform spheres and spherical shells are 
centrobaric. In elementary statics it is usual to assume that 
the attractions of the earth on the elements of a body are 
a set of parallel forces with a resultant acting through a point 
fixed relatively to the body, so that on this hypothesis every 
body has a centre of gravity and its position is defined by 

Mx = %mx, My ^Lmy, Mz = ^mz ...... (1). 

We shall now prove the following theorem: // the resultant 
attraction of a body at all external points passes through the 
same fixed point, that point must be the centre of gravity, as 
defined by (I), and the body must be such that every axis through 
the centre of gravity is a principal axis and all principal moments 
of inertia are equal. 

Take the fixed point as origin 
O. Let n be an element of mass 
<rf the body at Q whose co- 
ordinates are #, y, z. Let P be 
an external point at distance r 
from O and let I, m, n be the 
direction cosines of OP. Then 
for the^ potential at P we have 




'PQ- 



I 



S/i, (Ix + my + nz) 1 
_ _ 



__ 2)2}* 



a; 2 + y* 4- z 2 
__ 



3 S/x (Ix + my -f- nz) 2 

I 7T Q "i 

2 r 3 



4'8] POINTS OF EQUILIBRIUM 91 

where A, B, (7, D, E 9 F are the principal moments and the 
products of inertia. 

The resultant force at P is to pass through for all positions 
of P, so that, if Rewrite Z = cos0, m = sin#cos</, 7fc = sin#sin</>, 
we must have 97/90 and 9F/9</> both zero for all values of 0, </>. 
This can only be so if the coefficients of the different powers 
and products of I, ra, n vanish separately, so that 

= = = 0, A = B=C and D = E = F = 0, 
which proves the proposition. 

4 8. Points of equilibrium. If in the presence of attract- 
ing bodies the resultant force on a free particle is zero, the 
position of the particle is a point of equilibrium. If V denotes 
the potential at the point P of co-ordinates x, y, z, then P 
is a point of equilibrium if 

8 !-?Z-?Z-o m 

~~ ..................... ( h 



for these are the conditions that the resultant force at P is 
zero. 

Regarding F as a function of x, y, z, the equation 

F (a, y, s) = const ................... (2) 

represents a family of equipotential or level surfaces. Using 
suffixes to denote partial differentiation, the direction cosines 
of the normal at (x, y, z) on such a surface are proportional to 
V x ,V y9 V a , which shews that the resultant force at the point 
is normal to the surface (2-26). It follows that if the point is 
a point of equilibrium, it is a singular point on the equipotential 
surface; i.e. there is not a unique tangent plane at the point, 
but every straight line through the point cuts the surface (2) 
in two coincident points. 

Let V be the potential at a neighbouring point P' of co- 
ordinates x 4- , y + 77, z -f , then 

V = 



(B), 
where the derivatives denote values at (x, y, z). 
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The result (3) incidentally affords another proof of the 
theorem 4-3 (i) that the potential cannot have a maximum 
or a minimum value at any point unoccupied by matter. For 
a first condition that V may be greater (say) than V for all 
small values of , 77, (positive or negative) is the vanishing 
f V x , V y , V g , and the further conditions require that V xx , V yy , 
V zz shall all be positive; which is not possible since V 2 F = 0. 
At a point of equilibrium there will therefore be some directions 
in which V > V and others in which V < V, so that equilibrium 
would be stable .for some displacements and unstable for 
others [4-3 (i)], and these will be separated by the surface for 
which F' = F; and from (3), since V x = V v = V e = Q, this is, 
for small values of f , 77, , a quadric cone 

F,,P + F^ + F^ 2 + 2F^^ ...... (4). 

Further, since V xx + V y y + Vgg = Q, it is a cone with three 
perpendicular generators. 

If two sheets of the same level surface intersect one another 
along a line, every point of this line is a point of equilibrium, 
because there are two normals to the surface at such a point 
each of which would be the direction of the resultant force 
if any. In such a case the tangent cone to the surface, repre- 
sented by (4), breaks into two planes; so that (4) is equi- 
valent to 



where IV + mm' + nri = 

so that, if two sheets of a level surface intersect, they intersect 
at right angles. 

49. Examples, (i) A particle is placed on one of the plane faces of 
a uniform circular cylinder at a small distance from the centre of that 

face ; prove that it will make small oscillations of period 2 ( * ) (a 2 -f h 2 )*, 

\i*pn>/ 

where a is the radius of the cylinder, h its height and p its density, and p 
is the attraction between two unit masses at unit distance apart. 

[M. T. 1889] 

The attraction of the cylinder at an external point on its axis is 
found from 2*52 to be 

2irw{h-V[(h + *) z + a 2 ]+V(* 2 + a*)} ............. (1), 

where z is the distance of the point from the centre of the nearer end. 
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Taking the centre of the plane face as origin and the axis of x 
through the particle and the mass of the particle as unity, its equation 
of motion is 



OTr 



where dV/dx is the force of attraction at (x, 0, 0). By MacLaurin's 
theorem 

(dV\ 



And at the origin ( -- ) = 0, by symmetry, and ( -^-r-) = ( -^-y ) . But 
\OX/Q \ox /o \oy /o 



where the expression (1) represents dV/dz at (0,0,2). Hence, on 
the axis 



1 

)/ ' 



, , ., . . 7rap 

and at the origin ^- i ^ = _^_ ........................ (4). 

Hence, from (2), (3) and (4), the equation of motion is 



giving oscillations of the period stated. 

(ii) A solid bounded by the plane xy and the half of the ellipsoid 
x 2 /a* -f i/ 2 /6 2 -f z 2 /c 2 = 1 on the positive side of the plane has a volume 
density ij(I/a 2 + l/6 2 + 3/c 2 ). Prove that it produces the same field of 
attraction outside itself as would a layer of surface density 



over the plane boundary and ^zjp over the curved surface, when p is the 
central perpendicular on the tangent plane. [C. 1928] 

To prove the equivalence of the two fields of attraction in external 
space, it is only necessary to shew that if we reverse the sign of one 
distribution then the two distributions together produce a zero field 
at external points.^ t 

The potential in this combined field may therefore be assumed to 
be zero outside the solid and some function inside the solid which 
vanishes over the surface. We therefore write 

V = Az ( 1 - o; 2 /a 2 - y 2 /^ 2 - * 2 /e 2 ), inside, 
and V = 0, outside. 
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Hence we get 

- A 



a 



6 2 ' 



ox* a* cy* t 

If p is the volume density, we have 



or, putting A = 2/xTr, p = ^z ( + r + - a ) inside, 

^ - b cl (i). 

and p = - T- V 2 F X = outside 

4?r J 

Again differentiating along the outward normal, on the plane 
boundary we have 



and, putting A 



Also, on the curved boundary 

dV 3V 



which reduces to 2Az/p, so that, putting A = 2/^Tr, we have 

a^-p,z[p (3). 

Since the distributions (1), (2), (3) give a zero field outside the body, 
it follows that when the sign of the surface distribution is changed it 
produces the same external field as the volume distribution. 

EXAMPLES 

1. If the potential in the interior of an attracting mass bounded by 
two concentric spherical surfaces be proportional to l/A/ r where r 
is the distance from the centre, find the law of density, and the potential 
exterior to the attracting matter. 

Shew that at any point of space where the density is p and potential F, 
the magnitude of the force is given by the expression 

[C. 1892] 



EXAMPLES 95 

2. Find the distribution of matter which produces a potential 
ft (# 2 -f y a 4- 2 2 2ax) at all points inside the sphere # 2 -f y* 4- z 2 2a# = 
and zero potential at all points outside it. fC. 1908] 

3. Find the distribution of matter, none of which is external to the 
sphere r = a, which produces a potential 2-rrpcc 2 at points within or 
on the sphere, the origin being the centre, x the abscissa, and p a 
constant. [C. 1914] 

4. What distribution of matter will produce potential 

Mf. ax\ 
r TV 1 -**)' 
according as r is less or greater than at ' [C. 1902] 

5. The potential outside an infinite elliptic cylinder 



is zero; inside it is V = x z /a* -K2/ 2 /6 2 1. Find the distribution of matter. 

[C. 1916] 

6. If the potential within an attracting medium having spherical 
symmetry is given by F = (c 2 -f r 2 )"*, where r is the distance from the 
centre and c a constant, prove that the density is given by 



and that the potential energy of the whole medium is 37r/32c. 

[M. T. 1915] 

7. Find the distribution of matter that gives rise to potentials 

2 - y* - z 2 ), r < a 



and ^2 = ^ + A^(2^ 2 ~2/ 2 -2; 2 ), r>a. 



8. Prove (by taking x* -f # 4 -f z* 1 as a potential function) that, if the 
space enclosed by the surface # 4 + y* -f z 4 = 1 is filled with attracting 
matter of density proportional to the square of the distance from the 
origin, its attraction at all outside points will be the same as if its mass 
wore distributed over the boundary surface with surface density 
inversely proportional to the distance of the tangent plane from the 
origin. [M. T. 1890] 

9. Shew that the family of surfaces defined by 

x 2 + y z = const. 

can be a family of equipoteiitial surfaces in free space, and find the 
law of potential. [C. 1906] 

10. Shew that the system of co -axial cylinders 



can form a system of equipotential surfaces. [London Univ. 1934] 
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11. Shew that the cylinders x* -f- y 2 = 2Arc are a possible set of equi- 
potential surfaces in empty space, but that the spheres x* + y 2 -f z 2 = 2Aa? 
are not. 

In the former case find the potential as a function of A. 

[London Univ.] 

1 2 . Two equal circular cylinders, with their axes in the same straight 
line and their centres of gravity at a great distance r apart, attract 
one another according to the law of the inverse square of the distance. 
Shew that the resultant attraction is approximately 



where M is the mass of one of the cylinders, C is its moment of inertia 
about its axis and A its moment of inertia about a perpendicular line 
through the centre of gravity. 

13. If two distributions of matter, A and B, have throughout 
a region T, to which both are external, the same equipotential surfaces, 
the attractions due to A and B throughout the region T are in the 
same direction and in a constant ratio. 

14. If two distributions of matter, A and B t have the same closed 
surface S, enclosing both, as an equipotential surface, every surface 
outside S which is an equipotential surface for A is also an equipotential 
surface for B. 



15. A finite solid of revolution is cut at right angles by a plane, 
and a particle is placed on the plane, a smooth surface, near to its 
centre O. Shew that the time of a small oscillation is 



where V is the potential of the solid, the axis of the solid is the axis 
of z and d*V/dz * denotes the value at O. 

16. A gravitating solid of revolution is cut by a plane perpendicular 
to the axis. A particle is fastened by a fine string of length I to a point 
in the prolongation of the axis, so that when the string is vertical the 
particle just does not touch the plane face at its centre O. Assuming 
the conditions such that when the particle is slightly disturbed the 
motion is that of a simple pendulum, shew that the time of a small 
oscillation is 



where E is the force exerted by the solid on unit mass at and R' 
is the derivative of that force at 0, taken outside the solid, along the 
axis. [C. 1892] 
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17. A uniform circular disc is of mass M and radius a and has its 
centre of gravity at O. A particle of mass m is situated at P, at a great 
distance r from O, and OP makes an angle with the axis of the disc. 
If the attractive force of the particle on the disc is reduced to a force 
through O and a couple, shew that the moment of the couple is 
($ywMa 2 sin 26) /r 3 , and find the magnitude of the force. 

[London Univ. 1938] 

18. Prove that the ^-component of attraction of a gravitating 
medium on a small sphere of itself with its centre at (x, y, z) is approxi- 
mately 



where a is the radius of the sphere, p the density and F the potential 
at (x,y,z). [C. 1932] 

19. Shew that, if the curves /(#, y, A) = form a system of equi- 
potential lines in freo space for a two-dimensional system, the surfaces 
formed by their revolution round the axis of x cannot bo a system of 
equipotential surfaces in free space unless 



ySy dx \9y 

is a constant or a function of A. [C. 1905] 

20. Shew that the family of cylinders 

(x* + / 2 ) 3 - 2a 3 (a; 3 - 3xy*) + a 6 = const. 

is a possible form of equipotential surfaces, and find the corresponding 
potential. [London Univ. 1938] 



21. Shew that (x-c) 2 + y 2 = X[(x + c)* + y 2 ] represents a family of 
equipotential surfaces, and that 



where V is the potential function and A and B are arbitrary constants. 

[London Univ. 1939] 

22. Determine a distribution of matter whose potential shall be 
constant over confocal elliptic cylinders. [C. 1915] 

23. Shew that the family of ellipsoids 

** y* 



where A is a variable parameter and a, 6, c are constants, is a possible 
form of equipotential surfaces, and express the potential in terms of 
a, 6, c and A. 
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24. The potential functions in the three regions bounded by the 
spheres r = 6 and r = a ( > 6) are 



r>a; 

a>r>6; 



and F 3 

Find the volume and surface distributions of matter and the total 
mass of the distribution. [London Univ. 1939] 

26. In the interior of a spherical distribution of matter of radius a 
the gravitational potential at a distance r from the centre is 

) , , 
+ const. 



Prove that the density at any point is ~~~ 2 - , - , and that the mass 

of the whole distribution is AajO, G being the constant of gravitation. 

Prove also that the negative potential energy (compared with a state 

of infinite diffusion) is %A 2 a/G. [M. T. 1925] 

26. Find the attracting systems whose potential is 
a 2 -7r a ) + 47rrxl(x z + y 2 ) for 



a 9 
and ^-nxyz ^ -f 47ra#/(# 2 + y 2 ) for r > o, 



= x 2 Hr y 2 + z*. [C. 1927] 

27/ The axds being rectangular, determine what volume and surface 
distributions ojf matter will give rise to potential Kxyz(a x y z) 9 
where K and # are constants, at all points within the tetrahedron 
bounded by the co-ordinate planes and the plane x + y -f- z a, and 
potential iero/at all points outside the tetrahedron. 

Shew that ^he total mass of the matter contained within the tetra- 
hedron isKc^/SOTry. [London Univ. 1938] 

28. An jnnnite uniform straight line occupying the positive part 
of the axis of z attracts according to the law of nature. Shew from 
elementary considerations that its potential is a function of r z, 
where r is distance from the origin, and hence prove that its potential is 

c log (r- z) + c'. 

29. Prove that, if the matter be distributed symmetrically about 
bhe axis Oz, d*V/dzdr~Q, at every point of thi^ axis, where r is the 
distance from this axis. Prove also that if the origin be a point of 
equilibrium, tfte lines of force near it are given by the equation 

zr a = const. [C. 1901] 
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30. Compare roughly the couples exerted on the* earth by the sun 
and moon respectively, supposing them to be at the same angular 
distance from the equator and the orbits of the moon and earth to be 
circles. Take the ratio of the mass of the moon to that of the earth 
to be 1/80. [M. T. 1922] 

31. Obtain the mutual potential energy of two gravitating bodies 
placed at a considerable distance from one another, neglecting powers 
of the inverse distance higher than the third, and obtain the specifica- 
tion of the moment, about the centre of mass of either body, of the 
attraction exercised upon that body by the other. 

Shew that the couples exercised by the sun and the moon on the 
earth neutralize one another at an instant when they are on the same 
meridian if 



sin2ai~" 

where a/ and o> are the angular distances of moon and sun from the 
equator, E, M and S are the masses of the earth, moon and sun, 
T' and T being the lengths of the sidereal month and year respectively. 
The orbits of the moon and earth may be considered circular. 

[M. T. 1937] 

32. O is the centre of a uniform cube of mass M and P is an external 
point. The distance OP is large compared with the edge a of the cube, 
and the direction cosines of OP referred to three concurrent edges of 
the cube are Z, m, n. Shew how to express the gravitational potential 
at P in a power series in OP" 1 , and evaluate the potential in terms of 
OP, M, I, m, n and a, rejecting all terms in OP -1 of higher order than 
the fifth. [P- 1936] 

33. Shew that the equipotential surfaces of a uniform equilateral 
triangular lamina (of side 2a), which are at a great distance from the 
lamina, approximate to the surfaces of revolution obtained by rotating 

the curves 

1 a 2 (l-3cos 2 0) , 

r + "l53 -- ' = const. 

about the initial line. [London Univ. 1938] 

34. Prove that, in the case of two equal cubes at a great distance 
apart, the third order term in the expression for their mutual gravita- 
tional potential energy disappears and obtain the next non- vanishing 
term. [M. T. 1929] 

35. Prove that for two gravitating solids, if one is a fixed sphere 
and the other has unequal principal moments of inertia and is free 
to turn about its centre of gravity, the only positions of equilibrium 
are those in which a principal axis at the centre of gravity of the 
second body is directed towards the centre of the sphere, [P. 1933] 

7-2 
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36. Prove that a distribution of matter of surface density 

over an ellipsoid will produce the same potential at every external 
point as a homogeneous concentric sphere of mass m within the 
ellipsoid, r and p being the central radius and the perpendicular from 
centre on the tangent plane. What is the potential of this distribution 
inside the ellipsoid ? [C. 1892] 

37. In a region of space which does not intersect the plane 2/ = 
a function <f> satisfies V 2 </> = 0, and the surfaces with the equations 
<f> = constant are portions of the elliptical cylinders 



where A is a variable parameter. By transforming to co-ordinates 
f , r), defined by 

x = a cosh f COST;, 2/ = asinhf sin??, z=f, 



prove that </> is of the form A + B, where A and B are constants. 

[P. 1937] 
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5*4. Green's equivalent layer. Let S be a closed equi- 
potential surface of a given distribution of attracting matter 
of which some or all lies inside S. For the sake of brevity 
we assume that the matter inside S can be represented by 
a volume density p. Let r denote distance from a point Q 
outside S. Apply Green's theorem to the region inside 8 in 
the form 



- (Vv 2 F'<fo= (v'^dS- (v'VWdv ...(1), 



where V denotes the potential of the whole distribution of 
matter and V I jr. 
This gives 

(V^da- fFV'lcto- ( I d8- fi W* ...(2). 
J dn J r J rdn J r v/ 

Let V 8 be the value of V on the equipotential surface S, 
then the first integral in (2) is V s -dS. But Ifr is the 

potential due to a particle of unit mass at Q, and Q is outside $, 

so that by Gauss's theorem ~-^-d8 = 0\ therefore the first 

J dn 

integral vanishes. The second integral also vanishes since 
V 2 - = throughout the range of integration. Also V 2 F = 4rrp 
throughout the same range, so that (2) reduces to 



But I - is the potential at Q due to the matter inside 8, so 
that denoting this by (V) Q we have 



Therefore the part of the potential at any point outside S 
which is due to the matter inside 8 is the same as would be 
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produced by a layer of matter of surface density 

1 3F 



spread over S, where V is the potential due to the whole of 
the attracting matter. 
Further, the amount of matter in the layer is given by 



(ads 

J 



, 5-2(1), 



= mass contained within S. 

We have therefore proved that the matter contained within 
any closed equipotential surface 8 in a given field can be 

spread over the surface S, with a surface density (dV/dn) 

at any point, without altering the potential at any point of 
the field outside 8. This distribution is known as Green's 
equivalent layer. 

If we suppose that this distribution over the surface S of 
the matter inside S has been effected, S is now an equipotential 
surface containing no matter, so that the potential at any 
point inside S is now V s the value on the surface [4-3 (ii)]. 

This can also be established analytically thus: let r' be the distance 
from a point Q' inside S and let p denote the density of the matter 
outside S. We apply Green's theorem in the form (1) to the region 
between S and the infinite sphere, taking V as before to be the potential 
of the whole distribution of matter and F'= 1/r'. The normal deriva- 
tives in (1) are now directed along inward normals to dS, so if dn is 
to retain the same meaning as in (3) and (4) we must change the signs 
of the surface integrals and write 



The first term is V 8 I *' ' dS, and since r' is measured from a point 
Q' inside S, ^ dS = - 47r, (Gauss), and the first term in (5) is 
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&irV g . The second term vanishes as before, and the fourth term is 

4 7r //L. 7 _. so that 

1 ' 



Or, if we put a for the surface density of Green's layer, as in (4), 



We thus verify that the potential inside S due to the joint effects 
of Green's layer and the original matter outside S is the constant 
potential of the equipotential surface. 

5'5. Dirichlet's problem.* There exists a solution of 
Laplace's equation which with its first and second derivatives 
is single-valued, finite and continuous throughout a region 
bounded by a surface S and takes arbitrarily assigned values over 
the boundary S. 

Dirichlet's method of investigation is as follows: let V be 
any function which satisfies all the conditions save that it is 
not necessarily a solution of Laplace's equation, then we can 
prove that, of all such functions, that one which gives a 
minimum value to the integral 



through the bounded region is a solution of Laplace's equation. 
Let V l be the function which gives a minimum value to W. 
If V be any function which satisfies the same conditions of 
continuity and takes the same values over S, then V = V l + V 2 , 
where V z = on 8. Conversely, let V 2 satisfy the same con- 
ditions of continuity, and let V 2 = on S. If is any constant, 
then TI 4- e! takes the assigned values on S and satisfies the 
conditions of continuity. Thus the set of functions to be con- 
sidered are comprised in the formula \ + eT, where V 2 = on 8. 
Then ri/TW r)F\ 2 \ 

+ * 



* Peter Gustav Lejeune Dirichlet (1806-1869). German mathe- 
matician. 
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Since W l is the minimum value of W, W W 1 must be 
positive whatever be the sign of e, so that the second term 
on the right of (1) must vanish. But by Green's theorem the 
coefficient of 2e 



and Kj = on $, so that the first integral vanishes and 
therefore 



f 



for all functions K> . We must therefore have V 2 V l = through- 
out the bounded region; for if V 2 V l had a sign throughout 
any portion of the region we could take K> to have the same 
sign and the integral would not vanish. 

The remainder of the argument is an assertion that among 
the different functions which satisfy the conditions of con- 
tinuity and take the assigned values over 8 there must be 
one which gives W a smaller value than the rest, and, from 
what precedes, this is a solution of Laplace's equation. But 
the argument is invalid, because Weierstrass has shewn that 
we are not entitled to assume the existence of a minimum, 
but only the existence of a lower limit to the possible values 
of W and we cannot assert that the limit is attained. 

Many mathematicians have produced solutions of Dirieh let's 
problem but they are too lengthy for reproduction here.* 

The theorem is also applicable to the region between the 
surface $ and the infinite sphere if we impose the condition 
that the required function shall vanish like l/r at infinity. 

We are entitled to conclude from the argument given above 
that if there is a solution it is the only one. For the above 
argument shews that, if V 2 ^ and 1^ on $, then the 
coefficient of in (1) is zero, so that IF == B^ + 2 PK>. Similarly 
if V 2 K = 0, then since Ji= V dK,, we have, by the same 
argument, "H = PF 4- e 2 W. Hence e^W^^O, and so eTg =*=(), i.e." 
V = \l and the solution is unique. 

* See H. Poincare, Theoric du Potential Ncwtonien^ and, for a 
bibliography of the subject, P. Appell, Traiie de Mecanigue ffi&tiQnelle, 
t. in, p. 98. 



5*6] GIVETST BOUNDARY VALUES 113 

5'6. To find the solution of V 2 F inside a region bounded 
by a closed surface S when the values of V and dV/dn over 8 
are known. 

Take Green's theorem in the form 



Put F' = l/r, where r is the radius 
vector drawn from a point P inside S. 
Take a sphere H of small radius c 
with P as centre and apply the 
theorem ( 1 ) to the region bounded by 
8 and 2. Then, since V 2 (1/r) - and 
V 2 F = throughout the region, the 
right-hand side of (1) is zero, and 
therefore 




f 

J 

Taking account of the inward direction of dn over S we see 
that d (l/r)/dn 1/e 2 , and since we may put eZS = e 2 dco, we have 

V~~- r/S = Vdco, and this being independent of we may 
J 2 &% J 



Again 
reduces to 



take as small as we please and write for the integral 
where V P denotes the value of F at P. 

r i d v r i d v 

in, - = -TT 2 ao> -> as c -> 0; so that (2) 
J 2 r on J % e dn 

S (3), 



F P = 



giving the value of F at any point inside /S in terms of the 
values of F and dV fin on the surface. 

When the region is bounded by the surface 8 arid the infinite 
sphere, it is easy to shew that the same theorem holds good 
provided we impose the condition that RV and R 2 dV/dn 
remain finite as R tends to infinity, where R denotes distance 
from an origin in the neighbourhood of $. 
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5*61. Green's function. We learn from Dirichlet's problem 
that there should be a unique solution V of Laplace's equation 
which with its first and second derivatives is finite through 
a bounded region and takes a prescribed form over the 
boundary, but the procedure of 5' 6 seems to imply that we 
need to know not only the value of the function V but also 
the value of dV/dn over the boundary in order to determine V 
at points inside the region. The problem can however be 
reduced if we can find a function of a more restricted form 
known as Green's function and defined as follows: 

Let if be a solution of Laplace's equation inside the boun- 
dary $, which takes the value l/r on S, where r denotes the 
distance from an assigned point P inside S. 

Let G = H + l/r; then is called Green's function for the 
point P and surface S; and by definition we see that G vanishes 
on 8 and satisfies Laplace's equation at every point inside 8 
except at the point P where it becomes infinite. 

In 5-6 (1) put V' = H, then since V 2 F = by hypothesis and 
= 0, we have 



adding this, with a factor 1/47T, to the right-hand side of 
5-6 (3), we get 



But G vanishes on S so that 



>=-N> 

477" J 



By this means Dirichlet's problem is reduced to the deter- 
mination of Green's function. 

The case where the region is bounded by a closed surface 8 
and the infinite sphere does not need special consideration 
if we impose like conditions as before to the way in which the 
functions vanish at infinity. 



5-63] GREEN'S FUNCTION 115 

5*62. Green's function for a sphere. Let be the centre of^a 
sphere S of radius a and P an internal point at a distance / from O. 
Let P' be the inverse of P, so that OP. OP' = a 2 . Let r, r' denote the 
distances of any point M from P and P'; and for a point M f on the 
sphere let PM / = r 1 and 



rriu i_ i j. i T \ OM' a 
Then by similar triangles ~ = -^-p = > , 



so that 




Now consider the function H= a/fr'. It is a potential function 
due to a mass ( a/f) placed at P', and therefore it satisfies Laplace's 
eqiiation inside the sphere; and because of (1) it takes the form 1/r 
on S. 

Hence O = j, (2) 

^ / 

is the required Green's function for the sphere S and point P. 

The same formula can be shewn to hold good when the point P 
lies outside the sphere at a distance / from the centre, P' then being 
inside and the symbols having the same meanings. 



5'63. Reciprocal property of Green's function. // P, Q 

are two points inside a region bounded 
by a surface S and G (P, Q) denotes 
the value at Q of Green's function for 
the point P and surface S, then 



Apply Green's theorem in the 
form 5-6 (1) to the region bounded 
externally by S and internally by 
two spheres S 1? S 2 of small radii 
e 1} 2 having their centres at P, Q. 




8-2 
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Put F= = ~+#, and V' = G'=^+H', where r, r' denote 
r r 

distances from P and Q and G, 0' denote Green's functions 
for P and 8 and for Q and S respectively. Then we have 



By hypothesis H, H' are solutions of Laplace's equation, 
and since P, Q have been excluded from the region of volume 

integration V 2 ~ = and V 2 ~ = 0, so that V 2 (? = and V 2 (?' = 

and the right-hand side of (1 ) is zero. 

Again, by definition, G and G' vanish on 8, so that the 
first integral on the left is zero. 

In the second integral on the left we put ^2^ = e^dw so that 



which clearly -> as e -> 0, since H is assumed to be finite; and 

f ^ 3G iy (V 'l l ^ dH \ ^ 
- G ^~dL l = - \Gp \~i+^~ 2 ao>, 

Js, 3^ J \ 2 dn) 

and as e -> this reduces to - 4TrG P ', where Gp denotes the 
value at,P of Green's function for Q and 8, i.e. G (Q, P). 

In the same way the third integral in (1) may be shewn to 
reduce to ^GQ or knG (P, Q); so that (1) reduces to 

G(P,Q) = G(Q,P) .................. (2). 



The theorem is also true for a region bounded internally 
by a closed surface 8, and externally by the infinite sphere, 
assuming that H and H f are at least of order R~ l at a great 
distance jR. 
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1. Use Green's theorem to prove that, if F, V are solutions* of 
Laplace's equation such that F = V at all points of a closed surface S, 
then F= F' throughout the interior of S. 

2. Prove that, if F, V are solutions of Laplace's equation such that 
6V/dn = dV'/dn at all points of a closed surface S, then F-F' is 
constant at all points inside S. 

3. Shew that, if V lt F 2 are the potentials of two equipotential 
surfaces which surround bodies of mass M which produce a gravita- 
tional field, the contribution of the space between these equipotential 
surfaces to the total energy of the field is \M \V F a ). 

4. Prove that there cannot be a function F which satisfies Laplace's 
equation inside a closed surface S and makes dF/0n= kV at all 
points of S, where dn denotes the normal drawn towards S from the 
inside and A; is a positive number. 

5. A number of bodies are formed of gravitating matter. Prove 
that the resultant force, resolved parallel to the axis of x, due to the 

attraction on one of the bodies, of uniform density p, is p I nx x) NdS; 

where the integral is taken over the surface of the body, x is the abscissa 
of the centre of gravity of the body, and N is the normal attractive 
force at each point of the surface. [C. 1890] 

6. Prove that, if dm, dm' denote elements of two different distribu- 
tions of matter which have the same external level surfaces, and <f> 
any function which satisfies Laplace's equation, then 



Udm 



where M , M' denote the total masses of the two distributions. 

7. Prove that, if two distributions of matter have the same external 
level surfaces, they have the same centre of gravity and the same 
principal axes of inertia at the centre of gravity. 

[In Ex. 6 put </> = x, and <j>~yz.~\ 

8. Prove that, if F is the potential of a given distribution of matter, 
S a surface surrounding it, r the distance of an internal point from 
an element dS, dn an element of the outward drawn normal and 6 the 
angle between r and dn, then 



where the integration is over the whole surface. 
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9. Shew that, if V is a function which with its first derivative is 
finite and continuous throughout a region bounded by a surface S 
and r denotes distance from a point P within this region, then 



= fvg-(l 
J dn\r 



~ 

r dn 

where dn denotes the outward drawn normal and the integrations are 
over the surface and through the volume. 

10. Two gravitating solids M and M' have potentials V and F'; 
the inward normal force across any element dS of the surface of M 
is the sum of N and R due separately to M and M' ; and the inward 
normal force across any element dS' of the surface of M' is the sum 
of N' and 22', due separately to M and M' '. Shew that the potential 
energy of their relative position /differs by a constant from either 

j- ((RV-NV')dS or ~ ((R'V'-N'V)dS'. [M. T. 1897] 

11. If V 2 F 1 = 0, V 2 F = and F 1 = F = Fat a closed surface, while 
F > F! and their differential coefficients are everywhere finite without 
and within the surface respectively, and F vanishes at infinity, find 
the potential in all space of the surface distribution or whore 



and dn is outwards. 

Shew that IV ~-\-]dS is equal to / ^dS when r is, 
J dn \r/ n f r dn 

* [I dV 

measured from an external point, and to I - dS when r is measured 

from an internal point. [C. 1908] 



12. Prove that, if < denotes any function which with its first 
derivative is continuous, its mean value over a sphere of radius a 
exceeds its value at the centre of the sphere by 



integrated through the volume of the sphere, where r denotes distance 
from the centre. [Minchin.] 

v 13. Assuming Gauss's theorem and Green's theorem and that a 
distribution of matter consisting of volume density />, surface density a 
and separate particles has a potential function F which is infinite at 
each particle and satisfies the relations V 2 F= 47r/> wherever there 

is volume density, and ( -A ) ( ~ ) = ^ira across each surface layer 

Van/! \tftt / 2 
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and is otherwise continuous, prove that the potential must be of the 
form 

F = E + [? ds 
r J r 

where m denotes the mass of a particle. 

14. Prove the theorem of the equivalent layer [54] by applying 
Green's theorem to the region bounded externally by the infinite 
sphere and internally by the equipotential surface S and a small sphere 
round the point Q, and excluding the singularities of the potential 
function. 

15. Prove that, if M be the mass of a given distribution of matter, 
U the potential of a closed equipotential surface entirely surrounding 
it, and E the exhaustion of potential ejiergy in collecting*the matter 
from a state of diffusion at infinity into its actual distribution, then 
\MU<E. [M.T.I 895] 

16. A function V vanishes at all points outside a surface S. Shew 

that a distribution of matter p= V 2 F within S and 0-= 4 A upon 

r ^ 



S will have the function F for its potential. Here dV/dn means the 
rate of increase of V in the direction approaching the surface from the 
inside. [C, 1931] 

17. Prove that, if F is the potential at points outside a surface S 
clue to masses within the surface, and V l is the potential at points 
inside the surface due to masses outside it, the expression 



1 {(V V\ 8 ( l \J? l /Y 8F 
4^ ] ( F ~ Fl) 8n \ r ) d<S " 4* J \8n ~ 



has the value F at external points and the value V l at internal points, 
where r denotes distance from the point and dn an element of the 
outward normal to dS. [M. T. 1910] 

18. F is the potential and R the attraction at any point due to 
a given distribution of matter none of which lies between two level 
surfaces F = V l , F = F 2 . The space between these surfaces is filled 
with matter of density R*f ( V) . Prove that its potential at any external 

/T 2 
point P is 4t7rU I /(F) dV, where U denotes the potential at P due 

J Vi 

to that part of the original distribution which lies inside the inner of 
the two level surfaces. 

19. If F is the gravitational potential of a distribution of matter 
inside a closed surface #, shew that the value of F at a point P outside S 
is given in terms of its values on S by the formula 
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where the differentiation is along the outward normal to S, O vanishes 
on S, and, if r is distance measured from P, Gl/r is a harmonic 
function on and outside S. Give a physical interpretation of O and 



Shew that, if V is the potential of a finite distribution of matter on 
one side of an infinite plane, and P is a point on the other side, 



where d& is the solid angle subtended at P by an element of the piano 
at which the value of V is V l . . [M. T. 1934] 

20. Use 5*61, 562 to shew that, if V be a function whose value is 
prescribed over a sphere S of radius a, and P a point inside the sphere 
at a distance / from the centre, then a solution of Laplace's equation 
inside S is given by 

a*- 
y P j~ 



where r denotes the distance of dS from P. 

21. If < t and </>$ are any scalar functions whoso derivatives exist 
and are continuous at all points of a volume V and its boundary surface 
S, and if 



prove that , Ai^a^Aa 8 - 

If fa is a given harmonic function in V and if < 2 = <f> on S y shew that 



Hence shew that for such a <f> z the necessary and sufficient condition 
that / 22 should be a minimum is that <f> z should be harmonic in V. 
What is the physical interpretation of this result? [M. T. 1935] 



Chapter VI 
HARMONIC FUNCTIONS 

6-1. Any solution of Laplace's equation 

927 3 2 F 8 2 F A 

__ I ___ I __ .__ Q 

dx 2 dy z dz 2 

which is homogeneous in x, y, z is called a harmonic function 
or a spherical harmonic. The degree of homogeneity is called 
the degree of the function. We shall be concerned with the 
case in which the degree is an integer. 

The study of harmonic functions is an important branch 
of mathematical analysis. We shall only concern ourselves 
with so much of the theory as is useful for our immediate 
purpose and be content with arguments based upon physical 
ideas rather than attempt a rigorous treatment of the subject. 

6-11. If V is a harmonic function of degree, n, then 



is a harmonic function of degree np q t. 

For if we differentiate the equation V 2 F = p times with 
regard to a?, q times with regard id y and t times with regard 
to z, we get 



6-12. Surface and solid harmonics. In polar co-ordinates 
Laplace's equation may be written 

'-O [1.68,,)] 



(1). 
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If we substitute V = r n S n , where S n is independent of r 
we get 2-(r*-^-\~n(n + l)r n S n , so that (1) reduces to 



or, if cos = /x, 



A solution 8 n of equation (2) is called a Laplace's function 
or a surface harmonic of order n. Since n(n + I) remains 
unchanged when we write (n -f 1) for n, there are two solu- 
tions of (1) of which 8 n is a factor, namely r n S n and r~~ n ~ l S n , 
and these are called solid harmonics of degree n and (n + 1 ) 
respectively. 

6*1 3. It follows that if U is a harmonic function of degree n, 
then U/r 2n+l is also a harmonic function. For we may write 

TJ Sf 
U = r 71 /^ , so that -^TJ = ~^TI > which is a harmonic function. 

For example, it is obvious that xyz is a solution of Laplace's 
equation of the third degree; therefore xyz/r 7 is also a solution. 

Similarly, if U is a harmonic function of degree (n + l), 
then r 2n+1 U is also a harmonic function. For we may write 
U = r~ nLl S n) so that r* n + l U = r n S n , which is harmonic. 

These theorems combined with 6-11 make it easy to build 
up a great variety of harmonic functions. 

Thus x or 1 and 9 or tan" 1 (y\x) are both harmonic functions 

1 1 
- and - 

ii X 

as harmonics of degree 1 . Differentiating these with regard 

1* "?/ 2 2 77 

to x, y or z. we find that -5 , ~ , -^ and -~ tan" 1 - are harmonics 
y r 3 r 3 r 3 r 3 x 

of degree 2, and so on. 

Again, if dh denotes an element of length whose direction 

fl fl 7) 7) 

cosines are Z, m, n, so that ^ r = Z-^-+m-^-H-n7r>it follows that 

<m ox oy oz 



1 1 V 

of degree zero. If we divide them by r , we get - and - tan" 1 - 
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*\ / 1 \ 

5TJ" I is a harmonic function and h may be called its axis. 
en \r / 

Similarly, if h l9 A 2 , A 3 , ... are any number of arbitrary direc- 

3 a a 



tions, we can form a harmonic function 

\r 

6'2. Legendre functions.* Let M be a fixed point (/, g, h) 
and P a variable point (x, y, z), then if MP = J?, we know that 
the inverse distance IjR is a solution of Laplace's equation, 
for it is the potential at P due to a particle of unit mass at M . 
Suppose that M lies on a sphere of radius a with its centre 
at the origin 0, so that/ 2 + 2 + A 2 = a 2 , and let 0P = r. Then 
we have 

y = JL 



Since # 2 + ?/ 2 + z 2 = r 2 , the expression for F can be expanded 
in powers of 1/r, by the extension of Taylor's theorem to 
functions of several variables, thus 



dy 



n ! 



the series converging if r > a, i.e. if the point P is outside the 
sphere. But, if we write JJL for the cosine of the angle MOP, 

we have 

y _. 



and this can be expanded in powers of r/a when r < a, and in 
powers of a/r when r > a. Thus if we denote the coefficient of 
r n /a n+1 or of a n /r n + l by P n) we have 



and 



* Adrien Marie Legendre (1752-1833). French mathematician of 
great distinction. 
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The P's are clearly polynomials in p,, for they are the 
coefficients of successive powers of h in the expansion of 



or 

..V'|;;^ 



...... (5). 



By picking out coefficients of successive powers of h we 
find that 



and 



, ">\'b l)\U*)\H, &) n _ 4 ^ 

2.4(2n-l)(2n-3) * 

the series ending with a term in fju or a constant according as 
n is odd or even. For many purposes it is sufficient to know 
that P n is a polynomial in ju, of the form 



2 ^~* + .................. (7), 

where A n , A n _ 2 , ... are constants. 

The function P n is called the Legendre function or Legendre 
coefficient of order n. It is written P n (p,) when it is desired 
to specify the variable. We observe that //, = cos#, where 6 
is an angle measured from an axis M , called the axis of the 
function. 

In particular, for 9 = or TT we have ^ 1, and 

1 
P n (1) = coefficient of h n in - ^ = 1, 

'; 

and P-l = ........................ = ~ 1 



If we write ju. = cos 9, then 

)-* = ( 1 - he* 6 )-* ( I - he- if> )-* 



...V 
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It follows that P is of the form 



n 



where n , B n _ 2 , . . . are positive. Hence the numerical value of 
P n does not exceed the number obtained by putting 6 = in 
the above expression. 

Thus |P n O*)|<P w (l);i.e.|P n O*)|<l (9). 

6*21. r n P n and r~ n ~ 1 P n are harmonic functions. Since 
ju, = cos MOP = (fa + gy + hz)/ar, from 6*2 (7) we have 

rP=^(h 

n a n \J 

A, 



expressing r n P n as a homogeneous rational integral function 
of the co-ordinates of P. 

Further, we observe that in 62, (2) and (4) represent the 
same expansion of V, so that by comparing the terms in r~ n ~ l 
we have 

n ~ l d d i 9 

+ff 



But l/r is a harmonic function, therefore by 6-11 every 
term on the right-hand side of (2) is a harmonic function of 
degree (n + l). Hence r~ n ~ l P n is a harmonic function, and 
by 6-13 r n P n is also a harmonic function. 

622. Rodrigues's formula for P n .* To prove that 

1 d n 



This formula may be verified by expanding (^ l) n by the binomial 
theorem and differentiating the result n times; or we may proceed 
as follows: 



Let l-hu=(l- 2^ + ^)4 ..................... (1), 

so that u = p. + %h(u 2 1). 

* Corresp. sur V 'A 'cole poly technique, vol. in (1814-16), pp. 361-86. 
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By a theorem of Lagrange* u can be expanded in powers of h in the 

form n , n-l 



But by differentiating (1) we have 



~ - =:SfcP n .................. (3). 

Then by differentiating (2) and comparing the coefficients of h n in (3) 
and in the result of differentiating (2), we obtain the required formula. 

6-23. The equation P n (/*) = 0. Since the equation (ju, 2 l) n = has n 

roots equal to 1 and n roots equal to 1, the equation -^- (/u, 2 l) n = 

djji 

has n1 roots equal to 1, n 1 roots equal to 1, and a root f be- 

ef 2 
tween 1 and 1. Then -=-5 (/x 2 l) w = has n 2 roots equal to 1, 

din 

n 2 roots equal to 1, a root between 1 and g and a root between 
f and 1. Continuing the process it follows that the n roots of the 
equation P n (/x) = are all real and lie between 1. 



6-24. Differential equation for P n . Since r n i^ is a har- 
monic function and P n is a function of 6 only, r n P n satisfies 
Laplace's equation in the form 



and as in 6*12 we have 



or putting cos 9 = /z, 



This is commonly known as Legendre's equation. It can be 
solved in series in the form 6*2 (6). 

6*25. Integral theorems. To prove that 



* Whit taker and Watson, Modern Analysis, 7*32. 
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From 6-24 we have 



and 
Thetefore 



and, integrating by parts, this 



r 1 

Again, to find P n 2 (p,) dp,, we have 
J-i 



Integrate between limits 1 to 1 and use the fact that the 
integrals of all products P n P m (n=m) vanish, and we get 



J 



-1 

1, l-h 



s 6-26. From the expression for P n (^) in powers of ^, it is clear that 
by solving a set of linear equations any positive integral power of p 
can be expressed in terms of a finite number of P's in the form 

/ i w = a P n -fa 2 P n _ 2 -fa 4 P n _ 4 +... 
and consequently, from 6-25, if n < m, 
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And similarly, if F (p.) is a polynomial of degree less than ra, 



6*261. If a polynomial in ^ is expressed in terms of Legendre 
coefficients, say in the form 



the coefficient A n can be found by multiplying both sides by P n (fi) 
and integrating between limits 1 and 1 . For then by 625 all the terms 
on the right except the term in A n vanish, and 

2n+l 



6-27. Recurrence formulae. These are a set of formulae con- 
necting successive P's or their derivatives. 
We begin with the formula 

(l-2/^ + fc a )-*=l + S&"P n (,K) ............... (1). 

Differentiation with regard to h gives 

l P n (iJL) ........... (2).* 



Multiplying by (1- 2/x/t + A 2 ), we get 



(/t - h) l + hP n (/.) = ( 1 

then equating the coefficients of h n gives 

(nH-l)P n+1 ~(2n+l) / *P n + nP n _ 1 = ............ [A]. 

Again, differentiate ( 1 ) with regard to /x and we have 



- , 
i dp, 

and therefore, from (2), 

00 oo /7P 

^ L nh~^P n = (iJL 

1 1 

Whence, by equating coefficients of 7z. n 



Again, by differentiating [A] we get 
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and eliminating dP n /dp. from this result and [B] gives 

(2n+l)P n -?^-^^ .................. [C]. 

a\L ap 

Subtracting [B] from [C] gives 



In [D] change n + 1 into n and subtract p times [B] and we get 

(^-l) d f- n = n,.P n -nP n ^ .................. [E]. 

Integrating [C] gives 

fp n (ridp = ^~^ + const ................ [F]. 

Lastly, by adding successive formulae of the form [C] we get 



the series ending with 3P! if n is even, and with P or 1 if n is odd 

...... [0]. 

6-271. Example. Use the theorems of 6-25 to prove 6-27 [G]. 

Since P n is a polynomial of the nth degree in ^, therefore dP n /dp, 
is a polynomial of the (n 1 )th degree and by solving linear equations 
is expressible in the form 

dP * 

jjf = 0.-1 Pn-l + 0-s Pn- 3 + 

To find any coefficient, say a mt multiply by P m and integrate 
between limits 1 and 1 ; then 



since by 6*25 every other term vanishes, 

., by 6-25 (1). 



But since m<n, dP^jd^ is of the form c m _ 1 P m _ 1 4-c m _ 8 P m _ 3 + . 
where the suffixes are all less than n. Therefore 



By adding (1) and (2) we get 
2m + 1~~ /-] 



since n m is odd, 6*2 (8). 

R N A 
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Therefore a m = 2m 4- 1, and 



6-28. Definite integral forms for P n . There are two ex- 
pressions for P n as a definite integral, due to Laplace, * which 
may be obtained as follows: we have 

d(j> 77 

/o " 



so that, by putting a=lp,h and b = h VC^ 2 1), we get 
1 ! 



V( l- 

_ i r* __ d<f> 



__ ___ 

Tro 1 fJihh VV 2 l)cos</>" 
Expanding both expressions in powers of A, we get 

S A n ^ = - f " S A n {/* + V(/^ 2 - 1 ) cos cf>} n dJ>, 

'"'JO 

and, by comparing the coefficients of h n , 



(1). 



\/(/>t 2 1) is imaginary for values of p, between 1 and 1, 
but if the binomial expression is expanded the integrals of all 
odd powers of cos^> vanish and a real result is obtained. 

By writing 



cos (f) = {fji cos i/j V(^ 2 ~" 1 )//* V^ 2 1 ) cos 
the result may be transformed to 



...... (2). 



6-29. After this lengthy discussion of the properties of the 
Legendre polynomials we revert to their introduction in 6-2 
as a special form of surface harmonic, namely a solution of 

* Mecanique Celeste, liv. xi, ch. n, p. 33 [1823]. 
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6*12 (2) independent of <f>. But we must point out that, if 
in 6-2 we take M as the point (a, 9', </>') and P as (r, 6, </>), then 



(1) 

(2), 

where A = u// + ^/(\-u?) V(1~M-' 2 ) cos W -('), 1 t , * 

r ...(3). 

This as before is expansible in powers of r/a or a/r, and it 
may be shewn, as in 6-22, that the coefficient of r n /a n+l or of 

a nj r n+i j s _^ {(A 2 l) n }, where A is given by (3). 



Expressed in this way the Legendre function P n (X) is 
a symmetrical function of 0, <f> and 6', </>', i.e. of the angular 
co-ordinates of P and M 9 and it is only when the polar axis 
of co-ordinates coincides with OM or OP that we have A = /x 
or A = IJL' and the Legendre function takes the form of a poly- 
nomial independent of </ and </>'. 

6-3. Applications of the Legendre functions. Since the 
Legendre function P n (p,) is a function of p or cos 9 only and 
r n P n and r~ n ~ l P n are solutions of Laplace's equation, it is 
evident that these functions may prove to be of use in connec- 
tion with the potentials of bodies with an axis of symmetry. 
In this connection there is a theorem of Legendre: // the 
potential of a body with an axis of symmetry is known at points 
on a finite length of the axis outside the body it can be calculated 
in terms of the Legendre functions at all points that can be reached 
from this length of the axis without traversing matter. It is to be 
understood that there is symmetry both as regards the density 
. and the shape of the body. 

Taking the axis of symmetry as the axis Oz, let Q be a point 
on the axis outside the body at a distance z from the origin. 
And suppose that the potential at Q is expressible in a series 
in the form 



9-2 
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and that this series converges for the z's pertaining to a certain 
length of the axis. 

Also let V 'denote the potential at a point P, where OP = r 
and QOP = 9; and let 



...... (2). 

Then V is a harmonic function [6*21], so that V V is also 
a harmonic function. Now let P move up to Q, then V becomes 
U and V becomes the series on the right of (1) (since at Q 
r = z and = 0), so that at Q, and other points on Oz for which 
(1) converges, VV vanishes. Therefore by 4*3 (v) the 
potential function VV vanishes at all points that can be 
reached in passing from Q and other neighbouring points on 
the axis without passing through the attracting matter or 
through a region in which (1) ceases to converge. We qonclude 
therefore that the potential at any such point PTm given by 



...... (3), 

and this series can clearly be written down, whenever (1) 
has been calculated, by writing r for z and introducing the 
appropriate P in each term. 

It is clear that (3) is absolutely convergent for all values 
of r which make (1) absolutely convergent when r is written 
for z, since | P n \ ^ 1. 

6-31. Potential of a circular ring. ,The potential of a 
circular ring of radius a and line density m at a point on its 
a,xis at a distance z from its centre is 27ram/(a 2 + z 2 )*. This can 
be expanded in powers of z/a or a/z according as z is less or 
greater than a, so that on the axis the potential is given by 

either z z 

z<a, 






la 3 1.3a 5 
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Therefore the potential at any point (r, 6) is either 

\ 

4 -...) ) r<a, 

l.Sa* 

2---P + -..., r>a. 



632. Potential of a circular disc. With the same nota- 
tion the potential at a point on the positive axis of z is, by 
2-51, 27Tm{'\/(z 2 +a 2 ) z}, where m is the mass of unit area. 
This can be expanded in the forms 

lz 2 l.lz 4 1.1.3s 6 
-s+j-- _+_- 

a 2 1.1 a 4 1.1. 3 a 6 



Therefor^ the potential at (r, 0) is given by 
1 r 2 1 1 r 4 1 . 1 . 3 r 6 



and 

a 2 l.la 4 1.1. 



r>a ...... (2). 

If r > a, the expression (2) is valid for ^ 6 ^ TT, since we can 
move from the point (r, 0) to the point (r, 6) keeping r > a 
without crossing the disc. 

If r < a, and > |TT, we cannot move from (r, 0) to (r, 0) 
keeping r < a, without crossing the disc- Thus (1) only holds if 



<d^7T, then, by symmetry, V^ (r, 6) = V l (r, TT - 0). 
Now P^TT- -0)= -^(0), and P 2 , P 4 , ... are unaltered if 
is replaced by TT 0. 

Hence when r < a and |TT < ^ TT 



Thus, as in 2*51, the form of the potential is discontinuous 
in crossing the disc. 
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6*33. Example. Find, in the form of an infinite series, the potential 
due to a flat circular annulus (e.g. Saturn's ring] at points moderately 
near the centre of the ring. 

If a small but massive planet occupies the centre of the ring, shew that 
the maximum deflection of the plumb line on its surface (on account of the 
attraction of the ring) will occur in latitude 45, approximately. 

[M. T. 1916] 

Let a, b be the inner and outer radii of the annulus. Its potential 
is the difference between the potentials of discs of radii a and 6. Taking 
the axis of the annulus as axis of z with the centre as origin, the potential 
at the point (0, 0, z) is 



where m is the mass of unit area of the annulus. 
For z < a, this gives 



Therefore, at points near enough to the centre for higher powers of 
r/a to be neglected, the potential is 



If the planet be of radius c, the deflecting force on the plumb line 
at a point on its surface at an angular distance 6 from Oz is dV/cdd, 
and if c is small compared with a, this is 

2<rrmb c 2 dl\ 2rrma c 2 dP 2 /I INdA 

, __ _ __ __ _ _______ _ _ _ __ * I _ I _ - 

~ 



, __ _ __ __ _ _______ _ _ _ __ Tr _ _ - 

c 2b* d~e c '2a 2 dO - W7/ * c/ ^& a ) dO ' 
But P 2 = (3 cos 2 6 I), so that the deflecting force is 

\Tfmc ( T j sin 26, 
and has its greatest value when 0= 45. 

6-4. Harmonic functions in general. We now turn back 
to the subject of harmonic functions m general. From 6-1 
H n will denote a harmonic function of degree n if it is homo- 
geneous of degree n in x, y, z and satisfies Laplace's equation. 
If it is a rational integral function of degree n in x, y, z, it 
contains %(n + l)(n + 2) coefficients; but since V 2 H n = Q, and 
this relation is of degree n 2, there are \n(n\) relations 
between the coefficients. Therefore \ (n + l)(n + 2) \n (n 1), 
or 2/1 + 1, of the coefficients are independent, so that there 
are 2n + 1 independent rational integral harmonic functions 
of degree n. 
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6*41. Integral theorem for surface harmonics. To prove 
that \lS m S n da> = Q, where S m> S n are surface harmonics of 

different orders, doj is an elementary solid angle, and the integral 
extends to the whole surface of the unit sphere. 

Apply Green's theorem 5*1 (1) to the region bounded by 
a sphere of radius a with its centre at the origin, putting 
V = r m S m and V = r n S n , so that V, V 1 are harmonic. It follows 

that 



when the integration is over the sphere r = a, so that dS = a 2 da). 
The relation reduces to 

(n - m) a m + n + l S m S n dco = 0, 
and, since n 4= m, therefore 

m S n dco = (1). 



Particular cases: (i) Since P n is also a surface harmonic, 

therefore 

(2). 



(ii) Since P = 1, therefore 



.................. (3). 



(iii) Since P m , P n may be taken as special cases of S m , S n , 
and we may write da> = sin#d0d</= d/zd^, the integral over 



[1 f2ir 

the sphere in tljis case becomes P m P n dp,dc/) J or 

J~iJo 

r 1 

2?r P m P H rfjLt, so that (1) reduces in this case to the theorem 
J-i 

, f P m P n dfjL = Q, 



of 6-25, 

6*42. Tesseral harmonics. Reverting to the equation 
6-12 (2) for surface harmonics, viz. 
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we proceed to examine solutions of the form S n = 0O, where 
is a function of 6 alone and O a function of < alone. On multi- 
plying the equation by sin 2 0/0O, it becomes 



We must therefore have 



and 

The latter equation has single-valued solutions when 
c=2> 2 , viz. 



Bin^ ............... (3); 

and is then given in terms of 6 by 



or, in terms of p,, by 



For ^9 = this is Legendre's equation with a solution P n (^), 
where 



Differentiate (5) ^p times with regard to /* and multiply the 
result by (lp*)**. This gives 



fl.nP 
Let z = (l-u 2 )^ 
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Therefore 
d 



By subtracting (6) from the right-hand side of (7), we get 



= -*(+ 1)- 

so that ^ is a solution* of equation (4). It is convenient to 
denote solutions of this form by X n p and, using Rodrigues's 
formula for P n , write 



It is clear that X n p is a rational integral function of sin and 
cos 0, and that it vanishes if n +p> 2n, i.e. if p > n. Also since 
is a surface harmonic, therefore 



X o x lc B cf> X 2 >c08 7)rA X ncOB n(f> 

^n 9 -A-n sinr? -^n sin^r* ^n sin ^r 

are 2n+l surface liarmonics. They are called Tesseral har- 
monics. That they are linearly independent follows from the 
fact that if we had a relation 

n l cos(f>+ ... +A n X n n cosn<f> 



it would be sufficient to multiply by p<f> and integrate from 
to 277, to shew that A p = 0, B p = for p = 1, 2, ... n\ and then 
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6*43. H n expressed in tesseral harmonics. From 6'4, every 
surface harmonic of order n can be expressed as a linear com- 
bination of 2n+l linearly independent surface harmonics, 
so that 



X n "smn<f> ...... (1), 

where the A's and B's are constants. 

Also since H n = r n S n , therefore H n is expressible in terms of 
harmonics of the form r n X n p ^p</), and it is easy to see that 
these are polynomials in x 9 y, z. 

Thus 6-42(8) may be written 



and taking for P n its expansion in powers of p, [6-2 (6)] we see 
that X n v is the sum of terms of the form sin*> (cos 6} n - p ~ 2k ; and 

r n smp<f> sin*> 9 (cos 0)"-*- 2 * = r 2k (r cos e) n - p " 2k r? sin*> sinp^. 

Here the first two factors are polynomials in x, y, z\ and for 
the last factor, if p is odd say and =25+1, 

sinpcf) = sin </> {sum of terms of the form (sin <) 2s ~ 2m }, 
and 
r p sin p 8 smp(f> = r sin'# sin <f> (sum of terms of the form 

(r sin 0) 2 (r sin 9 sin <) 2 *~ 2 } 
- y {sum of terms of the form (z 2 + y 2 ) y 2s ' 2m }. 

Similarly if p is even; and thus r n X n p sinp</> is expressed as 
a polynomial in x, y, z, and H n is expressed in terms of 2n + 1 
linearly independent polynomials. 

6*44. The integral S n P n dw. To prove that, if 8 n be any 

surface harmonic of order n of the angular co-ordinates 0, </> and 
P n the Legendre function of the same order having (#',</>') as its 

/* A 

axis, then \S n P n da) = - -- r$ n ', where the integration is over 
J 2n + 1 

the whole surface of the unit sphere and S n f is the same function 
ofe',<f>'thatS n isofe,<f>. 



6-45] SURFACE HARMONICS 139 

Take the axis of P n as the axis from which 6 is measured 
so that ft = cos 6. Then 

do> = sin 6d0d(f> = 



and 



(s n P n da> = f r 
J J-iJo 



where ^4 p , JS p are constants. 

/*2?r f27r 

But cosp(f>d<f> = 0, suip(f>d(f> = and J rt = P n , so that . 

Jo Jo 



where A Q is constant. 

But from 6*42 (4), when /* = 1 , X n p vanishes save when p = 0. 
Hence, when /n= 1, 







But, when /x=l, -P n =l, so that the constant A is the 
value of S n at the point on the sphere for which 6 is zero, 
i.e. on the axis (0' 9 <f>'). Hence J. = S n ' and 



6*45. Expansion of a function in surface harmonics. A 

polynomial in cos0, sin 6 cos </>, sin sin can be expressed in 
terms of surface harmonics. For this purpose we use the 
following theorem : // H n is a solid harmonic of degree n, then 

V 2 (rff n ) - m (2ri+ m + I) r m ~ 2 H n , 
where r 2 = x 2 + y 2 -f- z 2 . 

r) ?)TJ 

We have ^- (r m H n ) mr m ~ 2 xH,, + r w -^ 

C7X 3# 

and 

^-s (r m H n ) = m(m-2) r m ~ 
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By adding the similar equations in y and z, we get 



[6'45- 



But since H n is homogeneous and of degree n 



dy 



therefore 



(1). 



Now let/ n (x, y, z) be a homogeneous polynomial of degree n 
n x, y, z, and suppose that 

t + ......... (2), 



rhere H n , H n _ 2 , . . . are solid harmonics of the degrees in- 
dicated by the suffixes and the series ends with r n ~^H^ or 
r n H Q according as n is odd or even. By using the result (1) we 
can shew how to determine the H's. 

Thus by applying (1) successively, we get 



or 



(n odd) 
\n even) 



...(3). 



From the last of equations (3) we can find H 1 or H Q according 
as n is odd or even, then from the preceding equation we can 
find H 3 or T 2 and so on until the H's up to H n _ 2 are determined 
and then H n is given by (2). 

By changing into polar co-ordinates, writing rcos0, 
rsin#cos</>, r sin 6 sin </ for x, y, z, and then dividing by r n , 
it follows that any polynomial in cos#, sin 9 cos </>, sin 9 sin <f> 
can in this way be expressed in a series of surface harmonics. 
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6*5. Surface density in terms of surface harmonics. If 
as in 6'2 we find the potential at P due to a number of particles 
situated on the sphere of radius a, we get results of the form 



CO QU 

and K = S TI^> r>a (2), 

o r n+1 

where U n denotes the sum of a finite number of surface 
harmonics (one for each particle) and therefore itself a surface 
harmonic. 

We infer that the same result would be true for any arbitrar j 
distribution of surface density on the sphere, since everj 
particle of the matter gives rise to potential functions of the 
form stated. 

Conversely, we can assume (1) and (2) to give the potential 
of a certain distribution of matter and proceed to find it. By 
hypothesis U n is harmonic, so that V 2 1^ = and V 2 P^ = 0; the 
matter therefore resides on the surface of the sphere, and its 
surface density a is given by 

4770* = I -TT ^ 






It follows that if we accept the physical argument that an 
arbitrary distribution of surface density on the surface of 
a sphere produces the same kind of field of potential as an 
aggregate of particles distributed over the sphere, viz. the 
potential given by (1) and (2), then the arbitrary surface 
density is expressible in the form (3). This implies that an 
arbitrary function of the two variables 9, <f> on the surface of 
a sphere is expressible in a series of surface harmonics. The 
validity of this development in series has been discussed by 
many writers; a bibliography is given in Encyclopedie des 
Sciences Mathematiques. * 

* Tome n, vol. v. n. 176. F. also H. Poincar^, Figures d'equilibre, 
p. 52 [1902]. 
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6*51 . Reverting to the case in which the field is due to a single 
unit particle at a point M on the sphere r = a (6*2), it is of interest 
to verify that the formula 65 (3) for the surface density gives zero 
surface density except at the point M. 

Taking the values of V from 6-2 (3) and (4), the formula 6-5 (3) is 
easily seen to give 



and it is convenient to take for P n the formula 6-21 (2) putting r = a, 
so that 



r* 

= 

Lo 



vhoref, g, h are the co-ordinates of M. 
This is equivalent to 



2n( 



\"i"| 
) 1-La 



Therefore 




_ 

PM 3 

and this is zero except at M. 

6-6. Potential of a thin spherical shell. In accordance 
with 6*5 we assume that any distribution of surface density 
on the surface of a sphere may be represented by a series of 
surface harmonics. We may therefore represent the surface 
density at any point M (a, 6' ', <//) by 

............... (1). 
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If P is the point (r, 0, <f>), the potential at P due to an element 
adS of the shell is (as in 6*2 and 6*29) given by 

Tr adS I 



r <a 

a 



art& n n 

. I . Z> _l_ P _l_ I P I 

And the potential at P due to the whole sphere is got by 
substituting the series (1) for a and integrating (2) over the 
sphere, and by writing a 2 da) for dS and using the theorem 
6-44, we find that 



r <a 
\ v za ~ ba* ~ zn+ia" " } 

and 

la la 2 1 a n ~ , 

o T OT juTl ~f~ 1 T 

(3), 

where $ ', /S/, /S 2 ', etc. denote the values of the surface har- 
monics S QJ S l9 $ 2 , etc. at the point on the sphere where the 
radius through P intersects it. 

It is easy to verify that the formula of 3*7 



leads to the expression (1) for the surface density. 

6*61* Examples, (i) The density of a thin spherical shell of radius 
a is given by a = Xz 2 , where z denotes distance from a diametral plane. 
Find the potential. t 

We must first express the density in terms of harmonics. We have 



........................ (1), 

since P =l and P 2 = (3cos 2 1). We can now either write down 
the potential from 6*6 (3), or we may assume a series of harmonics 
for V and find the coefficients by using 6-6 (4). Thus, let 



r<a\ 

<*> 

and F 2 = c ^P 4-c 8 - 3 P 2 , r>a| 
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Then substitute for V l , V 2 and a in the relation 



and we get P + 2 P 2 = 7rAa 2 (P + 2P 2 ) ............... (3). 

Since this relation must hold good at all points on the sphere, 
i.e. for all values of 0, the coefficients of P and P 2 must be the same 
on both sides, so that 

c = ^TrAa 8 and c 2 = f J T 



and V l = kir\a 3 ( 1 + s _ P 2 ), r<a\ 

.(4). 



We remark that in adopting this latter method we only insert in 
the series (2) terms to correspond to the harmonics that occur in (1). 
If any additional terms were inserted, the relation (3) would serve to 
shew that their coefficients were all zero. 

Further, the total mass M of the shell is given by 

M - \ " Aa 2 cos 2 6 . 27ra 2 sin 6d6 = TrAa 4 , 
Jo 

arid since 2r 2 P 2 = r 2 ( 3 cos 2 6 - 1 ) = 2z 2 - x 2 - y*, 

therefore the expressions for the potential may be written 



V, = + **(W-*-in, r>aj 

(ii) The density of a solid sphere of radius a is given by p \xyz. 
Find the potential at an external point. 

Consider a thin shell of radius a' and thickness da'. Wo may regard 
it as a layer of surface density 

a pda' = \xyzda' . .1 ( 1 ). 

Now xyz is a harmonic function of the third degree, so that by analogy 
from 66 this and its associated harmonic xyz/r 1 [6*13] are the only 
harmonics which can occur in the expressions for the potential. We 
therefore assume that for the potentials of the shell 

V l = Cxyz, r<a' 
and F 2 = Ca' 1 - , r > a', 

arranging the constants so that V l = F 2 when r = a', F x is finite at the 
origin and F 2 -> as r ~> oo. 
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Then substituting in 

dV, 8V t 
-' 



we get 

so that C = iTrAa' da'. 

Hence the external potential of this shell is f TrAa' 8 da' . -~ . And by 
integrating with respect to a' between limits and a, we get for the 
potential of the sphere ^TrAa 9 -~- . 






(iii) A mass M is distributed over a spherical surface of radius a so 
that the surface density at any point is proportional to t/ 2 z a -f z z x 2 + x 2 y 2 , 
where the origin is at the centre of the sphere. Find the potential inside 
and outside the sphere. 

It is necessary in the first place to express the surface density in 
a series of harmonic functions. Proceeding as in 6*45, we put 



...... (1), 

where H 4 , H z , H are harmonics of degree 4, 2, 0. 

Then V 2 /^2.7// 2 + 4.5/- 2 // ..................... (2) 

and V 4 /=4.5.2.3H ........................ (3). 

But, by differentiation, 

V 2 /=4(z 2 + 2/ 2 + z 2 ) 
and V 4 /= 24. 

Therefore, from (3), (2) and (1), 

# =J, # 2 = 
and // 4 = y*z 2 + z*x* + x*y* - Jr 4 . 

We may therefore express the surface density in the form 
a = A (?/ 2 z 2 + z*x* + x*y* - Jr 4 ) + J Aa 4 , 

and assume expressions for the potential containing the same har- 
monics, viz. 

F! = A f B (y*z 2 + z 2 ^ 2 + x*y* - r 4 ), r < a 



and V^A- + B (y*z* + zV + x*y 2 - Jr*), r > a. 



Then, by substituting in the relation 
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A. QB 

we get + *z* z*x* # - 



so that A faXa 6 and JS $ 

But M = I crd, integrated over the sphere, 

= J / Xa*dS, since the integral of a harmonic is zero, 



Hence it follows that 

V l = M g + Q J 5 (2/ 2 * 2 -f ^ 2 ^ 2 4- ^ 2 2/ 2 - 
and F 2 = M JA + 1 ^ (2/2^2 + ^2^2 + ^2^2 _ i 

6-7. Nearly spherical bodies. Suppose that the strata of 
equal density in a body are nearly spherical, so that the 
equation of a surface of equal density is 



where the coefficients C n are so small that their squares and 
products can be neglected. 

The volume contained by such a surface is Jjr 3 dco [1*42], 
which 



= |7ra 3 , since (s n da) = 0, [6-41], 
= vol. of a sphere of radius a. 



We note as special cases: 
(i) The surface r = a 4- eP^ , 
or r = a + e cos 6. 

When e is small enough, this is a sphere of 
radius a with its centre C at a small distance e 
from the origin. 

(ii) The surface r = a + eP 2 
or r 2 = a 2 4- 2a JP 2 . 
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To the first order in powers of , this is the same as 



r 2 
~ 

In Cartesian co-ordinates this can be reduced to 



so that the surface is a spheroid of semi-axes a Je, a-f e. 

We also notice that the normal to such a surface makes with the 
radius vector r an angle x which is at most of order c. So that for 
differentiation along the normal we have 

dv_ dv a_F dv 

dn~ dr ^rdd r sin 0d</>* 

where A = cos x = 1 to the first order, while /*, v are of higher order than 
the first in . Consequently, to the first order in e, we have 



6*71. Potential of a nearly spherical body. Consider a 
body stratified as in 6^7. The surface 



(1) 



may be regarded as the boundary of a solid sphere of radius a 
and unit density surrounded by a layer of surface density 

c n s n . 

The potentials inside and outside the sphere are, from 6-6, 



r<a 



TT t _ n , 

and U 2 = - + ^--~S n ', r>a 



If we regard a as variable, the differentials of U^ , U 2 will 
give the potentials inside and outside the shell whose para- 
meters are a and a + da. The density p of the corresponding 
stratum of the given body is by hypothesis a function of a, 
and if c, c' are the Tr alues of a on the inner and outer boundaries 
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of the body, the potentials inside and outside the shell are 
given by 



da 

,. 477 r* ( 2 d ^ nn ,j 
and K = pla 2 +- r -I l -~^~ - \da\ 

2 r da 



The potential in the substance of the body, on the surface 
of parameter 6 say, will be found by adding together the first 
integral between limits 6 and c' and the second taken from 
c to 6. 

6*72. Example. The form of a nearly spherical uniform gravitating 
solid of mass M is given by 



where I, m, n are the direction cosines of r. Prove that neglecting 2 the 
exterior potential is 

3 



If the solid is surrounded by a thin layer of water just sufficient every- 
where to cover its surface^ shew that t neglecting the self -gravitation of the 

water, the maximum depth of the water is o~, 7 [C. 1931] 



Since e is small we may regard the solid as formed of a sphere of 
radius a and uniform density p surrounded by a thin layer of surface 
density & = aplmn; and by [6*7] ^7ra*p = M 9 the whole mass. 

Also, since xyz or rHmn is a solution of Laplace's equation, therefore 
[6'13J xyzjr 1 or Imn/r* is also a solution. 

The potentials mside and outside the sphere due to the surface 
layer may be taken to bo of the form 

__ r 3 __ a 4 

Q& |4 

so that F! is finite at the centre, F 2 vanishes at infinity and F x = F 2 
at the surface. Then from 3*7 and 6-7, 

(* 2 ) = 47ryo-=47ryae/)Zmn. 
or or / r=KO 

Therefore 7 A = 47rya 2 ep = 3yMc/a, 

and the external potential due to the whole body is 

(i). 
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The free surface of the water must be an equipotential surface 
because in equilibrium the surfaces of equal pressure are at right angles 
to the resultant force. 

Hence the free surface is given by an equation 



(2), 



where C is a constant to be determined from the fact that the water 
only just covers the solid. The surface of the water will touch the 
solid on its protuberances, i.e. where the surface r = a(l + lmn) has 
its maximum radii; and the water will be deepest over the hollows 
on the surface ofihe solid, i.e. where r = a (1-f clmn) has its minima. 
It is easy to shew that Iran has maxima and minima values subject 

to Z 2 -hw 2 + n 2 =l, where Z=TO w= 775, n= -TO maxima 

y3 y3 y 6 

when the signs are chosen to make Imn positive, and minima when 
negative. 

Hence C is to be found by substituting Q rs for Imn and a 1 1 -f tjrj* ) 
for r; so that 



21 V3/' 

and the free surface is given by 
3 



A minimum radius of the body is obtained by putting Imn ~ ^- , , 

o V o 

so that the minimum radius is r = a ( 1 Q~TQ) > anc * the radius of the 
free surface in this position is given by 



r 
To the first power o$ e, this gives 



21-V/3 



, 



N 

21\/3/' 



r=a 1 + 



"21V3J- 
So the greatest depth of the water 

e ^ afl <-} 

e\ I U- I 1 ij . o I 

3y3/ 
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6-8. Clairaut's Theorem.* Variation of gravity. As the 

result of the measurement of meridian arcs there is reason 
to suppose that the earth is a slightly oblate spheroid. If 
a, c are the major and minor semi-axes of such a spheroid 
and c = a(I~- e), e is called the ellipticity of the spheroid. 
Neglecting the square of the ellipticity the equation of the 

surface is 9,0 > 

x* + y 2 z 2 = 

" a 2 a 2 (l-2e)~ 
or x 

or 
or r = a(l-cos 2 0) .................. (1). 



For the earth, e is about 1/300. 

The earth rotates with a small angular velocity o> = 277/24 . 60 2 . 
Bodies regarded as at rest relative to the earth are therefore 
subjected to a 'centrifugal force' as well as to the earth's 
attraction; so that if V denotes the potential of the earth's 
mass and X, F, Z denote the whole force on a unit particle 
on the surface, 



As by far the greater part of the earth's surface is fluid, 
it is reasonable to argue that the surface must be one of con- 
stant pressure, and from hydrostatical considerations this 
requires that Xdx+7dy + Zdz = .................. (3), 

so that V + $a>*(x* + y*) = C f .................. (4) 

over the surface, where C is a constant. 

If E denotes the mass of the earth, its potential at a great 
distance is E/r, as a first approximation. We have to correct 
this because of the term ecos 2 in (1), and it will suffice for 
our purpose to add a solution of Laplace's equation containing 
cos 2 0, i.e. a term AP 2 /r 3 , where A is a small constant of 
order e. 

* Alexis Claude Clairaut (1713-1765). French mathematician. 
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V A P 
Thuswetake F = -+4* ........................ (5) 

T T* 

and determine A so that (4) is satisfied by (5) when r is given 
by (1); i.e. 



...... (6). 

r = a is a sufficient approximation in the second and third 
terms since A and o> 2 are of order e. 

Since (6) is to be satisfied for all values of 0, by equating 
to zero the coefficient of cos 2 0, we get 



and the terms independent of cos 2 serve to determine C. 



Hence V = ~+aj^- (cos 2 0- J) ...... (8). m 

r r \ a J 

The radius from the earth's centre to the point on its surface 
makes an angle with the axis Oz. The force of gravity g 
at the earth's surface is the resultant of the attraction and the 
centrifugal force and if it makes an angle v with the radius 
vector the radial force is g cos v. But v is so small that its 
square is negligible, so we have 

g = - =- ( F + I o> 2 r 2 sin 2 0) at the surface, 

or 



E So 3 



_ 6E \ 



where r = a ( 1 e cos 2 0) and o> 2 is of order . 
This gives 

E I eE\ 

^ = - 2 (l+)^|o> 2 a + U^at cos 2 (9). 

ct \ d j 

If (? denotes the value of gr at the equator JTT, and m 
denotes the ratio of the centrifugal force to gravity at the 
equator so that co 2 a = mG, we have 

(10), 
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so that E = (l + fm-<r)(?a 2 .................. (11) 

and (9) becomes g = G{l + (f m - e) cos 2 0} ............... (12). 

We observe that if (? is found by means of a pendulum, 
(11) determines the mass of the earth. 

The numerical value of m or a) 2 a/G is about 1/289, so that 
it is of the same order as , and in the small term in (8) we may 
substitute aj*a = mG = mE/a 2 , so that the potential of the 
earth at an external point is given by 

pi 



(13). 



6-81. The Moon's attraction. The forces exerted by the 
moon on the earth are equal and opposite to the forces exerted 
by the earth on the moon, and are approximately the same 
as if the moon's mass M were collected at its centre of gravity. 
If r denotes the distance between the centres of the earth and 
moon, and 6 the moon's north-polar distance, the earth's 
attraction at the centre of the moon is composed of a radial 

'dV IdV 

attraction P = ^ - and a force = - - tending to increase 6. 

dr rd6 to 

Hence from 6-8 (13) 



Q = - 2 (lm - e) f sin cos 
r* 

The forces exerted on the moon are M'P radially and M Q 
through the moon's centre. The force system exerted on the 
earth by the moon consists of equal and opposite forces 
MP, MQ acting through the earth's centre and a couple 
of moment M Qr tending to cause rotation about an equatorial 
diameter, and since e \m is positive the sense of the couple 
is such as to decrease 6. The moment of this couple is 



2 (e -4m) -^3- sin cos ............... (2). 
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6*82. We can obtain some information by comparing the 
results of 6-8 and 6-81 with MacCullagh's expression for the 
potential in 4'6 and the couple deduced therefrom in 4*62. We 
notice that the potential as given by 6*8 (13) is independent 
of the longitude, therefore when the attracting body is the 
earth MacCullagh's expression 

E A+B + C-3I 

r + 2r 3 + "' 

should be independent of the longitude. This is satisfied by 
taking a principal axis of the earth as the axis of rotation and 
writing A = B, so that / = A sin 2 + C cos 2 0, and 



A comparison with 6-8 (13) now shews that 
C-A 



Ea* 



-(e- Jm). 



The same result follows from comparing the couples of 
4-62 and 6-81. 

6-9. General solution of Laplace's equation. It has been 
shewn by Whittaker* that there is a general solution of 
Laplace's equation of the form 



PTT 

== 

J-7T 



provided that / is such a function that differentiations with 
regard to x, y, z un4er the sign of integration are permissible; 
and it may be shewn that V can be expressed as a series of 
expressions of the types 

(z + ix cos u + iy sin u) n cos mu du, 

r 

(z + ix cos u 4- iy sin u) n sin mudu, 

r 

* Math. Ann. LVII (1902), p. 333. See also Whittaker and Watson, 
Modern Analysis, ck. XVHI. 
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where n and m are integers such that O^ra^n; also that by 
changing to polar co-ordinates these may be reduced to solu- 
tions of the forms X n p cosp<, X n p sinp< in the notation of 
6-42, 6-43. 

EXAMPLES 

1. Find, in terms of /u, P 4 (/x), P 5 (^), P 6 (/A). 

2. Prove that 

(x) dx = H 



[M. T. 1912] 
3. Prove that '""'" ' ^ ' * 



"2n+l* 
4. Prove that P w P n can be expressed in the form 

where m>n and the c's are constants. [C. 1908] 



e w n , n . ,, . 

.5. Express -~? , x , n and , 2 m the form 
ctx dx dx 

^-f^ 1 P 1 + ^L 2 P 2 + ^3P 3 + . 
n /f?p 

(-/* 

J -i\ /A 

, 
dp, 



6. Prove that 



7. Prove that P (1 -w 2 )^ 
J _! x ^ ' 



and 



8. Shew that, if m and n are integers, the value of 



r> 

is either 0, 2 or ^-^ . [M. T. 1908] 

* 1 ^" 



9. Prove that, if /(/x) is expressible in the form S a n P n (/x), then 





o 

10. Prove that I P m (u) P n (u) d/z is equal to zero if m and n are 

Jo 

unequal and both even or both odd integers, while if m is even and n 
is odd, its value is 

/ i)i<+n+i) 1.3... (m1) 3. 5... ft 



2.4...m 2.4 ..:(n- 1)' 
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11. Explain why, if F (cos 0) satisfies the equation 



so also does F (A cos 6 -f- /z sin 6 cos <j> -f v sin 6 sin </>), 

where A 2 + ft a +v a =l. [M. T. 1921] 

12. Shew that 

_ 1) ^ + 1 i ; (2n + 3) . [C.1905] 

13. If V n is a homogeneous function of x, y, z, of degree n satisfying 
= 0,then 

n satisfies V 2 V 2 F = 0. [C. 1891] 



14. Find the potentials of a surface distribution on a sphere in 
which the density at any point is proportional to the distance from 
a given piano. [C. 1904] 

15. A spherical surface of radius a is coated with matter of surface 
density proportional at each point to the cosine of the angle the radius 
vector from the centre to the point makes with the positive directign 
of the axis of x. Prove that the potential due to this matter at any 
external point is $7Toa 3 x/r 3 , and at any internal point \TTCTX; where the 
origin is. at the centre and a are the densities respectively at the 
points where the axis ofx cuts the surface. [C. 1890] 

16. The surface density on a sphere of radius a is key. Find the 
potential inside and outside the sphere. 

17. The surface of a sphere of radius a whose centre is at the origin 
is covered by a thin layer of attracting matter, the surface density at 
any point P being proportional to sin 2 POZ. Find the potential at 
any point inside or outside the sphere. 

Prove that inside the sphere the lines of force are given by 
y/x constant ; zx z = constant. [London Univ. 1925] 

18. Find the surface arid volume distribution of matter for which 



V = kr p (3 cos 2 e I ) when r ^ a, 

and V = 1ca,P+*r-*(3cos 2 e-I), when r>a, 

k, p and a being constants. [London Univ. 1936] 

19. If the mass, per unit area of a thin heterogeneous stratum of 
attracting matter placed on a sphere of radius a is pz 3 /a 3 , prove that 
the potential at a point within the sphere is 



where r 2 = # 2 4- 1/ 2 -f z 2 , and find the potential at a point without the 
sphere. [M.T. 1903] 
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20. Prove that Laplace's operator V 2 is invariant for any change of 
rectangular axes, and that a solid harmonic (or a surface harmonic) 
of degree n transforms into a solid harmonic (or a surface harmonic) 
of degree n in the new co-ordinates. 

21. Prove that, if x be real and positive, 

(1) f 1 p n(^ d P 2* 



2* 2* 



8 >* (2n+l)(2n-] 

[C. 1908] 

22. If P n is Legendre's coefficient of order n, prove that 

_ 2 (2n-4r+l){r(2n-2r+l)-2}P n _ ar 

when n is even, and that when n is odd the summation is from 1 to 
4(n-l). [M. T. 1904] 

23. Express xyz (x 3 + 2/ 3 + z 3 ) in the form 



where F 6 , F 4 , F 2 , V are solid harmonics of the orders indicated by 
the suffixes. [London Univ. 1931] 

24. Express in the form F 6 -f F 4 + F 2 , where F n denotes a spherical 
solid harmonic of degree n, a function which satisfies Laplace's equation 
within a sphere r = a, is finite at the origin, and is equal to x 3 y z z on 
the sphere r~a. [London Univ. 1926] 

25. Obtain a rational integral harmonic which has the value Az* 
at points on the surface of a unit sphere with its centre at the origin. 

26. A thin layer of matter is placed on a sphere of radius a whose 
centre is O, the density at any point P being proportional to cos 4 A OP, 
where OA is a fixed radius. Determine the potential at all points 
inside and outside the sphere. [London Univ. 1927] 

27. Matter is distributed in a thin layer over the surface of a sphere 
of radius a and centre at the origin, the density being proportional 
to (x + y) 2 and the whole mass being m. Shew that the potential at an 
external point is given by 



and find the value at an internal point. [7<ondon Univ. 1938] 
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28. If the surface density at the point (x, y, z) on a sphere of radius a 
is x 2 y xy 2 , prove that the potential is 



at an internal point, and 



at an external point. [London Univ.] 

29. The sphere x* + y 2 + z 2 = a 2 is covered by a thin layer of attracting 
matter, the total amount of matter being M and the surface density 
at any point (X, Y, Z) of the sphere being proportional to 



which may be written as 



Find the potential of the layer at any external point. 

Shew that the attraction on a small mass m, placed outside the 
sphere on the line x y = z at a distance r from the origin, is 

[London Univ. 1925] 

30. If V n is a rational integral harmonic function of degree n and 
r 2 = x 2 + y 2 -f z 2 , and U is a homogeneous function of x, y, z of degree 2n, 

n 
shew (by expressing U in the form S r 2n - 2p V 2j) , or otherwise) that 



the integral being taken over the surface of the sphere r = a. 

[M. T. 1939] 

3 1 . The density of a solid sphere of radius a, referred to axes through 
its centre, is x 2 yz. Shew that the potential of the sphere at an external 
point (x, y, z) is given by 



where H 2 = yz, H^ = 6x z yz - </ s z - yz*. [M. T. 1908] 

32. The potential over a sphere of radius a is equal to Axyz 2 , and 
is due to some surface distribution a. Find the potential at an internal 
point, at an external point and the value of a. [C. 1913] 

33. Find the potential in all space when the potential at any point 
of the surface of a -phere is y 2 z 2 . [C. 1906] 
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34. The density of a solid sphere of radius a at the point (x, y, z) 
is kz 3 , where A; is a constant. Prove that the potential at any point 
inside the sphere is 

[| (5z 2 - 3r 2 ) (9a 2 - 7r 2 ) + 7a 4 - 3r*], 

where r 2 = x 2 + 2/ 2 + z 2 . [London Univ. 1926] 

35. The density of a solid sphere of radius a is po + Pi^2 where 
p , p are constants and $ 2 is a spherical surface harmonic of the second 
degree with origin at the centre of the sphere. Prove that the potentials 
at internal and external points are 

(3 a -r') + f wy^r* (log" 



and $ Trypo ~~ + ./^ try Pl ^ S 2 . 

[London Univ. 1931] 

36. The density of a solid sphere of radius a is proportional to # 2 z 2 , 
the centre being the origin. Prove that the potential at an external 
point is 

ML a 2 -22/ 2 + z 2 a 2 7a 4 5 , 2 2 .a 4 
y r {1 + - ^9 - ^-9^-99 (a "^ +2) ^ 

[London Univ. 1926] 

37. The matter* of a heterogeneous solid sphere is arranged 
symmetrically about a diameter. It is found that the normal com- 
ponent of attraction due to the sphere at a point P on its surface where 
the radius makes angle with the diameter is A + B cos 4 6. Prove 
that the tangential component of attraction at P is 

s 2 ^). [C. 1907] 



38. A given mass is distributed in a thin layer over a spherical 
surface, symmetrically with respect to a diameter AB and to the 
corresponding diametral plane, so as to produce inside a field, such 
that the force at any point towards AB is proportional to the distance 
from it. Shew that the force parallel to AB is proportional to the dis- 
tance from the corresponding diametral plane. Shew also that if the 
surface density in the cjia^metral piano is zero, that at A is three times 
what it would be if thejin&ijss were uniformly distributed. [C. 1901] 

39. Shew that the potential of a uniform circular ring of radius a 
and mass M is 



( * C 2 

TT J 



where the origin is at the centre and the axis of z along the axis of 
the ring. 
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40. Verify, by using the formulae of 628 for P n , that the expressions 
Fj , F 2 of 6-32, for the potential of a circular disc, are equal when r = a, 
by proving them both equal to 



fir 

[!+{/!+ V(^ 2 ~l) cos <} 



2m 

41. Prove that, if the potential of a body symmetrical about an 
axis at a point on the axis at a distance z' from an origin is $(z'}, 
then the potential at any point outside the body is 



where r is the distance of the point from the origin, and z the projection 
of r on the axis. [M. T. 1887] 

42. Matter of mass M is distributed on the surface of a sphere whose 
centre is O and radius a, so that the density at any point is proportional 
to the square of its distance from an external point (7, where O(7 = 6. 
Prove that the potential at an external point P is 



where OP = r and z = r cos P0<7. [C. 1897] 

43. The density of a solid sphere at any point P varies inversely 
as OP, where O is an external point distant / from the centre. Shew 
that the potential at any external point is 



3M /a 2 \ n _ _ P n (cos 6) 
r n\fr (2 



r Vr 

where M is the mass of the sphere, a its radius, and r, polar co- 
ordinates whose origin is the centre and whose axis is the line from 
the centre to O. [M. T. 1909] 

44. The density of a thin spherical shell, centre (7, radius a and 
mass M, varies inversely as the distance from an external point 0. 
Shew that the potential at any external point P is 



IM (<#l 
2 a Jo 



and at any internal point P is 




2 a o. 
where r = CP, /= GO and** is the.Wgfe' ff& [M. T. 1895] 

45. A disc of radius a is of a; jnass M &ttd its surface density is 
proportional to the mth power of the distance from its centre (m + 2 > 0) . 
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Shew that the potential at points whose distance (r) from the centre 
is greater than a is 



[London Univ.] 

46. A thin circular disc of radius a is such that the mass per unit 
area at distance r from the centre is k(a 2 r 2 ). Find the attraction 
of the disc at any point on its axis. 

Henqe shew that the component attraction parallel to the axis of 
the disc at any point whose co-ordinates are r, 6, </> referred to the centre 
of the disc as origin and its axis as polar axis can be expressed in the 

form g 



ifr>a. [M. T. 1920] 

47. Prove that at a point in free space, where r>, the potential 
of a homogeneous gravitating hemisphere of radius a and density p is 



2 2 - 2 n! (n- 2) ! r*"-* 

where r is the distance of the point from the centre of the base and the 
axis of harmonics is the axis of symmetry of the hemisphere drawn 
from the centre into the interior. [London Univ. 1931] 

48. Shew that the potential of an oblate spheroid of density p 
and major and minor semi-axes a, c at large distances is approximately 



4 "| 
* J" 



1Q9Q , 
1929] 



49. A uniform plane annulus of mass M is bounded by circles of 
radii a, 6, large compared with the dimensions of a spheroid whose 
centre of gravity coincides with the centres of the circles. Shew that, 
if the normal to the plane of the? annulus makes an angle 6 with the 
axis of the spheroid, they exert on one another a couple of magnitude 

ZyM(A-C) . Q 
-JL\ -- ' sm cos Q 

ab(a + b) 

tending to increase 6, where A, A, C are the principal moments of 
inertia of the spheroid f and y is the gravitation constant. 

[M. T. 1931] 

50. A uniform nearly spherical solid of density p has the surface 
r = a(l -t-eP 2 ) as its boundary. It is surrounded by liquid of volume 
47r(6 3 a 3 )/3 and uniform density a. Shew that, provided the solid 
is completely covered with liquid, the equation of the free surface is 
r = 6 ( 1 -f T7JP 2 ), where 

3(p g)a 5 rlvr r 
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Put a 2 tan 2 d = u, and on reduction we get 

du 



f 

\ ~^ 



and similar expressions for B and C. 

7-32. Spheroids. We see from 73 (3) that 

*- (A - 



and 



da) 



-* + ^*+ c 
Two of the axes of a spheroid are equal, and if a = b then 



and ( 1 ) become 2 A + C = 4-Tr, 



a 2 + c 2 



o o 



where w = cos ^. 

There are two cases: 

(i) Oblate spheroid, a>c. 



where e is the eccentricity of the generating ellipse. 
(ii) Prolate spheroid, c> a. 



.(1). 



-o 

* ^ (2). 



where a 2 = c 2 (l e 2 ), so that e is the eccentricity of the generating 
ellipse. 

A and C can then be found in either case, since 2 A + C = 47r. 
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7*33. Attraction of an oblate spheroid of small ellipticity. 

When the ellipticity e is small, since c = a(l e) and C 2 = a 2 (l e 2 ), 
therefore e 2 = 2e approximately. 
Hence from 7*32 (i), 



-A Jl _ 1 __^_ ,]__?! __ I 

~ ^f 31-e 2 " i "5(l-e 2 ) 2 "'J' 
Whence, since 2 A + C = 4-rr, 

n ^2 i />4 i 
ryg d^ ? _ v 

06 - 4 H3r^~5(T^yj 






But 2,4 + (7= 47r, so that A = Jw (1 - f e). 
Hence in this case 



These results are of importance in connection with figures of equi- 
librium of rotating liquid. 

7*34. Attractions of a nearly spherical ellipsoid. Let 

b = a(l 17) and c a (I e). Neglecting higher powers of 77 and 
than the first, we may put A=A + Fy + F'e. 

If the ellipsoid were turned through a right angle about Ox, b and c 
would be interchanged, so we must have FF'. Also an oblate 
spheroid for which b ~ a and 7? = is a special case, for which, from 733, 



therefore 
Hence 

Similarly 



and C = Jw(l 

Therefore X= - jwp (1 - frj - f )/"j 

Y= _Vp(l + |ij-| e )jr ..................... (1). 
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We may obtain similar results in a symmetrical form, if we take k 
as the mean radius J (a -f 6 + c), then 



-6-c\ 
) 
a ) 



22a-6-c\ 

z 
5 



t-k\} 



Similarly, ^ _ Vll , A - 

\ O K 

/ 

and 



.(2). 



5 k 

7*35. Potential of a homogeneous solid ellipsoid at an 
internal point. From 7'21 (3) it follows that, since a solid 
ellipsoid can be divided into similar homoeoids, the innermost 
being of zero dimensions, by putting m = 0, the potential at the 
centre of a homogeneous solid ellipsoid is 

du 



To find the potential at an internal point (#, ?/, z), we make 
use of the fact that the attraction components are given by 

X,Y,Z=-Apx,-Bpy,-Cpz (7-3) ...... (2). 

Hence, if V denotes the potential, 



^- ^~ -^- 

ox oy oz 



= p(Axdx + By dy -\~Czdz). 
Therefore V = lp (D - Ax* - By* - Cz*) ............ (3), 

where the constant of integration %pD must be the potential 
at the centre. Hence by substituting for D from (1) above, 
and for A, B, G from 7-31, we get 



F = 



f 
J ^ 



_ ^_ 

( X 2 
l -tt- 
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It follows from (3) that the equipotential surfaces inside 
the ellipsoid are ellipsoids 

= const ................ (5), 



similar to one another though not similar to the boundary 
of the solid. 

7*4. Ivory's Theorem.* Attraction of a homogeneous 
solid ellipsoid at an external point. If x, y, z and x 1 ', y' ', z' 
are the co-ordinates of points P and P' one on each of two 
confocal ellipsoids 

aja/a a + y a /& a + * a / c8:=1 .................. (1), 

and x*/a'* + y*/b'* + z*/c'*=l ............... (2), 

then P and P' are called corresponding points if 

x/a = x'/a', y/b = y'/b', z/c = z'/c' ......... (3). 

Let US be an elementary strip of the first ellipsoid of cross- 
section dydz, parallel to the 
#-axis; and let R S' be the 
corresponding strip of the , 
second ellipsoid of cross- 
section dy'dz', so that 

dydzjdy'dz' bc/b'c'. 

Then if /' (r) denotes the 
law of force at distance r, 
and p the density of either 
ellipsoid, the component 
parallel to the #-axis of the 
attraction at P' due to the strip BS 



P' 




= -pdydz /' (r) cos P'QSdx, 
when Q is the position of an element dxdydz. This 

= -pdydz J/' (r) ( -^j dx = P dydz{f(P'S) -f(P'B)}. 
* Sir James Ivory (1765-1842). Scottish mathematician. 
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In the same way the as-component of the attraction at P 
due to the strip R'S' 

= P dy'dz'{f(PS')-f(PR')}. 

But it is easy to deduce from (3) that for any pairs of 
corresponding points PS' = P'S and PR' = P'R. 

Hence the ratio of the ^-component of the attraction of 
R8 at P' to the similar component of R' S' at P is 

dydz _ be ( . 

dy'dz'~~b'c' ..................... * '" 

If in like manner we take all such strips as R8 of the first 
ellipsoid and all such strips as R'S' of the second, we get 

X:X' = bc:b'c f ..................... (5), 

where X is the ^-component of the attraction of the first 
ellipsoid at P' and X' is the ^-component of the attraction 
of the second ellipsoid at P. 

This result is Ivory's theorem and is true for any law* of 
force. 

Now, from 7'3, X f = - A'px, where Pis the point (x, y, z) and 
A' is the same function of a', 6', c' that A is of a, 6, c. 

Therefore 

be A , abc A , x abc ., , 



where x', y', z' are the co-ordinates of P'. 

Similar relations hold for the y and z components of attrac- 
tion, so that the attraction at P', (x', y', z'), of the ellipsoid 
(a, b, c) has components 

* abc ' 



wt 

abc 



where a', 6', c' are the semi-axes of the confocal ellipsoid which 
passes through P', and A', B f , C' are the same functions of 
a', 6', c' that A, 9 B> C are of a, 6, c. 
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7-41. MacLaurin's Theorem.* The attractions of two con- 
focal ellipsoids at a point external to both are proportional to 
their masses and in the same direction. 

In 7-22 this was proved to be true of confocal thin homoeoids; 
that it is true for homogeneous solid confocal ellipsoids follows 
as a direct corollary from 7-4 (6), if we regard P' as a fixed 
point and consider the attractions at P r of two different 
confocaj ellipsoids, of semi-axes a, 6, c and density p and of 
semi-axes a l ,b l , c and density p t ; the #, t/, and z components 
are all in the ratio pabc : pi^b^ and are therefore propor- 
tional to the masses and have resultants in the same direction* 

7*42* Potential of a solid homogeneous ellipsoid at an 
external point. Let a, 6, c be the semi-axes of the given 
ellipsoid and x, y, z the co-ordinates of the external point P. 
Let a', &', c' be the semi-axes of a confocal ellipsoid which 
Y>asses through P, and F, V the potentials at P of the first 
and second ellipsoids regarded as of the same density. 

Then by MacLaurin's theorem (7'41) 

V:V' = abc:a'b'c'. 
But, from 7-35 (4), 



Therefore 



(ft n 9 \ 

x* y* z* \ i 

I _ .JL ----- _ d u 
a' 2 + u b' 2 + u c' 2 + u} 



( 
l -rf 



but since the ellipsoids are confocal, we have 

where A is the largest root of the cubic 

x 11 z 



* Colin MacLaurin (1698-1746). Scottish mathematician. 
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Substituting in (1) for a' 2 , 6' 2 , c' 2 , and writing v for 
we get 

f 1 

F = ^a&c J ^ ^ (a2 + |;)(62 



where the lower limit A is determined by (2) as stated. , 

7*5. Attraction of an infinitely long homogeneous solid 
elliptic cylinder. This can be deduced from 7*3 and 7-4 by 

making c -> oo. 

(i) Internal point. The components of attraction at a point 
(x, y, z) are given by 

X=-Apx, Y=- 
where, from 7-31, 



r 

J 



o 
By substituting for a 2 -f u, we get 

A A j, v dv 

A = 4:irab 



sf 

Jo 



or 



' a + 6" 

Similarly B = ~^"^ ; 

J a-h6' 



and z. ; 7 =- (1). 

a+6 a a+6 6 

(ii) External point. In like manner from 7-4 we find for 
the attraction components at an external point (x, y) 

4:7T P ab x lirpab y^ 
a'+b' m a" I ~~a'+b''b' ( >' 

where a', b' are the semi-axes of an ellipse confocal with the 
cross-section of the cylinder and passing through (x y y). 
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7'6. Equilibrium of rotating liquid. When a mass of gravi- 
tating liquid rotates uniformly about an axis, it is conceivable 
that for a certain form of the free surface the liquid particles 
may be in relative equilibrium. Since the resultant attraction 
on a particle depends on the form of the boundary surface, 
which is unknown, the problem does not admit of a complete 
solution. But assuming the liquid to be homogeneous, it can 
be shewn that certain forms are possible forms of relative 
equilibrium. 

Thus if p be the density and w the angular velocity, whether 
spheroidal or ellipsoidal forms are possible forms of relative 
equilibrium is found to depend on the numerical value of 
a) 2 /27rp.* We shall illustrate the theory by considering two 
simple cases. 

7*61. To shew that an oblate spheroid of small ellipticity 

is a possible form of relative equilibrium of a mass of uniform 
liquid rotating about an axis with a small uniform angular 
velocity. (MacLaurin's theorem.) 

The problem can be treated as a statical one if we compound 
with the force of attraction the reversed effective force, o> 2 r per 
unit mass, where r denotes distance from the axis of rotation. 

The equation for the pressure is then 

= Xdx + Ydy + Zdz + aj 2 (xdx + ydy) 
P 



or == Ap (xdx -\-ydy) Cpzdz + a> 2 (xdx + ydy), 
P 

assuming that the boundary is an oblate spheroid. 

But the pressure must be constant over the boundary, so 
that its equation must be of the form 

(a)*-Ap)(x 2 + y*)-Cpz* = const .......... (1). 

If a, a, a (1 e) are the semi-axes of the oblate spheroid and 
the ellipticity e is small, its equation is 



* Besant and Ramsey, Treatise on Hydromechanics, Part I, ch. vin, 
where a discussion of the subject with full references will bo found. 
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By comparing (1) and (2), we get 

a>*-Ap=-Cp(l-2) ............... (3). 

But for an oblate spheroid of small ellipticity, from 7-33, 
f) and C = $ir(l + c); 



IK 2 

and, by substituting these values in (3), we find that e = - , 

lo Tip 

and with this relation between and a) the condition is 
satisfied. 

7-62. Jacobi f s ellipsoid.* An ellipsoid with three unequal 
axes is also a possible form. 

For an eUipsoid the pressure equation is 



P 
where A, B, C are given by 7-31. 

The surfaces of constant pressure are therefore 

(a>*-Ap)x* + (Gj*-Bp)y 2 ~Cpz* = const ....... (1). 

And if we assume the free surface to be 



we get by comparing (1) and (2) 

a*(a>*-Ap) = b*(a>*-Bp) = -c*Cp ......... (3). 

By eliminating o>, we get 

a*b*(B-A) = (a*-b*)c*C ............... (4). 

If we put D = {(a 2 + u) (6 2 -f u) (c 2 + u)}* 9 then, from 7-31, 



r du rt , r 

7 - 5 - -=r , .B = 27ra6c 
Jo (a* + u)D J J 



and (4) becomes 

2a 2 =0 



-D 
* Carl Gustav Jacob Jacobi (1804-1851). German mathematician. 
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Hence, if a =f= b, the lengths of the axes must be such as to 
satisfy the equation 



Jo 



If a and fe are given, then (5) is an equation for determining 
c, and since the left-hand member is negative when c is small, 
and becomes positive as c increases, there must be a real value 
of c wnich satisfies the equation. 

Further, o> 2 is given from (3) in the form 



udu 



where M is the mass of liquid. So that when a 4= 6, 

udu 



foo 

5 = fJf 7-2 

Jo (a 2 



and since this expression is positive there is a real value of o>, 
thus establishing that the ellipsoid is a possible form. 
Further, from (3), 



udu 



so that for w to be real we must have c < a. Similarly c < 6, so 
that the ellipsoid must rotate about its least axis. 

7-63. Example. A rigid gravitating sphere of radius a and uniform 
density p is surrounded by a layer of gravitating liquid of volume 
$TT (6 3 a 3 ) and uniform density a and placed midway between two spheres 
of mass M whose distance apart (2f) is so great that powers of b/f above 
the third may be neglected. Prove that the boundary of the liquid is a prolate 
spheroid of small ellipticity 

45M6 3 



and determine the small angular velocity of the whole mass which, to this 
approximation, would make the boundary a sphere. [P. 1933] 

Let A, B be the centres of the two spheres of mass M and O that of 
the sphere of radius a. Let Q be a point in the liquid at distance r 
from 0, where the angle QOB0, and let OB be the axis of z. 



